m 


THE  TURBULENT  FLUXES  OF  MOMENTUM  AND 


SENSIBLE  HEAT  OVER  THE  OPEN  SEA  DURING 


MODERATE  TO  STRONG  WINDS 


WILLIAM  GEORGE  LARGE 


OCEANOGRAPHY 


THE  UNIVERSITY  OF 
BRITISH  COLUMBIA 


this  ■ 

\Aoi 


80  2  27  046 


REPORT  DOCUMENTATION  PAGE 


2  GOVT  ACCESSION  NO. 


4  TITLE  'vnd  Suhil:lr) 


the:  turbulent  fluxes  or  momentum  and  sensible 

HEAT  OVER  THE  OPEN  SEA  DURING  MODERATE  TO 
STRONG  WINDS 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  recipient's  Catalog  number 


S  TYPE  OF  REPORT  &  PE  Rl  00  COVERED 


Ph.D.  Thesis 


6.  PERFORMING  QRG.  REPORT  NUMBER 


.  AUTHOR'S., 

William  George  Large 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Department  of  Oceanography,  L 
University  of  British  Columbia, 
Vancouver,  B.C.,  Canada.  V6T  1W5 


it.  controlling  office  name  and  address 
Office  of  Naval  Research, 

Code  48 1 , 

NSTL  STATION.  Mississippi  39529. 


M.  MONITORING  AGENCY  NAME  &  ADDRESS^//  different  from  Controlling  Office) 


e.  contract  or  grant  number^ 

N00014-76-C-0446 
NR  083-207 


io.  program  element,  project,  task 

AREA  &  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

August,  1979 


13-  NUMBER  OF  PAGES 
180 


15.  SECURITY  CLASS,  (ol  thie  report) 


16.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Uni imi+ed 


Unclassi f ied 


15a.  DECLASSIFICATION  /  DOWNGRADING 
SCHEDULE 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  It  different  from  Report) 


19.  KEY  WORDS  (Continue  on  reverse  aide  If  necemaary  and  Idantlfy  by  block  number) 

momentum  and  heat  fluxes  at  high  wind  speeds, 
drag  coefficient:  moderate  to  strong  winds, 
air-sea  interaction. 


|20.  ABSTRACT  (Conlinu*  on  reverse  aid*  ** 


identify  by  block  number) 


Two  systems  ter  remote  measurements  of  the  air-sea  fluxes  of  momentum, 
sensible  her:  _;i0  moisture  during  moderate  to  strong  winds  are  described. 
One  emp  j  the  dissipation  method  and  the  other  the  Reynolds  flux  or  eddy 
•—  uiation  method.  A  modified  Gill  propeller-vane  anemometer  is  the  velo¬ 
city  sensor  and  a  method  of  resolving  the  vertical  velocity  component,  that 
accounts  for  the  propel ler's  non-cosine  behavior  and  avoids  its  non-linear 
op.  at i ng  region,  is  derived.  The  dynamic  responses  of  the  sensors  are 


DD  ,'STti  1473 


.  .•■iiON  OF  I  NOVeSlbur-lLETE 
S/N  t  •  •  \F  014  6601 


Unclassi f ied 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wh»n  Dtf  Snl*r*W, 


.'CLNITV  CLASSIFICATION  O'  THIS  PAGE^Wim  Dmlm  Enured, i 


»  found  from  measurements  in  the  actual  turbulent  conditions  of  the  flux 
measurements. 

The  results  of  an  experiment  on  the  Bedford  tower,  a  stable  platform 
moored  in  59m  of  water  10  km  offshore,  are  presented.  Spectra,  cospectra, 
turbulence  statistics  and  transfer  coefficients  are  calculated  from  the 
Reynolds  flux  velocity  and  temperature  data  and  found  to  be  comparable 
to  previously  reported  values.  Simultaneous  dissipation  and  Reynolds  flux 
estimates  of  both  the  momentum  and  sensible  heat  fluxes  in  up  to  20  m/s 
winds  are  shown  to  be  in  excellent  agreement. 

Also  presented  are  the  results  of  a  second  experiment  where  the 
systems  were  deployed  on  the  weathership  CCGS  Ouadra.  A  comparison  of 
ship  and  tower  drag  coefficients  from  the  dissipation  system,  demonstrates 
that  the  Bedford  tower  is  essentially  an  open  ocean  site.  The  neutral 
drag  coefficient,  CDN,  is  found,  on  average,  to  be  nearly  constant  at 
1.1 4x 1 0—  3  for  winds  between  4  and  10  m/s  and  to  increase  almost  linearly 
to  about  2.18x10-3  at  26  m/s.  No  variation  with  either  fetch  (greater 
than  10  km)  or  stability  is  observed.  Dissipation  estimates  of  the 
sensible  heat  flux  from  a  wide  range  of  conditions  are  presented.  The 
neutral  transfer  coefficient,  CTN,  is  found,  on  average,  to  vary  from 
about  0.69x10“3  in  stable  stratification  to  1.08x10-3  in  the  unstable 
case.  An  increase  in  CTN  with  increasing  wind  speed  is  suggested  by 
only  some  of  the  data. 

Time  series  of  the  fluxes  are  used  to  investigate  additional  sources 
of  variation  in  the  transfer  coefficients.  Their  statistical  variability 
about  a  running  mean  is  seen  to  be  about  10$.  Evidence  is  presented  that 
indicates  that  persistent  departures  from  average  values  are  related  to 
sea  surface  conditions.  CDN  is  observed  to  be  significantly  smaller, 
on  average,  during  rising  winds  than  during  falling  winds  or  after  a 
change  in  wind  direction. 


Unc  I  ass i f i ed 


;<^V 

V"  THE  TURBULENT  FLUXES  OF  MOMENTUM  AND 
_§EHSIBLE  HEAT  OVRB  THE  OPEN^SEA  DOBING 
HODEBATE  TO  ST BO NG  RINDS. 


L &  •  If, «  1  * 


MIT 


by 


£ILLIAH  ^EORGE/XARGE < 

B.  A. Sc.  ,  University  of  British  Columbia,  1972 

— - '  A  THESIS  SOBHITTED  IN  PARTIAL  FULFILMENT  OF 

THE  REQUIREMENTS  FOR  THE  DEGREE  OF 
DOCTOR  OF  PHILOSOPHY 


in 

THE  FACULTY  OF  GRADUATE  STUDIES 
DEPARTMENT  OF  PHYSICS 
INSTITUTE  OF  OCEANOGRAPHY^  " 


DTIC 

* 1  f^TE] 
-E3  2  7  1980 


Re  accept  this  thesis  as  conforming 


to  the  required  standard 


C  Rilliaa  George  Large,  1979 


ii 


/ — 


ABSTBACT 

'  V 

Tvo  systeas  for  reaote  measurements  of  the  air-sea  flaxes 
of  aonentun,  sensible  heat  and  aoistare  daring  moderate  to 
strong  winds  are  described,  one  employs  the  dissipation  aethod 
and  the  other  the  Reynolds  flax  or  eddy  correlation  aethod.  A 
modified  Gill  propeller- vane  anemometer  is  the  velocity  sensor 
and  a  aethod  of  resolving  the  vertical  velocity  component,  that 
accounts  for  the  propeller's  non-cosine  behavior  and  avoids  its 
non-linear  operating  region,  is  derived.  The  dynamic  responses 
of  the  sensors  are  found  from  measurements  in  the  actual 
turbulent  conditions  of  the  flux  neasureaents. 

The  results  of  an  experiment  on  the  Bedford  tower,  a  stable 
platfora  aoored  in  59a  of  water  10  ka  offshcre,  are  presented. 
Spectra,  cospectra,  turbulence  statistics  and  transfer 
coefficients  are  calculated  froa  the  Reynolds  flux  velocity  and 
temperature  data  and  found  to  be  coaparable  to  previously 
reported  values.  Simultaneous  dissipation  and  Reynolds  flux 
estimates  of  both  the  monentua  and  sensible  heat  fluxes  in  up  to 
20  a/s  winds  are  shown  to  be  in  excellent  agreenent. 

Also  presented  are  the  results  of  a  seccnd  ex  peri lent  where 
the  systeas  were  deployed  on  the  weathership  CCGS  Quadra.  A 


comparison 

of  ship  and  tower 

drag 

coefficients 

froa 

the 

dissipation 

systea,  demonstrates 

that 

the  Bedford 

tower 

is 

essentially  an  open  ocean  site.  ^The  neutral  drag  coefficient, 

\ 

CDH ,  is  found,  on  average,  to  be  nearly  constant  at  1.14x10-* 
for  winds  between  4  and  10  a/s  and  to  increase  alaost  linearly 
to  about  2.18x10-*  at  26  a/s.  Bo  variation  with  either  fetch 


ill 


(greater  than  10  ka)  or  stability  is  observed.  Dissipation 
estimates  of  the  sensible  heat  flax  fros  a  vide  range  of 
conditions  are  presented.  The  neutral  transfer  coefficient, 
CTV,  is  found,  on  average,  to'  vary  fros  about  0.69x10***  in 
stable  stratification  to  1.08x10***  in  the  unstable  case.  An 
increase  in  CTV  vith  increasing  vind  speed  is  suggested  by  only 
soae  of  the  data. 

Tine  series  of  the  fluxes  are  used  to  investigate 
additional  sources  of  variation  in  the  transfer  coefficients. 
Their  statistical  variability  about  a  running  aean  is  seen  to  be 
about  10%.  Evidence  is  presented  that  indicates  that  persistent 
departures  fron  average  values  are  related  to  sea  surface 
conditions.  CDV  is  observed  to  be  significantly  saaller,  on 
average,  during  rising  vinds  than  during  falling  winds  or  after 
a  change  in  vind  direction. 
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CHAPTER  1 


INTRODUCTION 

This  thesis  describes  an  experimental  program  designed  to 
measure  the  turbulent  exchanges  between  the  open  ocean  and  the 
atmosphere  in  moderate  to  strong  (5-50  m/s)  winds.  Heasurements 
of  the  most  important  exchanges  have  been  the  subject  of  recent 
reviews,  in  which  they  are  parameterized  by  non-dimensional 
transfer  coefficients.  Recent  deter minaticns  of  the  drag 
coefficient,  which  is  used  to  express  the  momentum  flux  in  terms 
of  the  square  of  the  mean  wind  speed,  have  been  reviewed  by 
Garratt,  1977.  The  sensible  heat  and  moisture  fluxes  are 
parameterized  by  Friehe  and  Schmitt,  1976,  in  terms  of  the 
surface  -  air  temperature  and  humidity  differences,  respectively 
and  the  mean  wind  speed.  There  are  several  obstacles  that  make 
open  ocean  measurements  in  high  winds  difficult  and  that  have, 
therefore,  restricted  the  majority  to  low  winds  and  to  near  or 
onshore  platforms.  The  most  common  methods  of  obtaining  the 
fluxes,  the  Reynolds  flux  (or  eddy  correlation)  and  profile, 
work  best  on  stable  platforms  with  minimal  flow  distortion,  but 
these  conditions  are  not  easy  to  satisfy  during  storms  at  sea. 
Adverse  conditions  accompanying  high  winds,  cause  towers  to 
collapse  and  many  sensors  either  to  fail  completely  or  to  lose 
their  calibration. 

The  air-sea  energy  exchanges  are  involved  in  a  number  of 
important  processes,  including  the  large  scale  circulations  of 
the  ocean  and  atmosphere  and,  at  smaller  scales,  thermocline 
development  and  wave  generation.  By  extending  the  existing  data 
set  to  the  open  ocean  and  to  higher  wind  speeds,  this 
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experimental  program  should  be  relevant  to  the  stad;  of  these 
processes. 

Hodelling  and  predicting  large  scale  features  require  the 
fluxes,  which  are  too  difficult  and  costly  to  measure  directly 
on  this  scale,  to  be  calculated  from  easily  measured  quantities 
through  parameterizations  based  on  relatively  few  direct 
observations.  For  this  purpose,  extension  of  the  measured 
transfer  coefficients  to  20  n/s  ought  to  suffice,  because  higher 
winds  rarely  contribute  very  much  to  the  fluxes  averaged  over  a 
month  or  more  (Fissel  e&  ql. .  1977).  A  further  extension  to 
about  25  a/s  should  clearly  reveal  any  wind  speed  dependencies 
of  the  coefficients.  #  At  present  there  is  an  opinion  that,  in 
view  of  the  scatter,  a  constant  drag  coefficient  up  to  about  a 
14  n/s  wind  speed  is  appropriate  (Stewart,  1974),  while  Smith 
and  Banke,  1975,  and  others  find  a  significant  increase  with 
wind  speed.  The  average  stress  computed  from  either  type  of 
drag  coefficient  formulation  should  be  nearly  the  same,  because 
the  trend  is,  at  most,  small,  but  the  curl  of  the  wind  stress 
could  be  affected  to  a  greater  degree.  The  large  amount  of 
scatter,  typical  of  turbulence  measurements,  suggests,  that  in 
order  to  arrive  at  a  representative  picture  of  the  open  sea,  a 
great  deal  of  data  from  all  possible  conditions  are  required. 
Vith  a  large  data  set  it  would  be  possible  to  exanine  the 
effects  of  stability,  the  wave  field  and  other  sources  of 
variability  in  the  transfer  coefficients  apart  from  the  real 
statistical  scatter  and  systematic  instrumentation  errors. 


Continuous  records  over  a  long  period  of  tiae  should 
include  a  variety  of  local  and  short-lived  phenomena,  such  as 
frontal  passages,  for  which  a  large  scale  parameterization  may 
not  be  applicable.  In  some  cases  it  may  be  possible  to  find 
appropriate  transfer  coefficients  and  in  others  direct  stress, 
estimates  nay  be  the  simpler  approach.  The  capability  to 
operate  in  winds  above  40  m/s  should  allow  the  entire  tiae 
histories  (winds,  temperatures  and  fluxes)  of  most  storms  to  be 
followed.  Such  tiae  series  should  be  useful  for  the 
investigation  of  many  small  scale  processes. 

A  modified  Gill  propeller-vane  anemometer  proved  to  be  a 
very  suitable  velocity  sensor  for  this  study.  The  momentum  flux 
and  drag  coefficient  were  successfully  measured  in  26  a/s  winds. 
Temperature  and  humidity  sensors  were  housed  in  protective 
enclosures.  The  aicrobead  thermistors  were  often  broken  by 
spray  and  contaminated  by  salt,  however  they  did  survive  some 
high  winds  and  heat  transfer  coefficients  corresponding  to  large 
fluxes  were  obtained.  No  useful  humidity  data  were  ever 
recorded,  because  of  the  failure  of  several  types  of  sensors. 
For  open  sea  work  a  ship  is  the  most  convenient  platform  and  its 
motion  and  flow  distortion  can  be  tolerated  by  the  dissipation 
method  of  measuring  fluxes  (Pond  and  Large,  1978) .  This  method 
was  first  tested  on  a  stable  offshore  platform  where  its  results 
conpared  favourably  to  the  Reynolds  flux  method.  It  was  then 
employed  on  a  ship,  allowing  open  sea  data  to  be  collected.  In 
order  to  gather  as  much  data  as  possible,  the  instrumentation 
was  designed  to  record  continuously  for  a  month  or  more,  while 
operating  remotely. 


CHAPTER  2 


EXISTING  TgEOHX  AND  EXPERIMENTAL  RESDLTS 


2.1  Air-Sea  Interaction 

Exchanges  between  the  atmosphere  and  ocean  are  nost  easily 
■easured  in  the  atmospheric  surface  layer  where  the  transfer 
processes  are  doninated  by  turbulence.  Viscous  and  diffusive 
■olecular  transfers  are  negligible  in  this  layer,  which  begins  a 
few  centimeters  above  the  surface  and  extends  up  to  a  level 
where  the  earth's  rotation  and  the  geostrophic  pressure  gradient 
become  important.  Detailed  treatments  of  the  turbulent  flow  in 
the  layer  may  be  found  in  Lumley  and  Panofsky,  1964,  Honin  and 
Yaglom,  1965  and  1967,  and  Kraus,  1972.  This  chapter  and  Busch, 
1977,  are  specifically  concerned  with  the  theory  related  to 
turbulent  flux  measurements  and  their  interpretation.  Here  the 
sensible  heat  flux  and  the  momentum  flux  are  treated  explicitly 
and  the  theory  is  also  extended  to  any  passive  atmospheric 
scalar  quantity,  R.  Following  Reynolds'  convention,  the 
turbulent  properties  are  partitioned  into  a  mean  (  <  >  denotes  a 
time  average)  and  a  fluctuation  (lower  case  symbols) .  The 

“A  *<B 

components  of  the  instantaneous  wind  vector,  V  =  Oi+Tj+Hk,  where 
i,  j”,  and  k  are  unit  vectors  of  an  x-y-z  coordinate  system, 
become  0=<D>+u,  V=<V>+v  and  W=<l>+w.  The  usual  orientation  of 
the  axes  puts  "k  vertically  up  and  i  along  the  mean  horizontal 
wind  vector  such  that  the  mean  cross-stream  and  vertical 
components,  <V>  and  <N>,  are  both  zero  (Burling  and  Stewart, 
1967).  Similarly  any  scalar  field  R  becomes  <R>*r  (and  the  air 
temperature  T=<T>+t  and  the  air  pressure  P»<P>+p).  By 
definition  <u>,  <v>,  <w>,  <t>  and  <r>  are  all  zero. 
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It  is  the  fluctuating  vertical  velocity  which  bodily 
transports  fluid  properties  up  and  down,  giving  rise  to  the 
Seynolds  fluxes  defined  by: 

Hoaentua  flux  <uv> 

Sensible  heat  flux  Hs  »  p  Cp  <wt>  (2. 1) 
any  scalar  flux  Hr  *  <vr>  , 

where  ^  is  the  sean  air  density  and  Cp  is  the  specific  heat  at 
constant  pressure.  Since,  <H>*<w>=0,  the  fluctuating  quantities 
in  2.1  can  be  replaced  by  their  instantaneous  values  0,  T  and  B. 
In  this  coordinate  systea  <vw>  should  tend  to  zero  with  a  long 
enou'gh  averaging  period,  so  ~C  represents  the  total  aoaentua 
flux.  Since  it  gives  rise  to  a  force  in  the  direction  of  the 
aean  wind  on  a  unit  area  of  underlying  surface,  ~C  is  also 
referred  to  as  the  Beynolds  stress.  Hs  is  a  turbulent  heat 
transfer  (positive  up) •  The  aoisture  flux,  also  an  iaportant 
air-sea  exchange,  is  expressed  by  siaply  substituting  absolute 
huaidity  for  B.  Siailarly  gas  fluxes  such  as  carbon  dioxide  aay 
also  be  considered.  Often  the  terns  aoaentua  flux  and  sensible 
heat  flux  refer  to  the  kinenatic  fluxes  <uw>  and  <wt>. 

The  Beynolds  fluxes  arise  in  the  eguation  of  notion  and  in 
scalar  conservation  equations  where  their  surface  values  becoae 
iaportant  boundary  conditions  for  both  the  ataosphere  and  ocean. 
In  the  equation  for  <U>  in  the  boundary  layer,  the  Coriolis 
force  due  to  the  cross-strean  coaponent,  v,  is,  on  average, 
zero.  The  equations  for  the  aean  flow  and  aean  teaperature  in 
the  surface  layer,  assuaing  horizontally  hoaogenous  turbulence 


i 


(Bosch,  pp  74  and  75)  ,  bat  retaining  terns  with  the  largest  aean 
horizontal  gradients,  are: 


e 


a<p> 

~3x~ 


d<T>  ♦  <D>  d<T>  ♦  3  (Hs  *  BT) 

3t  5x  5z  ?  Cp 


0 


(2.2) 


where  BT  is  the  vertical  heat  transfer  doe  to  radiation.  Bith 
negligible  horizontal  pressure,  horizontal  teaperatnre  and 
vertical  radiative  flux  gradients  and  a  steady  aean  state,  the 
turbulent  fluxes  are  constant  throughout  the  layer.  It  is  then 
possible  to  neasure  the  surface  flaxes  above  wave  influences  at 
a  convenient  height,  Z.  However,  away  fron  the  surface, 
rotation  and  the  large  scale  horizontal  gradients  eventually 
becone  influential. 

In  steady  flow,  the  neasured  stress,  T"(Z) ,  is  less  than  the 
surface  stress,  Z~o.  if  the  difference  at  a  height,  he,  is  10X 
(Lusley  and  Panofsky,  1964,  arbitrarily  use  20X) ,  then  the  flow 
below  can  be  regarded  as  being  a  "constant  flux"  or  "constant 
stress"  layer.  Effectively,  he  is  taken  to  be  the  upper  lisit 
of  the  ataospheric  surface  layer.  In  aid- latitudes  winds  above 
the  surface  layer  are  governed  by  the  geostrophic  balance, 

f  Og  »  1/£  3<P>/  3n  , 

where  f  is  the  Coriolis  paraaeter  (about  1x10-*  s-*),  Og  is  the 
geostrophic  wind  and  n  is  a  horizontal  coordinate  perpendicular 
to  the  Og  direction.  Observations  have  shown  Og  to  be  about 
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1.3  <H>  and  about  16°  to  tbe  right  of  the  i  direction  (Deacon, 
1973).  Substituting,  d<P>/  dx  S’  sin  16°  d<P>/  8n,  into  the 

first  of  2.2  yields, 

3<0>  Sf  i  i?  ♦  1.3  f  <0>  sin  16®  . 

St  5  3z 

In  steady  flow  a  10%  reduction  is  found  when 


0. 1  *  To  -  TV 


Ta-Z.J.LZL 

To 


C  <P>  sin  16®  Z  , 


so  he  *  0.1  <uw>  *  2790  seconds  <uw>.  (2.3) 

1.3  f  <0>  sin  16®  <0> 


Hea sore aents  over  the  sea  at  10a  have  alios t  always  shown 
<uw>/<0>  to  average  aore  than  10~*  <D>  (drag  coefficient  > 
1x10-*,  Garratt,  1977),  setting  14a  as  a  lower  liait  of  he,  when 
<0>»5  a/s.  In  unsteady  flow,  d<U>/  St  can  easily  be  the  sane 
order  of  aagnitude  as  (1.3  f  <0>  sin  16°)  *  0.13  <D>  /hour.  On 
the  rising  wind  the  acceleration  is  down  the  pressure  gradient 
and  he  rises,  because  a  snail er  stress  gradient  is  sufficient  to 
balance  2.2.  On  the  falling  wind  the  loss  of  flux  with  height 
is  enhanced  by  the  deceleration  and  in  this  situtaion  he  aay  be 
considerably  lower.  These  arguaents  are  possibly  good  to  a 
factor  of  2  despite  the  neglect  of  horizontal  advection  terns, 
relative  to  3<?>/dx.  They  serve  to  point  out  that,  with  the 
ease  surface  stress,  the  aeasured  10a  stress  aay  be  greater  on 
the  rising  wind  than  during  falling  winds,  when  it  could  be  less 
than  the  surface  stress  by  5  to  10%. 


The  height  to  which  the  sensible  heat  flux  resains  within 
10X  of  its  surface  value  is  not  obvious  and  it  say  sosetises  be 
below  usual  seasuresent  heights.  Since  ^<T>/dx  is  not  sinply 
related  to  the  large  scale  pressure  gradient,  it  is  not  possible 
to  scale  the  vertical  divergence  of  Bs  to  a  Coriolis  tern  as  was 
done  for  the  Beynolds  stress.  Neglecting  radiation,  a  10% 
change  in  Hs  is  found  at  a  height  he  given  by 

he  ■  0.1  <wt>  (fiT/fit)-*  , 

where  6t  is  the  change  in  tenperature  due  to  the  vertical  heat 
flux  divergence,  during  a  tine  interval  fit.  Ilea  sure  nents  over 
the  sea  show  that  <wt>  is  of  order  10-*  <U>  AT  (Friehe  and 

Schnitt,  1976) ,  where  AT  is  the  tenperature  difference  between 
the  sea  surface  and  atnosphere.  To  keep  he  above  lOn,  with 
<U>  AT  as  low  as  10  °Cn/s  reguires 

$T/  fit  <  0.36  °C/hour, 

which  nay  not  necessarily  be  satisfied.  Rhen  d<T>/dt  over  the 
sea  is  large,  presunably  horizontal  advection  is  the  najor 
contributor  with  fiT/  fit  hopefully  renaining  snail.  Rhen  the 
tenperature  is  steady  <0>  fi<T>/dx  nust  balance  dHs/dz  in  2. 2, 
giving 

he  «  10-*  AT  (  8<T>/fix  )-* 

and  inplying  that  he  is  only  above  10a  when  the  horizontal 
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teaperature  gradient  is  rather  saall,  less  than  0.01  °C/ka  for  a 
AT  of  1°C.  i  farther  coaplication  arises  daring  low  winds  when 
infrared  absorption  by  water  vapour  nay  prodace  a  radiative  flax 
divergence  that  tends  to  caase  Hs  to  increase  with  height 
(Bosch,  1977).  This  effect  enhances  the  positive  flax  gradient 
when  AT  is  negative,  bat  reduces  the  loss  of  sensible  heat  flax 
with  height  when  AT  is  positive.  It  will  be  assaaed  that  Bs 
anywhere  in  the  "constant  stress"  layer  over  a  temperate  sea  is 
eguivalent  to  the  sarface  flax  to  within  the  accuracy  of  the 
aeasureaent.  Hopefully,  this  assn apt ion  is  true  on  average,  but 
verification  would  require  a  direct  aeasareaent  of  Hs  at  two 
levels. 


2.2  Honin-Obukhov  Siailarity  Theory 

The  understanding  of  the  turbulent  atmospheric  sarface 
layer  is  largely  doe  to  flonin-Obukhov  siailarity  theory  (Honin 
and  Tagloa,  p  425  ff) .  The  theory  assuaes  that  turbulent 
characteristics  depend  only  on  a  few  physical  paraaeters,  which 
facilitates  the  application  of  diaensional  analysis,  ibove  the 
direct  influence  of  the  botton  boundary  the  iaportant  paraaeters 
are  the  height  (the  only  spatial  variable  left  in  the  assuaed 
horizontally  hoaogenous  turbulence),  the  air  density,  the 
turbulent  transports  and  the  stability  of  the  air  coluan.  The 
supposed  height  independence  of  the  fluxes  naturally  leads  to 
the  following  scales  which  incorporate  the  transports  through 
the  layer  and  the  density: 


friction  velocity  u*  *  (Tb/£)  V*  *  |-<uw>|  V* 

temperature  scale  t*  *  -<wt>/  Ku*  (2.4) 

scalar  scale  r*  *  -<wr>/  Ku*  , 

where  von  Kansan' s  constant,  K,  is  included  to  simplify  later 

equations,  but  other  scales 

T*  =  -<wt>  /  u* 
and  B*  =  -<wr>  /  u* 

are  sometimes  used.  &n  appropriate  stability  parameter,  Z/L,  is 
obtained  from  the  ratio  of  the  convective  or  buoyant  turbulent 
kinetic  energy,  B,  produced  in  a  non-neutral  air  column,  to  the 
purely  mechanical  production  in  the  equivalent  neutral  case,  Po. 
It  will  be  seen  in  section  2.5,  that  ^ 

B  =  g  <w  Tv>  /  To  and  Po  *  u*V  KZ  ,  (2.5) 

where  Tv,  the  virtual  temperature  in  degrees  Kelvin,  accounts 
for  the  fluctuating  temperature  and  moisture  contributions  to 
the  fluctuating  density,  g  is  gravitational  acceleration  and  To 
is  the  local  average  virtual  temperature.  The  important  scale 
is  the  Honin-Obukhov  length,  L,  whose  magnitude  gives  the  height 
at  which  Po  *  |B|,  thus 

-Z  =  B  and  L  *  -u»3  To  . 

L  Po  K  g  <v  Tv>  (2.6) 


In  neutral  stratification  <v  Tv>  and  Z/L  go  to  zero  while  L 


11 


approaches  infinity  and  the  sign  is  chosen  to  make  L  and  Z/L 
positive  in  stable  conditions.  Oinensional  analysis  predicts 
that  all  turbulent  functions  non-dimensionalized  by  Z,  u*,  r* 
and  L  should  be  functions  of  the  only  possible  dinensionless 
group,  Z/l.  Note  that  each  additional  scalar  adds  both  a  scale 
and  a  diaension  to  the  problen.  Use  of  spectra  introduces  a 
frequency,  f,  making  a  further  dinensionless  group,  fZ/u*, 
possible,  but  because  of  the  difficulty  in  obtaining  u*  the 
dinensionless  natural  frequency  n  *  fZ/<U>  is  usually 
substituted.  Nornalized  spectra  and  cospectra  can  be  regarded 
as  non-dimensional  turbulent  functions  and  therefore  should  be 
functions  of  both  n  and  Z/L,  whereas  their  integrals  such  as 
(ou) *  *  <u*>  and  (tft) 2  ®  <t*>  should  depend  only  on  Z/L.  Nith 
the  iaportant  paraaeters  coaaon  to  all  surface  layer  flows,  the 
structure  of  the  turbulence,  according  to  this  theory,  nust 
always  be  "similar”,  with  any  dependencies  cn  Z/L  and  n  being 
universally  applicable. 

in  iaportant  consequence  of  similarity  theory  is  the 
logarithaic  profile  of  the  aean  wind  and  aean  scalars. 
Oinensional  considerations  lead  directly  to  the  forms: 


KZ  a<0>  «  *n  (Z/L) 

u*  dz 


and  d<R> 

r*  dz 


*r(Z/L)  , 


(2.7) 


where  von  !Caraan*s  constant  sets  t*n*1  at  neutral  stability  and 
has  already  been  included  in  the  definition  of  r*  so  that 
fc(0)*1.  Heasured  values  of  K  vary  between  about  0.35  and  0.02 


(Busch,  1977),  so  at  least  a  5%  error  in  its  customary  value  of 
0.40  aust  be  allowed.  In  a  review  of  flux-profile 
relationships.  Dyer,  1974,  suggests  that  the  best  forms  of  the 
universal  functions  are 

0  <  Z/L  <0.2:  fa  =  fr  *  1  ♦  5  Z/L 

-1.0  <  Z/L  <  0  :  fa  -  (1  -  16  Z/l)-»/* 

f r  *  (1  -  16  Z/L)-»/2  , 

where  the  unstable  case  is  the  Businger-Dyer  representation 
provided  by  Dyer  and  Hicks,  1970.  The  mean  values  at  a  height 
Z,  OZ  and  BZ,  are  found  by  integrating  2.7,  (Paulson,  1970), 

OZ  =  (u */  K)  •  [  In  (Z/Zo)  -  Y4( Z/L)  ] 

bz  »  rsfc  ♦  r*  [ln(Z/Zor)  -  Vt( z/L)  ],  (2.8) 

where  V^Z/L)  =  £  1  -  f  (?)  3  /£ 

stable:  ¥m(Z/L)  =  >r(Z/l)  =  -5  Z/L 

unstable:  Ya  (Z/L)  *  2  ln[(1+X)/23  ♦  ln[  (1+X*)  /2  ] 

-  2  tan-*X  ♦  7T/2 
"Yr  (Z/L)  «  2  ln[  (1*X8)  /2] 

I 

with  X  ■  (1-16  Z/L)  . 

it  neutral  stability  the  integral  vanishes,  leaving  Vt0)«0. 
The  constants  of  integration  Zo  and  Zor,  assumed  to  he  such 


smaller  than  Z,  are  the  roaghness  lengths,  which  fall;  describe 
the  purface  as  "seen"  by  the  turbulence,  but  they  are  not  simply 
related  to  sea  surface  parameters  such  as  wave  height  and 
temperature.  They  need  to  be  included  in  the  dimensioned, 
analysis  only  near  the  surface  where  they  set  the  magnitude,  but 
not  the  structure,  of  the  turbulence  throughout  the  layer. 


2.3  Bulk  Aerodynamic  Parameterizations 

The  fluxes  are  parameterized  in  terms  of  the  mean  wind,  the 
sea  surface-air  temperature  difference  and  .  surface-air  mean 
scalar  difference,  AR,  with  the  bulk  aerodynamic  formulae: 

Z/g  =  u**  »  -<uw>  »-  CD  <0>2 

Hs/QCp  *  -  K  u*  t*  «  <wt>  »  CT  <U>  nr  (2-9) 

Hr  ■  -  K  u*  r*  *  <wr>  »  CR  <D>  AR, 

with  AT*TSCP-TZ  and  AR«RS7C-RZ,  where  TS7C  and  RSPC  are  the  mean 
surface  temperature  and  scalar  Talues,  respectively.  The 
non-dimensional  transfer  coefficients  CO  and  CT  are  the  drag 
coefficient  and  Stanton  number  while  the  corresponding 
coefficient  of  moisture  transport,  CS,  is  the  Oalton  number. 
Their  dependence  on  stability,  roaghness  and  height  is  evident 
from  2.8,  but  Zo  over  the  sea  has  a  complicated  functional  fora 
(Burling  and  Stewart,  1967).  However,  they  can  also  be 
determined  experimentally  from  2.9  using  measured  flaxes  and 
bulk  quantities  and  such  calculated  coefficients  provide  a 


convenient  leans  of  coaparing  flax  measure sent s.  To  eliminate 
the  variation  with  height  they  are  commonly  evaluated  at  10a  as: 

CIO  =  -<uw>  /  (010)* 

CT10  *  <wt>  /  [010  (TSFC-T10)  ]  (2.10) 

CB  10  =  <wr>  /  [010  (SSFC-B10)  ]  . 

Equation  2.8  shows  the  wind  speed,  temperature  and  scalar  aeans 
at  10m  (010,  T10  and  B10)  to  be: 

010  =  0Z  -(uV  K)[ln  (Z/10b)  -Yb(Z/L)  ♦Yb(10i/L)  ] 

T10  *  TZ  -  t*  [  In  (Z/IOm)  -  (Z/L)  +  Y£(10n/L)]  (2.11) 

B10  ■  BZ  -  r*  [In  (Z/IOn)  -  Yi  (Z/L)  ♦  Y£(10m/L)]  . 

For  conparative  purposes  it  is  convenient  to  eliainate  the 
stability  dependence  by  evaluating  the  roughness  lengths  froa 
2.8  and  using  thea  to  find  the  coefficients  in  the  equivalent 
neutral  case  at  10a  froa: 

CDN  =  X*  /  [In (lOm/Zo) ]* 

CTN  »  K*  /  [  ln(10a/Zot)  •  ln(10a/Zo)  ]  (2.12) 

CBN  *  K*  /  [  In  (lOn/Zor)  •  ln(10a/Zo)]  . 

The  neutral  coefficients  should  be  constants  over 
hoaogeneous  terrain  where  the  roughness  lengths  can  be  regarded 
as  constant.  This  prediction  has  been  verified  over  land  where 
Zo,  for  exanple,  is  exclusively  deterained  by  topography  and 
vegetation.  It  is  not  unreasonable  to  expect  this  concept  to 


15 


work  even  better  over  the  sea  where  there  is  only  one  type  of 
surface.  However,  direct  aeasureaents,  reviewed  by  Garratt 
(1977)  are  very  scattered  and  indicate  that  the  neutral  drag 
coefficient  varies  with  wind  speed  and  is  such  ssaller  than  that 
found  over  land.  This  result  iaplies  that  there  are  additional 
iaportant  paraaeters  deteraining  the  roughness  of  the  sea 
surface,  in  obvious  difference  between  a  land  and  sea  boundary 
is  surface  gravity  waves,  so  it  seeas  appropriate  to  add  the 
acceleration  due  to  gravity,  g,  to  the  problea.  This  leads  to 
the  Charnock  (1955)  diaensionless  group  for  flow  near  the  waves, 

Zo  g  /  u**  *  «<  (2.13) 

where  *0.0144  is  suggested  by  Garratt  (1977).  Stewart,  1974, 
notes  that  for  winds  below  10  a/s  the  Charncck  representation 
with  •<  constant  predicts  a  aore  rapid  increase  in  CDH  with  wind 
speed  than  is  indicated  by  the  results  of  Brocks  and 
Rruggeraeyer  (1970)  and  this  feature  is  also  present  in 
Garratt* s  review.  It  appears,  therefore,  that  aore  surface 
paraaeters  aay  be  iaportant  to  this  aspect  of  turbulent  flow. 
Stewart  discusses  the  possible  roles  of  surface  tension 
(capillary  waves) ,  the  length  and  phase  speed  of  the  longest 
excited  waves,  the  wave  slope  and  the  total  wind  generation 
force  which  is  proportional  to  (wind  speed  -  wave  speed). 
Kitaigorodskii  and  Zaslavskii,  1974,  consider  the  phase  speed  of 
the  doainant  wave  and  a  purely  viscous  aoaentua  flux.  Burling 
and  Stewart,  1967,  exaaine  the  iaplications  of  dependency  on 
various  aoaents  of  the  wave  spectrua.  flith  the  roughness 


16 


lengths  possibly  depending  on  many  parameters  it  appears  that  a 
rather  detailed  knowledge  of  the  sea  surface  would  be  reguired 
before  the  turbulent  fluxes  over  the  sea  could  be  found  froa  Zo 
and  Zor. 

The  paraaeterizations  can  also  be  regarded  as  eapirical 
foraulae.  in  experimental  formulation  of  the  neutral 
coefficients  at  10a,  2.12,  for  example,  would  allow  the  fluxes 
to  be  estimated  froa  mean  or  bulk  quantities,  DZ,  TZ  and  TSFC 
(4T  =  TSFC-TZ) ,  with  Z/l,  if  known,  providing  a  stability 
correction. 

The  momentum  flux  can  be  found  froa  2.9  by  finding  the  drag 
coefficient  at  the  measurement  height,  Z,  the  wind  speed,  OZ, 
and  the  stability,  Z/L,  from  CDN.  Elimination  of  Zo  froa  2.12 
and  2.8  and  substitution  into  2.9  leaves, 

CD  *  CDN  {1*  CDN*/*  K-*  (In  (Z/10 m) -Ym  (Z/L)  )  }-*  ,  (2.14) 

where  ln(Z/10m)  andVm(Z/L)  describe  the  variation  of  the  drag 
coefficient  with  height  and  stability,  respectively.  If  CDN  is 
given  as  a  function  of  the  10m  wind,  then  OZ  must  first  be 
shifted  to  10a  before  a  drag  coefficient  can  be  determined. 
Substituting  u*/0Z  =  CD*/*  =  C10V*  OIO/OZ,  into  2.11  and  solving 
for  OZ/OIO  leaves, 

OZ/010  *  1+  C10V2  K“*  [  In  (Z/1 0a) -  Y'm  (Z/L)  ♦  Vm  ( 10  a/L)  ]. 
The  term  in  square  brackets,  K(Z,Z/L),  is  usually  dominated  by 


ln(Z/10a),  which,  with  CIO V*/K  *  0.1,  wakes  010  aboat  1 0S 
larger  or  saaller  than  02  for  a  Z  of  3.7a  and  27a,  respectively. 
Throughout  this  height  range  neglect  of  the  stability  portion  of 
K (Z,Z/L)  introduces  an  error  in  010  of  less  than  2X  for  >1.0 
<Z/L<  0.08,  however  at  Z*27a  it  decreases  010  by  about  an 
additional  6%  with  Z/L*0. 2  and  3%  at  Z/L*0.1,  so  it  is  not 
ignored.  CIO  is  egnivalent  to  CD  at  Z*10a,  so  2.14  gives 

C10V*  *  CDHV*  {1  -X-1  CDH*/*  Va(10a/L) }-»  . 

Since  Ya(10a/L)  ranges  fron  about  1  at  Z/L*  -1  to  about  -1  at 
Z/L*0.2,  taking  the  tera  in  curly  brackets  to  be  1.0  introduces 
an  error  of  only  IX  in  010,  which  is  within  usual  aeasureaent 
error,  when  it  is  calculated  froa  0Z,  CDH  and  Z/L  using 

010  *  OZ  [1  ♦CDS*/*  *-»  K  (Z,Z/L)  J*1.  (2.15) 

Should  CDH  itself  depend  on  010,  2.15  aay  have  to  be  solved  with 
an  iterative  technique. 

If  the  stability  is  in  the  range  -1.0  <Z/L<  0.2,  but  is 
unknown  and  assuaed  to  be  neutral,  errors  arise  froa  an 
inaccurate  010,  in  finding  CDI  froa  010  and  through  2.14, 
because  CDH  is  not  shifted  to  the  proper  stability.  The  total 
error  in  010*s  found  froa  2.15  should  be  less  than  10X.  Saith 
and  Banke,  1975,  report  a  drag  coefficient  equal  to  0.00061  ♦ 
0.000075  010,  so  a  10X  error  in  010  at  20  a/s,  reduces  to  7X  in 
CDH  and  the  aonentun  flux,  which  is  less  than  the  error 
associated  with  direct  aeasureaent s.  With  Yn (Z/L)  ranging  froa 
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1.0  to  -1.0,  taking  it  to  be  zero  could  lead  to  a  20%  error  in 
the  momentum  flux.  Only  in  the  range  -0.25  <Z/L<  0.1  does  the 
total  error  in  assuming  neutral  stability  remain  less  than  10%. 

If  CD  and  CDN  are  available,  the  analogous  procedure  can  be 
followed  to  find  the  sensible  heat  flux  fron  TZ,  TSFC,  and  0Z, 
using  2.9.  The  Stanton  nunber,  CT,  is  expressed  as  a  function 
of  CTN,  CDN ,  CD,  and  2/ 1,  by  eliainating  Zo  and  Zot  froa  2.12 
and  2.8,  then  substituting  into  2.9: 

CT  *  CTN _ fCD/CDN)  & _ _ . 

[  1+CTN  CDN -*/*  (In  (Z/IOaj^nZ/l) )  ]  (2.16) 

The  major  error  coaes  froa  the  uncertainty  in  CTN.  The  errors 
in  CDN  and  CD  due  to  010  cancel  in  (CD/CDN)  and  a  10%  error  in 
CDN  and  CTN  introduces  a  2%  error  in  the  denoainator.  Since 
Y%(Z/L)  is  not  very  different  froa  YU  (Z/L) ,  assuning  neutral 
stability  causes  about  the  sane  error  in  CT  as  in  CD.  However, 
this  assumption  is  never  needed,  because  a  aeans  of  estimating 
Z/L  from  UZ  and  AT,  which  are  required  in  2.9,  is  developed  in 
section  2.6.  The  uncertainty  of  this  estimate  (section  4.2) 
should  result  in  about  a  5%  error  in  CT  and  CD. 

The  sensible  heat  flux  is  sometimes  parameterized  as 

<vt>  =  a  010  (TSFC-T 10)  ♦  b  , 

where  a  and  b  are  experimentally  deterained  functions  of 
stability  and  wind  speed.  The  air  temperature  at  10a,  T10,  can 
be  found  as  follows;  froa  2.11 


TZ  -  T10  ■  t*  [ In  (Z/10a) -  (Z/L)  ♦>£(10b/L)]  , 

where  the  tera  in  sguare  brackets,  to  be  denoted  as  Kt(Z,Z/L), 
behaves  as  K(Z,Z/L)  and  the  stability  contribution  can  again  be 
ignored  in  near  neutral  conditions.  Substituting  for  t*  using 
2.9  gives 

T10  ■  TZ  ♦  (CT/COV*)  (TSPC-TZ)  Kt(Z,Z/L).  (2.17) 


2.4  The  Beynolds  Flux  Method 

The  Beynolds  flux  or  eddy  correlation  nethod  is  the  nost 
direct  aeasureeent  of  the  fluxes  and  has  been  employed  over  the 
sea  by  Pond  si  al..  1971,  Ricks,  1972,  Snith  and  Banke,  1975, 
and  others.  It  involves  integrating  the  cospectra  of  a  and 
either  u  or  r  to  obtain  the  covariances  and  hence  fluxes.  The 
spectral  for as  of  the  covariances  are: 


<uw> 

»  jVuw(f) 

df 

<wt> 

-  /*vt(f) 

df 

(2.18) 

<wr> 

«  /fwr(f) 

df. 

In  practice  the  cospectra  are  deterained  by  digital  fast  Fourier 
transfora  technigues,  which  give  discrete  values  of  f(f)  at 
intervals  of  f  such  that  ♦  (£)  Af  gives  the  covariance  in  a  band 
centered  at  £  of  width  Af.  The  highest  freguency  coaputed,  the 
lyguist  freguency,  fny,  is  set  by  the  digitization  period  at: 
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(fny* 1/2st) .  Contributions  fro»  higher  frequencies  can  be 
allowed  to  alias  back  below  fny  so  that  the  effective  upper 
liait  of  the  integration  is  increased  (  to  about  2  fny  in  the 
systea  described  in  section  3.3) .  The  contributions  to  the 
covariances  from  natural  frequencies,  n=fZ/<U>,  greater  than  1 
are  only  a  few  per  cent  of  the  total.  Bith  20  a/s  winds, 
aeasurenents  at  10  a  will  therefore  include  nost  of  the  high 
frequency  contributions  if  At  is  no  longer  than  Z/<U>,  about  0.5 
seconds.  The  lowest  calculated  frequency,  fl,  is  the  reciprocal 
of  the  duration  of  the  aeasurenent  and  f^(fl)  includes  the 
covariance  down  to  fl  /2.  Some  cospectra  are  still  non-zero  at 
n-0.001,  so  in  order  to  be  able  to  neasure  fluxes  in  5  a/s  winds 
at  10  a  the  saaples  must  be  taken  for  about  Z/(0. 002<0>)  -  1000 
seconds  or  about  15  ainutes.  However  f(fl)  is  not  a 
statistically  well  deterained  quantity  and  in  practice  at  least 
3  sequential  deterainations  need  to  be  averaged,  requiring  the 
duration  of  a  flux  run  to  be  at  least  45  ainutes. 

Unfortunately,  aeasurenents  cannot  be  extended  to  nore  than 
about  1  hour  because  stationarity  begins  to  be  lost  as  a  new 
flow  situation  develops.  The  low  frequency  contributions  to  the 
fluxes  are  often  not  well  established  by  the  Beynolds  flux 
nethod.  It  is  also  clear  that  another  disadvantage  to  this 
aethod  is  the  large  aaount  of  data  required.  For  exanple,  a 
single  <uw>  estinate  fron  a  45  ninute  run  with  ot*0.5  seconds 
requires  about  five  thousand  digitizations  of  each  variable.  Of 
course  the  spectra  of  the  neasured  quantities  nay  also  be  found 
froa  the  sane  data  and  are  often  a  great  advantage  in  checking 
sensor  perforaance.  The  sensitivity  of  the  Beynolds  flux  aethod 
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to  instrument  orientation  is  a  groat  handicap  on  ships  and 
buoys,  whoso  notion  and  soan  tilt  offoct  tho  soasorod 
covariances.  A  ono  degree  error  in  the  Kean  tilt  cf  an 
aneaoaeter  nay  indace  errors  in  <aw>  in  the  order  of  10X  (Pond, 
1968) .  It  is  possible  to  neasare  the  instantaneous  platform 
notion  and  to  correct  the  velocities  point  by  point  (Bitsata  and 
Fujitani,  1974),  bat  this  greatly  increases  the  recording 
regoirenents  and  is  not  a  very  practical  neans  of  obtaining 
large  anoants  of  open  sea  flax  neasarenents. 

The  Beynolds  flax  nethod  is  not  very  applicable  to  renote 
open  sea  operation,  bat  it  has  becone  the  standard  to  which 
other  nethods  are  conpared  either  directly  or  through  the 
calculated  transfer  coefficients. 


2.5  The  Dissipation  Bethod 

Following  Deacon's,  1959,  suggestion  the  dissipation  nethod 
has  been  enployed  in  open  sea  conditions  by  Pond  ££  al. .  1971, 

Bucfcnitz,  1976,  Dennan  and  Biyake,  1973,  and  others.  It  is  a 
very  attractive  nethod  because  it  does  not  involve  an  explicit 
neasorenent  of  the  vertical  velocity,  allowing  noving  platforns 

to  be  used  and  reducing  neasurenent  errors.  In  addition,  flow 

* 

distortions  that  would  hinder  covariance  neasarenents,  can  be 
tolerated.  Instead,  the  najor  sources  of  error  arise  in  the 
uncertainty  of  various  constants  and  in  the  necessary 
assunp-tions.  in  the  cited  studies  its  results  have  been 
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compared  to  direct  Reynolds  flax  measurements  up  to  aoderate 
winds,  but  further  coaparison  at  high  wind  speeds  is  still 
necessary. 

In  the  case  of  the  aonentua  flax  the  aethod  steas  froa  a 
consideration  of  the  balance  of  turbulent  kinetic  energy  per 
unit  aass,  e  =  (u2  ♦  f*  ♦  w*) /2,  in  horizontally  hoaogeneous 
flow  (Busch,  1977), 

^<e>  =  u*2  d<0>  ♦  g  <w  Tv>  -  £  -  3  [<we>  ♦  J[  <wp>] 

5t  3z  ~To  Sz  £ 

*  P  ♦  B  -  6  -  D  ,  (2.19) 

where  Tv  is  the  virtual  teaperature  in  degrees  Kelvin  and  To  its 
local  average  and  p  is  the  fluctuating  pressure.  Turbulent 
fluctuations  are  chiefly  produced  by  aechanical  interactions  of 
the  Reynolds  stress  with  the  nean  flow  represented  by  the  first 
tera,  P=u*2  d<0>/dz  and  lost  at  small  scales  to  nolecular 
dissipation,  €.  Froa  Lumley  and  Panofsky,  1964  p95,  B  is 
recognizable  as  the  loss  or  gain  due  to  buoyancy  referred  to  in 
section  2.2.  D  is  the  sub  of  two  vertical  divergences:  the 
first,  of  the  turbulent  kinetic  energy  flux  <we>  and  the  second, 
of  the  work  done  per  unit  area  by  the  fluctuating  pressure, 
<wp>/g.  These  are  referred  to  as  the  turbulent  and  pressure 
transports  of  kinetic  energy,  respectively.  The  complete 
divergence  tera  has  been  investigated  by  BcBean  and  Elliot, 
1975.  In  this  work  <we>  and  <wp>  were  aeasured  over  land  at  one 
height  for  a  range  of  Z/L  values,  k  fit  between  -0.31  <Z/L< 
0. 12  gave 


2.3  Z/L 


0.20 
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<Wp>  /  {£  U*») 


Their  <we>  resalts  were  plotted  with  those  of  Garratt,  1972,  and 
Banke  and  Saith,  1973.  In  view  of  the  scatter  it  is  not 
unreasonable  to  assuie  a  relation 


<we>  /  a**  *>  -2.3  Z/L  *  constant. 


Extensive  aeasareaents  over  land  of  the  tarbalent  transport  tern 
have  been  aade  by  Wyngaard  and  Cot6,  1971,  over  a  wider  range  of 
stabilities.  In  unstable  conditions  their  results  can  be 
expressed  as 

<we>  /  u*a  ■  -2.5  Z/L  ♦  constant. 

The  coabined  experimental  evidence  in  the  range  -1  <Z/L<  0.1 

suggests  that  to  a  very  good  approximation 

<we>  ♦  <wp>/£  *  a  constant 

Differentiation  by  z  implies  that  on  average,  the  kinetic  energy 
gained  through  pressure  transport  nearly  balances  that  lost  by 
turbulent  transport,  that  is,  D*0.  Wyngaard  and  Cot6  also 
conclude  that  the  effects  of  horizontal  inhoaogeneities  and 
non-stationarity  are  negligible  by  more  than  two  orders  of 
magnitude,  that  is,  d<e>/dt  *  0. 


Combining  2.6  and  the  profile  equations,  2.7,  with  the 
remaining  terns  in  2. 19  results  in  the  siaple  set  of  eguaticns: 


p 

s 

(u*3  /  KZ) 

f  a  (Z/L) 

Po 

^m (Z/L) 

B 

* 

-(u*3  /  fez) 

Z/L  = 

-Po 

Z/L 

6 

= 

Po  [  (Z/L) 

-  Z/L] 

where  Po, 

introduced  in 

section  2.2, 

is 

the  mechanical 

production 

in 

the  equivalent  neutral 

case. 

Thus  u*  can  be 

simply  expressed  as  a  function  of  €  and  Z/L: 

u*’  *  K  Z  €  /  if  n  (Z/L)  -  Z/L  ]  .  (2.20) 

In  a  similar  fashion  the  problem  of  evaluating  r*  and  the 
scalar  fluxes  can  be  simplified  to  finding  u*  and  the 
dissipation  rate  of  scalar  fluctuations.  Nr.  The  analogue  of 
2.19  is  the  simpler  scalar  variance  budget,  (Busch,  1977), 

J  ckr*>  =  -<wr>  3<R>  -  Hr  -  1  d<wr*>  . 

2  dt  dz  2  ~bz  (2. 21) 

The  study  of  Wyngaard  and  Cot€  (1971)  investigates  this  equation 
thoroughly  for  temperature.  The  vertical  divergence  term  turns 
out  to  be  an  order  of  magnitude  smaller  than  the  production  term 
and  the  time  rate  of  change  and  inhomogeneities  are  again 
negligible.  Substituting  for  d<B>/dz  from  2.7  gives  the 


straightforward  relationship 


r**  »  Hr  Z  /  [K  u*  ^r(Z/L)  ] 
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where  Taylor's  hypothesis  is  used  to  replace  the  downstream 
radian  wavenumber  with  (27t£/<0>).  The  fora  of  these  equations 
is  also  based  on  diaensional  analysis  so  the  1-dimensional 
Kolaogoroff  constants  K'  and  Br'  aay  be  functions  of  stability, 
but  they  are  not,  as  yet,  well  enough  established  far  any 
dependency  to  be  observable.  Seasonable  values  are  K'  =  0.55 


and 

Br'  * 

0.80 

for 

both  teaperature  and  moisture 

(Baquin 

and 

Bond,  1971 

and 

Busch, 

1977) ,  with  a  possible  10% 

error. 

In 

ter  as 

of 

the 

natural  frequency,  n»fZ/<0>,  the  -5/3 

region 

has 

been  found  to  be  well  developed  by  n*1  so  that  dissipation 
estiaates  aay  be  obtained  froa  relatively  low  frequencies  (about 
2  Hz,  for  20  a/s  winds  at  10a  height). 

Several  aethods  of  calculating  the  aoaentua  flux  froa 
measure lent s  of  6  and  Z/L  are  feasible.  The  siaplest,  method  1, 
used  by  Denaan  and  Hiyake,  1973,  is  to  eaploy  the  neutral  fora 
of  2.20 
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<uw>DISS1  *  (  X  2  6  )*A»  .  (2.24) 

This  equation  is  also  valid  in  non-neutral  conditions  providing 
there  is  an  overall  balance  between  the  vertical  divergences, 
the  buoyant  production  and  the  stability  modification  of  the 
mechanical  production.  The  experimental  evidence  over  land 
suggests  that  the  complete  form  of  2.20  should  be  tested  over 
the  sea.  Hethod  2,  therefore,  assumes  only  that  the  vertical 
divergences  balance  and  uses 

<uv>DISS2  =  (  K  Z  6  )  &  •  l+n  (Z/L)  -  Z/L)-^  (2.25) 

This  is  the  method  used  by  Khalsa  and  Businger,  1978,  and  by 

Wucknitz,  1976,  but  the  latter  uses  a  Richardson  number 

formulation  and  the  former  assume  a  balance,  following  Syngaard 
and  Cote,  of  buoyancy  and  turbulent  transport  with  a  pressure 
transport  of  order  -Z/L  to  arrive  at  2.25.  It  has  often  been 
assumed  that  local  production,  P,  balances  dissipation  (Smith 
and  Banke,  1975) ,  implying  an  overall  balance  between  the  two 
divergences  and  buoyancy.  This  is  the  assumption  of  method  3, 
which  states 

<uw>DISS3  *  (  K  Z  6  )  &  •  [fm(Z/L)  .  (2.26) 

Pond  et  a 1. ,  1971,  found  that  in  unstable  conditions  the 

momentum  flux  from  the  dissipation  and  Reynolds  flux  methods 
were  in  the  best  agreement  if  they  assumed  that  the  reduction  in 
mechanical  production  due  to  stability  modification  of  the 

I 
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profile  was  coapensated  by  the  net  gain  in  turbulent  energy  froa 
vertical  divergences,  that  is,  €  *  Po  ♦  B.  This  fourth  aethod 
is  expressed  by 

<UV>DISS4  *  {  K  Z  6  )&  •  (1-  Z/L)-*/s.  (2.27) 
The  general  expression  of  the  four  dissipation  aethods  is 


<UW>EISS 

s 

(  X  Z  6  )  •  FI 

Hethod  1 

FI 

*  1 

Method  2 

F2 

*  £*n(Z/L)  -  Z/I 

Hethod  3 

P3 

=  tMz/D  ]-# 

Hethod  4 

F4 

«  £1  -  Z/L  J-*/3 . 

(2.28) 


The  functions,  FT,  are  plotted  in  figure  1  over  the  range  of 
stabilities  likely  to  be  encountered  over  the  sea  at 
aid-’latitudes.  It  is  apparent  that  a  reasonable  aeasure  of 
stability  is  iaportant  to  all  but  aethod  1.  -hny  one  of  the 
aethods  aay  be  valid  over  an  individual  run,  but  over  any 
stability  range  there  should  be  one  dissipation  aethod  that  is 
the  aost  appropriate,  on  average. 

Equations  2.23  and  2.28  indicate  that  <uw>  froa  all  the 
dissipation  aethods  is  proportional  to 


[  K  Z]*/»  K«-»  (Mf)  <0>-^»). 


(2.  29) 


Even  with  no  error  in  the  aeasureaent  of  J*u  (f )  and  <U>,  the 
uncertainties  in  /C,  5%,  in  the  aeasureaent  height  Z,  say  0.5a  in 


10a  and  in  K',  10%,  coaid  coabine  to  produce  a  15%  error  in 

<uw>.  These  errors  and  assuaption  errors  are  likely  to  be 
soaevhat  systeaatic,  bat  fortunately  they  are  not  the  saae  in 
the  Beynolds  flax  aethod.  Intercoaparisons  with  Beynolds  flax 
aeasareaents  are  therefore  essential,  in  order  to  establish  the 
"best**  dissipation  aethod  and  to  ensure  that  there  are  no  aajor 
systeaatic  errors. 

There  are  fewer  ways  of  "juggling"  the  terns  of  2.21  to 
arrive  at  Ht,  and  only  two  aethods  of  calculating  the  sensible 
heat  flux  are  practiced.  At  neutral  stability  and  when  the 
vertical  divergence  tera  is  balanced  by  stability  modification 
of  the  temperature  profile,  aethod  1,  gives 

<wt>DISS1  *  £  K  U*  Ht  Z]l£  .  (2.30) 

Again  the  experimental  evidence  over  land  suggests  using  2.22 
from  which  aethod  2  assumes 

<wt>DISS2  =  £  K  u*  Ht  Z]*fc  •  £  f t  (Z/L)  J- *£.  (2.31) 


The  general  fora  of  calculating  the  sensible  heat  flux  is  siaply 


<wt>DISS 

* 

£  K  u*  Ht  Z 

Method  1 

?1 

»  1.0 

Method  2 

?2 

»  1+tlZ/L) 

(2.32) 


FI  and  F2  are  shown  in  figure  2.  The  aethods  differ 


considerably  even  near  neutral  stability,  therefore,  the  "best" 


I 


FIGOBE  2 


Stability  adjustaents  of  the  two  dissipation 
aethods  of  estiaating  the  sensible  heat  flax. 
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■ethod  should  be  easy  to  establish*  depending  on  the  accuracy  of 
the  stability  aeasureaents,  since  aethod  2  is  sensitive  to 
errors  in  Z/L. 

Equations  2.32*  2.23  and  2.20  indicate  that  <vt>  froa  both 
dissipation  aethods  is  proportional  to 

(  K  [ft  fu]*/a  .  (2.33) 
[Bt»  K*  ]*/* 

Again  uncertainty  in  X*  Z  and  the  Rolaogoroff  constants  together 
could  produce  a  15%  error  in  <vt>DISS.  However*  it  should  be 
possible  to  substantially  reduce  systeaatic  errors  through 
Beynolds  flux  intercoaparisons. 


2.6  Estivating  The  Stability  Paraaeter  Z/L 

In  section  2.2*  the  stability  of  the  air  coluan  is 
characterized  by  a  stability  paraaeter,  Z/L,  which  plays  a 
fundaaental  role  in  the  theory  and  aeasureaent  of  turbulence. 
Unfortunately*  it  is  difficult  to  obtain  with  2.6,  so  three 
aeans  of  estiaating  it  froa  incoaplete  data  will  now  be 
investigated.  The  coaplete  expression  froa  2.6  is* 

2  -  “  KJL  9  <JL_U>  • 

L  u**  To 

Tv  and  To*  the  instantaneous  and  local  average  virtual 
teaperatures*  are  defined  as  the  teaperatures  reguired  to  give 


dry  air  the  sane  density  as  actual  aoist  air  at  the  sane 
pressure.  Lualey  and  Panofsky,  1964,  show  that; 


Tv  *  T  (1  ♦  0.61  a)  , 

where  Tv  and  the  air  tenperature,  T,  are  in  degrees  Kelvin  and  a 
is  specific  huaidity.  Lualey  and  Panofsky  (p  96)  also 
approxiaate  the  virtual  temperature  flux  by 

<w  Tv>  *  <wt>  ♦  0.61  TZ  <wa>  , 

where  TZ  is  the  aean  air  tenperature  at  the  height  Z.  The 
conversion  to  absolute  huaidity,  Q,  in  g/n3,  is  accomplished 
with  (Phelps,  1971), 

Q  a  1298  ( 27 3°K  /  To)  a 
so  Tv  S  T  [1  ♦  To  Q  1.72  X  10-*]  . 

Over  temperate  seas  at  Z  more  than  about  10a  there  is  not  a 
large  tenperature  gradient  and  the  virtual  and  air  temperatures 
differ  by  less  than  2%,  making 

To  =  <Tv>  =  TZ  ♦  TZ*  QZ  1.72  X  10-* 

Z  =  -  g  <w£>  •  [1+  To*  1.72x10-*  <w0>1 

L  U*3  To  <wt>  (2.34) 


reasonable  approximations. 


✓r. 


33 


Over  temperate  seas  the  moisture  content  of  the  atmosphere 
and  the  atmospheric  pressure,  affect  the  air  density  by  about  IX 
and  5%  respectively.  For  dynamic  fluz  calculations,  it  is 
enough  to  only  include  the  atmospheric  pressure.  Pi  in  JtPa,  by 
calculating  the  density  from 


1.29  (273/TZ)  (Pi/101)  . 


(2.35) 


Very  often  the  absolute  humidity  and  moisture  fluz,  <wQ>, 
are  unknown  and  Z/L  must  be  approzimated  from  measurements  of 
u*,  <wt>  and  the  mean  air  and  surface  temperatures.  The  ratio 
of  sensible  heat  fluz  to  latent  heat  fluz,  the  Bowen  ratio,  G, 
can  be  used  to  give 


<2 Q>  * 
<  wt> 


¥8 


0.534  (273/  TZ) 

G 


in  g/a3/°C,  where  the  pressure  and  moisture  effects  on  the 
density  have  been  neglected.  Phelps  and  Pond,  1971,  report  a 
value  of  0.24  for  G  from  their  San  Diego  results.  Substituting 
into  2.34  gives 


1  *  -SJg  <wt>  •  [1  ♦  0.001  To], 

L  u**  To 


which  shows  that  the  moisture  fluz  nay  contribute  about 
one~tbird  as  much  to  the  stability  as  does  the  sensible  heat 
fluz. 
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An  in.  situ  Bowen  ratio  can  be  estimated  from  the  balk 
aerodynamic  parameterizations,  section  2.3,  viz: 

G  -  o  Cp  <wt>  -  P_Cfi  <P>  CT  at  . 

L  <wQ>  L  <0>  CE  AQ 

The  Stanton  number,  CT,  and  Dalton  nuaber,  CE,  have  soaetiaes 
been  found  to  be  nearly  equivalent  (Pond  £t  ^1.  ,  1971),  but 
Francey  and  Garratt,  1978,  find  CT  to  be  30%  lower  than  CE. 
Because  the  huaidity  is  not  the  najor  contributor  to  Z/L  and 
because  of  the  large  error  in  AQ,  CE  and  CT  will  be  assuned 
equal.  The  saturation  huaidity  as  a  function  of  tenperatuxe,  is 
given  by  Hertzaan  e£  aj. ,  1974,  as 

QSAT(T)  =  Cl  ezp  (C2/  T) 

with  Cl  =  6.4038  x  10*  and  C2  *  -5107.4.  This  expression  is  a 
fit  to  a  table  of  saturation  humidities  over  pure  water  at 
various  temperatures  so  that  the  surface  huaidity  over  salt 
water,  QSFC,  is  0.98  QSAT(TSFC).  In  order  to  reduce  systematic 
errors,  a  relative  huaidity  of  75%  will  be  assumed,  which  is  in 
the  middle  of  the  humidity  range  expected  over  temperate  seas. 
Since  huaidity  contributes  only  about  30%  to  Z/L,  this 
assumption  should  introduce  a  random  error  in  Z/L  of  about  20% 
at  worst  and  usually  of  less  than  10%.  The  sea-air  huaidity 
difference,  4Q,  and  hence  a  Bowen  ratio,  can  now  be  found  froa 
TSFC  and  TZ  using 

G  (AT)  S  (p  Cp/L)  AT  [0.98  QSAT(TSFC)  -0.75  QSAT(TZ)]-». 
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Substituting  into  2.34  estimates  Z/L  f r on  u* ,  <wt>,  TZ  and  TSFC: 


£(u*,<wt»  S  -U9  <Si>  •  [1  ♦  lq  2.5x10-0  ] 
t  u*3  lo  G(AT)  (2.36) 

In  cases  where  <wt>  is  also  unknown,  bulk  pa raaete citation 
replaces  it  with  CT  OZ  AT.  The  review  by  Friehe  and  Schaitt, 
1978,  indicates  that  CT  is  about  IxlO-3  in  unstable 
stratification  and  about  0.86x10-3  in  stable.  Therefore,  a 
siapler  estimate  of  Z/L  is  given  by 

Z/L  (u*,aT)  2  -  X_Z  g  CT  OZ  AT  £1  ♦  To  2.5x10-»1. 

U*3  To  G  (AT) 

(2.37) 

Following  Deardorff,  1968,  the  bulk  f  omul  a,  2.9,  replaces 
u*3  with  CD^f*  OZ3  so  that  a  stability  parameter  Z/L (aT)  can  be 
deter ained  soley  froa  the  bulk  parameters,  OZ,  TZ,  and  TSFC.  i 
portion  of  the  Z/L  expression  is  identified  as  a  bulk  Bichar dson 
nuaber 


Hi  (AT)  «  -g  J5 _ AT  (1  ♦  To  2.5x10-M. 

OZ3  To  G  (AT) 

so  Z/L  (AT)  »  X  £2 _  Hi  (AT)  . 

CD  CD  V3 

k  reasonable  average  CD  is  1.25  x  10-3  (Garratt,  1977) ,  so 

Z/L  (AT)  a  11  Hi  (AT)  (CT/  CD)  (2.38) 


with  CT/CD  »  0. 70 

■  0.80 


for  AT  <  0 
AT  >  0 


is  a  jrery  practical  estiaate  of  the  stability  parameter.  This 
differs  'from  Deardorff's  final  fora  of  2/1  *  12  Bi  (AT)  for 
unstable  conditions,  as  it  reflects  more  recent  deter Binations 
of  the  bulk  coefficients.  However,  it  ought  to  be  compared  to 
the  more  exact  expression,  2.36,  whenever  <vt>  and  u*  are  both 


available 
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CHAPTER  3  H5  IHSTROBEHT ATIOK  AMD  EXPEBIBEHTAt  BBOGBAB 

3. 1  Introduction 

In  order  to  collect  the  desired  anoant  of  high  wind  speed 
data  a  Reynolds  flax  system  and  a  dissipation  system  have  been 
designed  for  unattended  operation.  Additional  aspects  of  both 
systems  including  error  analysis,  design  criteria,  and  sensor 
response  requirements  are  given  in  Pond  and  Large,  1978,  and  the 
detailed  analysis  of  the  velocity  measurement  is  also  in  Pond 
et  al..  1979.  The  essential  considerations  are  low  power 
consumption  and  a  large  recording  capability  to  keep  the 
servicing  period  long  and  sensors  able  to  function  in  the  hoped 
for  30-40  m/s  winds  and  accompanying  spray,  ffhen  operating  in  a 
hostile  environment  for  long  periods  of  time,  sensors  and 
electronics  are  likely  to  fail  periodically,  so  whenever 
possible  the  important  measurements  are  either  duplicated  or 
their  sensors  are  calibrated  ia  situ. 

The  results  from  two  field  operations  are  to  be  presented 
in  this  study.  The  first  was  conducted  on  the  Bedford  tower 
near  Halifax  Hova  Scotia,  which  provided  a  stable  enough 
platform  to  allow  meaningful  Reynolds  flux  measurements  to  be 
used  to  "calibrate"  the  dissipation  system.  Intercom parisons 
are  also  possible  with  the  air-sea  interaction  system  from  the 
Bedford  Institute  of  Oceanography,  BIO,  which  was  also  installed 
on  the  tower.  In  a  preliminary  experiment  on  Sable  Island  all 
systems  were  found  to  be  compatible  when  operating  on  the  same 
platform.  The  results  of  that  intercomparison  and  of  a  previous 
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BIO  experiment  on  the  island  are  reported  in  Smith  et  al. ,  1976, 
and  Smith  and  Banke,  1975,  respectively.  The  dissipation 
measurements  were  extended  to  higher  wind  speeds  and  sore  open 
sea  conditions  in  a  second  experiment  conducted  from  the  CCGS 
Quadra  during  its  patrols  at  ocean  weather  station  "PAPA". 


3.2  The  Sensors 

The  velocity  measurements  are  based  on  the  Gill 
propeller-vane  anemometer  (R.H.  Young  Co.) ,  whose  propellers 
are  carefully  constructed  helicoids  that  turn  a  precise  number 
of  revolutions  for  each  meter  of  passing  air  (Baynton,  1S70). 
This  number  was  checked  in  a  wind  tunnel  and  found  to  be  within 
2%  of  the  factory  calibration.  This  accuracy  is  maintained  from 
propeller  to  propeller  and  is  not  affected  by  considerable 
ablation  of  the  leading  edge.  &  problem  does  arise  when  the 
axial  wind  falls  below  about  1  m/s,  because  the  inertia  and 
friction  begin  to  produce  a  non-linear  output.  Another  problem 
is  that  when  the  wind  vector  makes  an  angle  (angle  of  attack) 
greater  than  about  20°,  to  the  propeller  axis  the  apparent  axial 
velocity  component  is  less  than  the  expected  cos  (O')  times  the 
magnitude  of  the  wind,  by  a  factor  /5(-fr)  .  For  angles  of  attack 
between  35  and  75  degrees,  this  non-cosine  behavior  is 
approximated  by  /${*)  *1. 103-0.27  %,  for  in  radians  (Pond  and 
large,  1978).  Although  these  problems  pose  no  serious 
difficulties  in  determining  the  horizontal  velocity  components. 


they  complicate  the  measurement  of  the  vertical  velocity.  B  is 
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derived  froa  a  propeller.  Gill- a,  whose  axis  is  tilted  at  an 
angle,  -<®60o,  to  the  axis  of  a  standard,  Gill-a,  propeller, 
which  is  generally  tilted  at  a  saall  angle  6  froa  the  horizontal 
(figure  3).  At  an  average  wind  speed  greater  than  about  4  a/s 
the  axial  coaponent  of  the  tilted  propeller  always  contains 
enough  of  the  horizontal  wind  to  avoid  its  non-linear  regiae. 
The  propeller  axes  and  instantaneous  wind  vector  are  kept  in 
essentially  the  sane  plane  by  the  vane  and  in  this  way  the 
geoaetry  of  figure  3B  is  always  aaintained  and  corrections  for 
the  non-cosine  behavior  are  possible. 

The  twin  propeller- vane  aneaoneter  of  figure  3  is  described 
fully  in  the  references  cited,  so  only  an  outline  of  how  the 
velocity  coaponents  are  resolved  follows.  The  defined  angles  of 
figure  3C  and  the  following  notation  confora  to  these  previous 
publications.  The  Gill-w  (taking  the  ncn-cosine  behavior  at 
into  its  calibration)  and  Gill-u  signals  froa  the  Beynolds  flux 
systea  supply  the  velocities: 

Gill-u  »  71  *  Q  cos  6  *  w  sin  S 

and  Gill-w  *  V2  =  (Q  cos  sin  K+SH »f  1-0.27  (6-tan-1  (w/Q) )  1 

/3(-0  /?<*) 

with  6  -  tan-1  <w/Q)  ,  (3. 1) 

where  Q  and  w  are  the  horizontal  and  vertical  velocity 
coaponents,  respectively.  Because  72  contains  a  considerable 
contribution  froa  the  horizontal  wind,  the  low  freguency 
variations  in  the  Gill-u  and  Gill-w  signals  track  one  another 
very  closely,  providing  a  check  that  everything  is  working 


FIGURE  3  A:  The  Gill  twin  propeller-vane  aneioieter. 

B:  The  HAT  sensor  housing. 

C:  Definition  of  angles  used  in  resolving  the  velocity 
components  and  calculating  the  tilt  angle  o. 
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properly.  The  tilt  angle  8  needs  to  be  evaluated  before  Q  can 
be  removed  from  the  Gill-w  signal  and  v  calculated.  Since  the 
average  vertical  velocity  must  eventually  go  to  zero,  a  good 
estimate  of  £  is  that  angle  of  rotation  needed  to  make  the 
calculated  <R>  exactly  0,  where  <  >  denotes  the  averaging 

period.  An  average  8  is  measured  over  at  least  15  minutes  to 
give  good  averages  of  the  ratio  I  =  <V1>/<V2>  from  which  it  is 
derived.  Assuming  tan-*(w/Q)  =  w/Q  and  <w>=0, 

<V1>  =  <Q>  cos  8  and 

<V2>  *  <Q>  {cos  {■*♦£)  [  1-0.27  5/^0  ] 

+<wZ>  <Q>-z  sin  {<♦  S)  0.27/ /$(*<)} 

are  correct  to  second  order.  The  term  in  <w*>  <Q>-2  is  only 
about  C.2 %  of  the  previous  term,  leaving 

Y-  {1-0.27  8/  }  [cos«C  -  sin«C  tan  S  ]  , 

which,  assuming  8  =  tan  8,  gives  a  quadratic  in  8.  It  is  the 
negative  square  root  of  the  quadratic  formula  that  is  needed,  so 
8  is  estimated  from 

6  *  b  -  {  b*+  /0. 27)  (T-cos  c< )  /sin*  )  *1* 

b  *  0.5  [cot  *♦  /3(*)/0.27].  (3.2) 

This  expression  shows  that  offset  errors,  which  enter  1,  and 
errors  in <*  and  the /3  relation  produce  apparent  tilts  that  can 
seriously  affect  Reynolds  stress  measurements.  Row  some 


straight  forward  algebra  yields: 


T2-I71  *  B  w  ♦  (1/71)  w* 

1  ■  £0.27/^(<*)  ]  sia  (•<♦&)  cos  &  (3.3) 

B  «  cos  («♦  S)  0.27/^M  ♦  {1-0.27  £/^(*) }  sin «< /cos  £, 

where  Q  has  been  replaced  by  Vl/cos  £  and  a  tern  in  <w*>/<Q>* 
and  teras  of  order  (w/Q)  3  have  been  neglected,  on  a  fixed 
plat fora  the  instantaneous  tilt  actually  depends  on  the  wind 
direction,  bat  the  ase  of  the  "average1*  £  introduces  very  little 
error  (Bond  and  Large,  1978).  The  gnadratic  in  w  is  solved 
using  the  positive  square  root  of  the  gnadratic  for aula,  to  give 
an  estiaate  of  w  for  each  pair  of  recorded  Gill-w  and  Gill-u 
values  in  the  averaging  period.  The  instantaneous  horizontal 
velocity  coaponents  are  then  found  froa: 

Q  »  71/cos  £  -  w  tan  £ 

0  *  Q  cos  (an)  (3.4) 

V  *  v  ■  Q  sin  (an)  , 

where  the  instant an eons  wind  direction,  AH,  equals  <AN>+an  and 
the  nean  direction,  <AH>,  is  chosen  such  that  the  average 
cross-streaa  velocity  <▼>  *  <Q  sin (an)  >  >0.  The  possible  errors 
associated  with  resolving  U,v,w  and  Q  in  this  Banner  are 
suaaarized  in  table  I  (reproduced  froa  fond  and  Large,  1978), 
together  with  their  effects  cn  the  calculated  aoaentua  flux  and 
drag  coefficient.  Soae  errors  should  tend  to  cancel*  so 
hopefully  there  is  no  aore  than  a  *10%  error  in  the  average  CD. 


souses 


1  error  of  ±1' 


| <ow > |  |  CD 


|  Offset  at  5m/s  |  ±2% 
terrors  at  10a/s  j  ±1* 


I  ±2%  in  calibration  |  ± 


relation 


COMMENTS 


|  The  *<  error 
I  believed  to 
|  within  1.5° 

I  2—356  effects 


(Offset  errors  partly! 
|  cancel  rather  than  | 
|  add  giving  about  | 
j 1/2  the  effects  i 


|  6  fluctuations  j Negligible! Negligible! 


|  Non-cosine  |  fi=0  |  ±1* 

Iresponse  of  |6=±10°j  ±10X 

I VI  propeller!  | 


±1*  |Fcr  |S|  <  2-3°  errorj 
±10X  | is  similar  to  £-0.  | 

I  For  | £ |  -5°  it  is  | 

j likely  within  ±5X.  | 


TAELE  I  Summary  of  possible  errors  in  the  velocity 
measurement  and  their  effects  on  the  Beynolds  flux 
method. 


Only  the  dissipation  method  is  practicable  on  a  moving 
platform,  because  the  8  derived  from  T  may  be  very  different 
from  the  instantaneous  tilt  needed  to  find  w.  U  is  also  not 
calculated  from  this  type  of  data  because  it  is  not  a  simple 
matter  to  separate  the  wind  effect  from  the  platform  motion  in 
the  vane's  signal.  The  Gill-u  signal  from  the  dissipation 
system  does  give  an  average  Tl/cos  8  (from  which  a  ship's 
velocity  is  removed  vectorally)  and  approximate  average  values 
of  Q  (and  U,  since  <Q>?1.005  <0>) .  In  this  method  the  velocity 


measurement  is  not  the  major  source  of  error  (section  2.5).  A 
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♦2%  calibration  error  only  leads  to  2.7%  and  -1.3%  errors  in 
|<a«>|  and  CO  respectively  and  at  5s/s  they  are  affected  by  the 
offset  error  by  only  0.3%  and  1.3%.  However*  only  frequencies 
above  those  contaainated  by  the  platfora  notion  aay  be  utilized. 
The  Gill-w  signal  functions  as  a  check  that  the  Gill-u  propeller 
and  electronics  are  working  properly  and,  if  necessary,  as  an 
input  to  the  dissipation  aethod. 

The  enclosure,  HIT,  of  figure  3,  serves  as  a  radiation 
shield  and  offers  protection  against  rain  and  spray  for  the 
teaperature  and  huaidity  sensors  that  are  aonnted  in  it.  Glass 
coated  aicrobead  theraistors  {Victory  Engineering  Corp.)  aeasure 
both  the  aean  and  fluctuating  air  teaperature  vhile  glass  rod 
theraistors  potted  in  epoxy  aeasure  both  the  sea  teaperature  and 
a  aean  air  teaperature.  All  these  transducers  fora  part  of 
siailar  bridge  circuits  whose  non-linearity  balances  that  of  the 
theraistors,  aaking  the  output  of  an  operational  anplifier 
detector  linear  over  about  a  25°C  range.  All  probes  were 
initially  calibrated  in  a  water  bath  against  a  standard  aercury 
theraoaeter  but  later  the  rod  theraistor  provides  an  4a  situ 
calibration  check  of  the  aicrobead.  The  two  teaperature 
aeasureaents  should  not  differ  by  aore  than  0.1°C,  when  both  the 
rod  and  aicrobead  are  working  properly. 

Although  no  latent  heat  flux  data  are  as  yet  available,  a 
brief  description  of  the  atteapted  huaidity  aeasureaents 
follows.  On  Sable  Island  the  huaidity  fluctuations  were  taken 
with  an  aluainioa-oxide  sensor  (Panaaetrics  corp.)  and  a  Brady 
array  (Thunder  Scientific  Corp.),  but  these  failed  because  the 
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sensors  deteriorated  in  the  salt  air  environment.  This  occurred 
less  rapidly  in  the  case  of  the  Brad?  so  an  attempt  was  Bade  on 
the  Bedford  toner  to  provide  an  in  sit a  calibration  by  replacing 
the  aluminium- oxide  probe  with  a  second  Brady  covered  with  a  60 
aicron  stainless  steel  sintered  filter.  The  calibration  drift 
was  auch  reduced,  but  still  serious,  asking  both  Brady  arrays 
unsuitable  for  long  unattended  operation,  for  up  to  a  week  or 
two  the  drift  of  the  filtered  Brady  was  not  too  bad,  but  its 
calibration  was  complicated  by  what  nay  have  been  hysteresis 
effects  and  temperature  sensitivity.  In  addition,  the  response 
of  an  open  Brady  is  aarginal  at  best  (Saith  et  al. .  1976).  It 
therefore  seeaed  best  to  abandon  the  Brady  array  and  to  regard 
all  the  data  froa  it  as  unreliable. 

For  the  ship  operations  where  power  requirements  are  not 
restrictive  a  Lyaan-alpha  hnmidiometer  (Electromagnetic  Research 
Corp.)  has  been  employed  to  give  the  fluctuating  humidity,  while 
a  Cambridge  Systems  (Bodel  2000)  dewpoint  system  provided  the  in 
situ  calibration  and  average  .  The  latter  worked  properly  for 
over  a  month  before  needing  servicing  and  is  proaising.  It  also 
gives  an  aspirated  aean  air  temperature  for  checking  the 
microbead.  The  Lyman-alpha,  however,  required  constant 
attention  as  its  windows  quickly  became  so  dirty  that  its  signal 
went  off  scale  before  providing  any  useful  data. 
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3.3  The  Reynolds  Flos  System 

This  system  includes  the  sensors,  an  electronics  package 
and  6  digital  cassette  tape  recorders  to  sanple,  digitize,  and 
record  the  data  needed  to  detersine  the  turbulent  flaxes  by  the 
Reynolds  flax  sethod,  section  2.9.  It  processes  each  of  its  6 
channels  in  the  sane  nanner  as  shown  in  figure  4,  which 
illustrates  the  data  flow  and  paraaeters  considered  in 
converting  stored  data  back  to  the  original  physical  quantities 
sensed  by  the  transducers  (V1,V2,hB,B,  where  R  is  any  scalar), 
lltbough  the  data  processing  is,  for  convenience,  shown  for  the 
spectra,  it  is  actually  performed  on  the  Fourier  coefficients 
fron  which  the  spectra  and  cospectra  are  derived.  1  Reynolds 
flux  record  consists  of  KG  sequential  groups,  each  forned  by 
saapling  first  the  prewhitener,  Vp,  at  3Hz,  IF  tiaes,  then  the 
law-pass  filter,  YL,  at  1/SSP,  IS  tines.  The  slow  saapling 
period,  SSP,  is  aade  as  long  as  possible  to  conserve  power  and 
tape.  The  high  frequency  variance  and  covariance  lost  because 
the  low-pass  filter  prevents  full  aliasing  is  recovered  by  the 
fast  subsaaples.  The  3  Hz  rate  is  fixed  because  it  is  as  fast 
as  6  channels  can  be  recorded  by  a  pair  of  parallel  recorders, 
each  receiving  3  channels.  HF,HS,HG  and  SSP  are  prograaaed  by 
aeans  of  thuabvheel  switches  and  together  they  detersine  the 
record  length,  the  subsaapling  scheme  and  the  portions  of 
frequency  space  covered  by  each  saapling  rate.  In  additional 
switch  sets  the  tine  interval  between  the  start  of  records.  In 
order  to  prevent  the  cassettes  filling  with  low  wind  speed  data, 
a  wind  speed  liait  can  be  set  by  another  switch  so  that  a 
scheduled  record  is  not  taken  if  the  average  wind  speed  in  the 
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FIGDRE  4:  Signal  processing  in  the  Reynolds  flux  system, 

showing  the  sampling  scheme  and  all  parameters 
considered  in  the  analysis. 
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previous  six  sinates  is  less  than  the  set  li ait.  The  system 
continues  to  collect  data  until  the  three  pairs  of  cassettes  are 
fall. 


The  low  frequency  u,v,w  and  r  Fourier  coefficients  and 
spectra  ♦uvwr(f),  are  available  froa  the  slow  saaples.  The 
low-pass  filters  have  a  1  second  tiae  constant  and  transfer 
function,  HL(f)  (3db  down  at  fL*  0.16  Hz).  SSP  has  always  been 
set  to  3  seconds,  giving  a  Hyguist  frequency  fsn,  (1/6  Hz), 
alaost  equal  to  fl,  so  the  filter  loss  and  aliasing  should 
nearly  cancel  one  another.  The  sensor  transfer  function,  Bs(f), 
is  ignored  because  it  is  typically  3db  down  at  nore  than  5  tines 
fsn.  &  fast  Fourier  transform  of  the  slow  saaples  produces 
spectra,  +L(f),  which  aust  be  nearly  equivalent  to  C*  f2(f)  for 
f<  fsn,  where  the  gain,  C,  includes  the  diaensianal  conversion 
froa  physical  units  to  voltage.  Because  of  the  non-linearity  in 
the  w  calculation,  the  velocity  data  aust  first  be  converted  to 
VI,  V2  and  1H  ■  C-»  (VL-BL-B)  and  then  to  0,  v  and  w  with 
equations  3.3  and  3.4,  before  transforaing  and  obtaining 
the  ^uvw(f)  spectra.  HS  has  always  been  set  to  256,  so  one 
group  of  slows  lasts  12.8  ainutes  and  its  Fourier  cofficients 
occur  at  frequencies  froa  0.0013  Hz  to  fns  and  contain  the 
variance  and  covariance  froa  0.00065  Hz  to  0. 167  Hz. 

Finding  the  high  frequency  cofficients  fron  the  fast 
saaples  is  nore  conplicated  because  the  prewhitener,  Bp(f),  and 
sensor  responses  are  explicitly  involved.  The  prewhitener 
circuits  behave  as  tiae  differentiators  at  low  frequencies  to 
increase  the  signal  levels  and  elininate  spectral  distortion 


froa  non-saapled  frequencies,  then  roll  off  as  B-C  low-pass 
filters  at  the  hiqher  frequencies.  By  design  their  response  is 
aaxiaua  and  nearly  flat  at  the  fast  saapling  Hyguist  frequency, 
1.S  Hz,  so  that  aliasing  occurs  without  auch  loss  of 
contributions  below  a  few  Hz.  The  Fourier  coefficients  of  the 
fast  saaples  (corresponding  to  fp  (f) )  are  easily  converted  to 

+2  If)  *  hif)  [  C  |  Hs  (f )  |  |Hp  (f)  |  y*  . 

Hith  the  velocity  channels,  the  ^2(f)  are  inverse  Fourier 
transforaed  and  a  dc  level  is  added  to  Bake  the  last  fast  saaple 
of  VI,  V2  and  iN  equal  to  the  first  slow  saaple  of  each  group. 
0,  v,  and  w  are  then  calculated  using  the  Bean  coordinates  axes 
as  deterained  by  the  slow  saaples,  and  it  is  the  Fourier 
coefficients  of  these  tine  secies  that  are  used  to  produce  the 
desired  high  frequency  ^uvw(f)  spectra  and  cospectra.  lithout 
the  non-linear  w  calculation  only  one  transforn  need  be 
perforaed  and  this  possibility  is  discussed  in  the  cited 
references.  NF  has  always  been  set  to  128  aaking  the  lowest 
fast  frequency  0.0234  Hz,  so  all  the  variance  and  covariance  of 
the  lower  6  and  part  of  the  seventh  frequency  bands  are  already 
contained  in  the  slow  saaples. 

Flux  estiaates  and  statistical  quantities  are  calculated 
froa  runs  of-HGBP  sequential  groups  of  a  record.  The  total  run 
tiae  is  used  as  the  averaging  period  in  the  6  calculation  and 
coordinate  deteraination.  The  group  spectra  and  cospectra  are 
averaged  together  over  the  run.  in  atteapt  has  been  aade  to  use 
the  overlapping  frequency  bands  to  adjust  the  high  frequency 
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portion,  however,  this  has  not  been  successful  because  the  lower 
6  estiaates  iron  the  fast  sanples  are  not  statistically  very 
certain  and  often  radically  different  fron  one  another  and  fron 
the  other  high  frequency  estiaates.  The  128  low  freguency 

values  span  0.00065  <f<  0.167  Hz,  but  the  eighth  high  freguency 
band  does  not  begin  until  f*  0. 176  Bz  so  an  interaediate  value 
at  f»  0.1715  Bz  (bandwidth  0.009  Bz)  is  foraed  by  interpolating 
between  the  seventh  and  eighth  values. 

The  resulting  186  point  spectra  and  cospectra  are 

integrated  in  three  different  ways.  The  first,  the  II  aethod, 

is  to  siaply  multiply  each  value  by  its  bandwidth  and  sua,  that 

is,  to  integrate  fron  f»  0.00065  Bz.  There  are  two 

« 

disadvantages  with  this  aethod:  first,  the  lowest  natural 
freguency,  n»  fZ/<0>,  included  in  the  integration  decreases  with 
<0>  and  this  could  give  rise  to  apparent  wind  speed 

dependencies;  second,  the  lowest  freguencies  are  not  as 
statistically  certain  as  one  would  like  and  failure  to  converge 
over  a  flux  run  could  create  a  great  deal  of  variability  in 
observations,  h  second  integration  aethod,  12,  alleviates  the 
first  problea  by  including  the  variance  and  covariance  cf  the 
individual  group  aeans  about  the  overall  Bean  of  a  run.  This 
effectively  adds  a  further  BGHP  saaples  at  a  13.5  ainute  period, 
which  contain  the  contributions  fron  0. 00062/HGBP  to  0.00062  Bz. 
Bowever,  the  statistical  uncertainty  in  this  additional 
contribution  is  very  high,  enhancing  the  second  II  problea. 
Bote  that  a  saall  spectral  gap  fron  0.00062  to  0.00065  Bz  is 
present  in  the  12  integration,  because  the  slow  saapling  is 
suspended  during  the  fast  saapling.  The  13  aethod  always  begins 
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its  integration  at  the  sane  natural  frequency,  n=  0.004,  where 
the  spectral  values  are  reasonably  well  established.  The 
integral  is  then  multiplied  by  a  constant  factor  E  to  compensate 
for  the  excluded  low  frequency  contributions.  The  value  of  B 
for  each  quantity  is  found  from  the  normalized  spectrum  or 
cospectrua,  Nf  (n) ,  which  is  determined  in  section  4.3  from 
averages  over  all  available  runs.  The  non-dimensional  (n) , 
plotted  against  log{n)  for  convenience,  should  display  a 
universal  form,  depending  only  on  Z/L.  In  section  4.3,  B  is 
calculated  from  the  normalized  wt  and  uw  cospectra  and  is  found 
to  be  a  function  of  stability.  Although  the  implied  low 
freguency  contribution  may  not  be  exact,  the  results  should  give 
representative  averages  and  be  subject  to  minimal  variability. 
The  three  methods  are  expressed  by: 
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The  13  method  is  the  most  attractive  and  is  compared  with  the 
others  in  section  4.5  to  determine  if  it  is  universally 
applicable. 
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3.4  The  Dissipation  Systea 

This  systea  eaploys  the  saae  sensors  and  preanplifiers,  but 
has  its  ovn  electronics  package  and  two  digital  cassette 
recorders.  It  provides  the  estiaates  of  €  and  Hr  for  the  flux 
calculations  (section  2.5)  and  the  nean  sea  and  air 
teaperatures,  huaidity  and  vind  for  paraaeterization  (section 
2.3)  and  stability  calculations  (section  2.6).  The  architecture 
of  the  data  flow  is  shown  in  figure  5.  Between  recording  at 
intervals  of  Dt,  each  band-pass  filter  output  is  saapled  at  20 
Hz,  digitized,  squared  and  sussed  HI  tiaes.  The  suss  are  stored 
in  internal  nenory  units  until  they  are  written  onto  a  cassette, 
at  which  tiae  all  the  low-pass  filter  channels  are  also  saapled 
and  recorded.  Switches  prograa  HI  to  be  as  large  as  possible 
while  shill  allowing  the  snaaation  to  be  coapleted  before  the 
start  of  a  tape  write.  Dt  is  usually  set  to  4  or  5  ainutes. 
Baking  this  recording  systea  very  econoaical  of  tape. 

The  low-pass  filters  are  single  pole  B-C  circuits  with  a  25 
second  tine  constant,  followed  by  a  unity  gain  operational 
aaplifier,  hence  they  are  well  suited  fcr  providing  the  aean 
velocity  and  scalar  values  over  the  averaging  period  used  for 
the  &  calculation  and  velocity  coaputations  of  section  3. 1. 
These  averaged  results  aay  be  coabined  to  give  the  aeans  over 
longer  tiae  intervals  if  desired. 

The  band-pass  filters  consist  of  double  pole  high  and  low 
pass  stages  centered  noainally  at  fc*  0.4,  0.8  and  1.6  Hz.  The 


prevhitener  is  a  differentiator  at  low  frequencies  that  rolls 
off  as  a  single  pole  B-C  filter  above  10-15  Bz.  The  two  filters 


together  fora  band-pass  filters  represented  by  combined  transfer 
functions  He  (f)  with  center  frequencies  at  about  0.55,  1.05  and 
2.1  Hz.  Since  the  input  spectrum  falls  off  rapidly  (f-*/3)  the 
frequencies  of  greatest  power  output  are  near  the  fc  values. 
The  stored  data  words  divided  by  NX  give  the  average  power,  <P>, 
passing  through  each  band-pass  filter,  k  -5/3  spectrum  across 
He  is  assumed  in  order  to  get  a  discrete  spectral  value  at  each 
fc  for  all  selected  channels  from 

<P>  a  EHO  ^  (f  C) 

where  HHO  -  C*  J  (f/f c)  -#/»  |Hs(f)|>  |  He  (f )  |  *  df 

with  the  integral  taken  over  the  range  of  frequencies  passed  by 
the  filters.  Because  BHO  contains  Hs(f) ,  a  function  of  0 
(section  3.5),  the  integral  must  be  evaluated  at  each  wind  speed 
for  each  band-pass  filter  channel. 

Equation  2.23  relates  ^(fc)  to  the  dissipation  of  scalar 
fluctuations.  Hr,  the  molecular  dissipation,  €,  and  the 
one- dimensional  Kolmogoxoff  constants,  Br'  and  K' .  In  the  case 
of  scalars  ^r(f)  *  f 2(f)  and 

*  <P>  S^3  (2n  /<U>)  */*  fc*/3  (3.6) 
EHO  Br* 

is  available  from  each  scalar  band-pass  filter  once  6  has  been 
determined.  Calculating  £  from  the  velocity  signals  is  more 
complicated  because  72  always  contains  some  w  and  a  non-zero 
tilt  introduces  some  w  into  71  as  well.  Some  error  is  thus 


introduced  into  BHO,  because  at  the  frequencies  utilized  only 
4u(f),  and  not  #w(f),  is  expected  to  be  proportional  to  f-*f*  and 
hence  neither  71  nor  72  should  have  exactly  a  -5/3  spectrun. 
Since  the  spectral  values  of  the  horizontal  velocity  coaponent 
are  nearly  equal  to  those  of  the  dovnstreaa  coaponent,  ^o(f), 
(Pond  and  Large,  1978),  3.1  gives 

f  2  (f )  =  i  fu(f)  S(S,f) 


for  the  71  signal  A  =  cos2  S 

S(S,f)  *  [1+tan2S  ^w(f)/^u(f)  +2tan$  fuw  (f)  /^u(f)  > 

S  (S,f)  should  really  be  placed  inside  the  integral  of  the  BHO 
expression  and  integrated  over  the  band-passed  frequencies,  but 
because  lv(f)/fu(f)  and  tuv(f)/lu(f)  can  only  be  approxiaated, 
they  are  taken  to  be  constants  over  each  filter  and  placed 
outside  the  integral  to  give  an  approxinate  S •  (fi) •  In  the 
inertial  subrange  ♦v(f)/fu(f)  *  4/3,  but  this  behavior  is  not 

observed  near  the  fc*s  and  instead  it  is  taken  to  be  0.81,  1.11 

and  1.29  at  fc  *  0.4,  0.8  and  1.6  Hz  respectively.  These  are 
sinply  the  averages  of  14  sonic  aneaoaeter  observations  froa 
Sable  Island  which  also  show  t uw(fc) /+n(fc)  to  average  -0.11, 
-0.16  and  -0.15.  S'  (5)  is  only  a  saall  correction  (about  1.03 
at  S»-5°  and  0.98  at  6=+5°) ,  but  not  accounting  for  it 
introduces  a  systeaatic  error,  which  on  a  leaning  tower  could 
turn  out  to  be  a  function  of  wind  direction  and  hence  fetch. 
Therefore,  S' (&)  is  applied  to  reduce  the  error  and  sake  it  Bore 
randoa.  The  correction  is  larger  and  the  error  sore  serious  in 
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the  72  case  where  3. 1  gives 

1  «  £  cos  K+6)  ^J/^W  ]* 

S(6,f)  »{1«-tan*  **♦£)  +2tan  &+&)  *uw(f)/^u(f)) 

An  approxiaate  s • (6)  can  again  be  used,  but  ❖  =  S*  <x  aust  be 
assuaed  since  the  data  to  calculate  the  instantaneous  angle  of 
attack  at  each  digitization  is  not  available.  Therefore, 
whenever  possible  it  is  far  aore  desirable  to  calculate  the 
dissipation  froa  the  Gill-u  rather  than  the  Gill-w  data.  Hith 
these  approxiaations  and  reservations  each  velocity  band  pass 
filter  gives: 

e2*3  *  f27r/<0>l»fr  <P>  fc^fr  (3.7) 

A  K*  Bio  S  •  (5) 


3.5  Sensor  Besponse 

When  calculating  €  and  Nr,  sensor  response  corrections  are 
essential  and  they  are  of  soae  iaportance  to  the  Beynolds  flux 
aeasureaents.  Fortunately,  the  dissipation  systen  provides  a 
aeans  of  establishing  Hs  (f )  under  the  actual  turbulent 
conditions  encountered.  Assuaing  the  sensors  behave  as  an  B-C 
filter  (3db  down  at  fo  :  Hs(f)*  (1+jf/fo)”*,  j*/-T  ),  either  Nr 
or  6  values  froa  any  two  band-pass  filters  can  be  Bade  egual  by 
an  appropriate  choice  of  fo.  Of  course  randoa  departures  of  the 
input  spectrua  froa  its  average  -5/3  slope  create  scatter  in 
these  fo's,  but  on  average  nearly  the  saae  response  is  indicated 
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by  all  three  pairs  of  filters  (the  0.4  and  0.6  Bz,  the  0.4  and 
1.6  Hz  and  the  0.8  and  1.6  Bz).  The  success  of  this  technique 
establishes  that  the  assumption  of  an  H-C  response  is 
appropriate  and  that  there  is  a  consistent  average  -5/3  region 
throughout  the  range  of  frsguencies  passed  by  the  three  filter 
cos binations. 

the  response  of  mechanical  sensors,  such  as  the  Sill 
propellers,  is  characterized  by  a  distance  constant,  D  *  the 
wind  passage  required  for  a  63*  recovery  free  a  step  change  in 
velocity.  The  B-C  filter  analogue  gives  27 r  fo«<0>/D.  It  is 
more  important  to  establish  fo  during  low  winds  where  the 
response  is  poor.  In  figure  6,  the  hourly  averaged  wind  speed 
is  plotted  against  the  27rfo  required  to  make  the  Gill-u 
calculations  of  6  from  the  0.8  and  1.6  Bz  filters  egual.  Bow  if 
D  were  truely  a  constant  this  plot  would  be  a  straight  line 
through  the  origin  of  slope  0.  The  data  imply  that  below  12 
m/s,  D  is  about  0.65m,  a  significantly  better  response  than  the 
quoted  value  of  0.8m.  At  higher  speeds  the  response  appears  to 
improve  even  sore  than  expected  (D  decreases)  and  this  is 
incorporated  by  using  the  solid  line  of  figure  61,  D  *  0.56m 
0/ (0-1. On/s) .  The  line  is  a  fit,  for  0>  4m/s,  to  comparisons  of 
all  three  band-pass  filter  outputs  over  the  entire  time  that 
this  particular  anemometer-propeller  (19cm,  2-bladed) 
combination  was  in  use.  She  response  changes  only  slightly  with 
a  different  anemometer,  but  depends  strongly  on  the  type  and 
weight  of  the  propeller.  Beavier  propellers  of  similar 
construction,  used  on  Sable  Island,  give  D  *  1.0m  0/(0-0.7m/s) 
and  19cm  4-bladed  ones  show  0  *  0.79m  0/(0- 1. 8n/s) . 


27TfB  IHZ) 

FIGOBE  6:  Deteraination  of  the  distance  constants  fros  the 

0.6  :  1.6  Hz  band-pass  filter  ratios  of 
i:  the  Gill-a  horizontal  propeller, 

B:  the  Gill— v  tilted  propeller,  59.5°. 
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The  data  of  figures  61  and  6B  coae  froa  the  saae  tiae 
periods,  but  evidently  aore  air  aust  pass  in  the  horizontal 
(longer  distance  constant,  DT)  before  the  tilted  propeller 
responds.  This  nay  reflect  the  non-cosine  behavior,  because 
assuaing  DT*D gives  the  solid  line  of  figure  6B  and 
DT*1.22  0  for°<*60o,  which  is  an  acceptable  fit  for  0<12  a/s. 
The  technique  fails  at  higher  speeds  because  the  w  spectrua 
which  foras  a  large  part  of  the  input  signal,  72,  is  no  longer 
-5/3  throughout  the  Gill-w  0.8  Hz  band-pass  filter.  This  aeans 
that  useful  Gill-w  data  are  less  plentiful  as  well  as  aore 
uncertain  than  that  of  the  Gill-u,  so  it  seeas  preferable  to 
take  the  distance  constant  of  the  tilted  propeller  as  D/ (°<) . 
Such  an  approach  turns  out  to  lie  between  Hicks  (1972  B)  and 
Gill  (1975) ,  whose  results  for  a  60°  angle  of  attack  give 
DT* 1.41  D  and  1.15  D,  respectively. 

1  siailar  procedure  was  used  to  establish  the  in  situ 
response  of  each  aicrobead  theraistor  aounted  in  the  HAT. 
fiegrettably,  it  had  to  be  based  on  very. little  data  because  each 
bead  worked  for  only  a  relatively  short  tiae  and  a  lot  of 
recorded  data  has  been  rejected  because  of  calibration  probleas 
and  suspected  salt  contaaination  of  the  aicrobead.  Host  of  the 
reliable  teaperature  data  froa  the  Bedford  tower  is  included  in 
figure  7,  in  which  the  hourly  averaged  wind  speed  is  plotted 
against  the  sensor  response  required  to  aake  Ht  calculated  froa 
the  0.8  and  1.6  Hz  teaperature  band-pass  filters,  egual. 
Clearly  the  response  iaproves  with  wind  speed.  The  solid  line, 
of  slope  0.90a,  is  an  acceptable  fit  to  all  the  data,  iaplying 
that  the  aicrobead  response,  in  winds  up  to  at  least  20  a/s,  can 


FIGURE  7 


Deteriinati.cn  of  the  aicrobead  theraistor  response 
froa  the  0.8  :  1.6  Hz  band-pass  filter  ratios.  The 
solid  line  has  slope  0.90a. 
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also  be  described  as  an  B-C  filter  with  a  distance  constant,  DB. 
Birake  si  al.  (1970  B)  quote  a  26  Bz  response  for  sisilar 
aicrobeads  at  an  aircraft  spec*.  of  70  s/s  (distance  constant 
about  0.4a).  The  lower  wind  speeds  say  give  the  poorer  response 
indicated  by  figure  7.  It  is  suspected,  however,  that  the 
response  is  liaited  by  the  BIT  itself  and  the  success  of  a 
distance  constant  description  is  a  consequence  of  the  aaount  of 
air  needed  to  flush  the  enclosure,  lccordingly,  the  data  froa 
individual  aicrobeads  agree  with  a  DB  of  0.90a  to  within  ±5%. 

If  the  response  is  being  treated  properly,  the  Beynolds 
stress  derived  froa  each  of  the  three  band-pass  filters  should, 
on  average,  be  egual.  In  figure  8  the  ratics  of  €z t*  (froa  3.7) 
and  equivalently  (froa  equation  2.28)  <uw>  are  plotted  as  a 
function  of  wind  speed  for  the  saae  data  as  used  in  figure  6. 
Evidently  the  chosen  response  is  reasonable.  Very  few 
individual  ratios  differ  froa  1.0  by  aore'than  10X.  At  the 
lower  winds  the  averages  of  all  three  ratios  are  nearly  1. 0  and 
in  the  case  of  the  0.8/1. 6  ratio,  this  is  true  up  to  20  m/s . 
Between  12  and  20  a/s,  the  0.4  Bz  band-pass  filter  appears  to 
have  about  5%  saaller  stress  values  (less  output  power)  than 
expected.  This  aay  reflect  a  deviation  froa  a  purely  B-C 
response,  but  sore  likely  it  is  evidence  of  the  lower  portion  of 
the  filter  lying,  on  average,  below  the  f”*/*  frequencies  during 
these  higher  winds.  At  16  a/s  the  distance  constant  foraulation 
gives  D»  0.6a,  but  if  it  were  taken  to  be  a  constant  0.65a  the 
observed  0.4/0. 8  ratio  would  be  even  less,  because  the  response 
correction  at  0.8  Bz  increases  faster  than  at  0.4  Bz  as  the 
response  is  aade  poorer.  The  correction  is  highly  non-linear 
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PIGUBE  8: 


Ratios  of  <us>  calculated  fro*  different  pairs  of 
hand-pass  filters  as  a  function  of  wind  speed.  The 
solid  line  of  figure  61  gives  the  sensor  response. 
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with  wind  speed  and  very  different  from  filter  to  filter. 
Because  the  ratios  of  figure  8  vary  so  little,  on  average,  with 
■ean  wind  speed,  it  is  unlikely  that  errors  in  the  sensor 
response  will  produce  any  najor  spurious  wind  speed 
dependencies.  &  sinilar  check  on  the  nicrobead  response  is  not 
feasible'  because  both  the  0.4  and  0.8  Hz  band-pass  filters  are 
often  not  entirely  in  the  -5/3  region  of  the  tenperature 
spectrum.  The  inevitable  presence  of  contasinated  data  adds  a 
further  coaplication. 


3.6  The  Experimental  Program 

The  Bedford  tower  experiment  lasted  from  September  1976  to 
April  1977.  A  description  of  the  tower  and  the  results  from  the 
BIO  system  can  be  found  in  Smith,  1979.  The  tower  was  a 
floating  spar  buoy  moored  in  59a  of  water,  which  makes  the  site 
essentially  a  deep  water  wave  regime  (Smith,  1979).  The 
location  and  a  photogragh  of  the  tower  are  shown  in  figure  9. 
The  shortest  fetch  to  the  tower  site  is  from  the  vest  and  is 
about  10  km,  while  open  fetch  conditions  extend  over  a  170° 
range.  The  electronics  packages  can  be  seen  on  the  main  deck, 
about  3a  above  the  sea.  The  Gill  and  BAT  are  at  the  very  top 
alongside  the  BIO  thrust  and  aerovane  anemometers  and  aiczcbead 
thermistor.  The  tide  tables  for  Halifax  harbour  are  used  to 
find  the  phase  and  amplitude  (assumed  equal  to  half  the  tidal 
range)  at  the  beginning  of  each  run,  from  which  the  measurement 
height  Z  is  calculated  assuming  a  purely  H2  tide,  which  is  the 
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predoainant  tide  in  this  area.  The  sea  teaperature  sensor  was 
tied  to  the  tower  about  10a  below  aean  sea  level.  Surface 
aeteorological  observations,  including  the  ataospheric  pressure, 
FA,  for  air  density  calculations,  were  routinely  recorded  at  the 
Shearwater  "Aw  land  station  of  the  Ataospheric  Environaent 
Service,  located  about  15  ha  north  of  the  tower  site  (figure  9). 

While  on  the  tower,  the  Reynolds  flux  systea  recorded  the 
three  aneaoaeter  signals  plus  a  aicrobead  theraistor,  an  open 
Brady  huaidioaeter  and  a  filtered  Brady  froa  the  HAT.  The 
switch  settings  of  section  3.3  were  eaployed  so  that  flux  runs 
of  either  3  or  4  groups  corresponding  to  40.5  or  54  ainutes 
respectively,  can  be  processed.  Froa  September  15  a  new  record 
began  every  hour,  but  on  the  23rd  the  interval  was  increased  to 
3  hours  and  the  lower  wind  speed  liait  was  set  to  10  a/s.  On 
October  4  the  liait  was  changed  to  8  a/s,  then  increased  to 
12a/s  on  October  6  and  put  back  to  lOa/s  on  February  15  where  it 
reaained  until  the  end  of  the  experiaent.  Bean while  the 
dissipation  systea  (figure  5)  sampled  and  recorded  the  low-pass 
data  from  the  sane  six  signals  plus  a  rod  theraistor  froa  the 
HAT  and  the  sea  teaperature.  The  Gill-u,  Gill-w,  aicrobead  and 
open  Brady  signals  were  selected  to  be  band-pass  filtered.  The 
nuaber  of  suaaations,  HI,  was  set  to  4500  to  coaply  with  a  4 
ainute  Dt,  which  enabled  the  pair  of  cassettes  to  last  for  52.5 
days. 
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The  weather  ship  experiment  included  the  four  patrols  of 
CC6S  Quadra  between  July  1977  and  April  1978.  During  the  third 
and  fourth  patrols  a  single  propeller  Gill  anenoseter  was  used 
with  the  Gill-u  signal  also  being  processed  as  the  aissing 
Gill-w  signal.  The  typical  sode  of  operation  of  the  ship  when 
at  "PAPA"  (figure  10)  was  to  drift  with  the  wind  then  return  to 
station  by  steasing  into  the  wind  at  less  than  4  knots.  During 
the  latter  operation  the  great  bulk  of  good  data  were  collected, 
but  sose  useful  data  were  also  collected  as  the  ship  steaaed  at 
7  to  12  knots  while  en  route  to  "PAPA"  (figure  10) .  The 
location  of  the  sensors  on  the  ship's  foreaast  is  shown  in  the 
photograph  of  figure  10.  Cables  were  run  to  the  electronics 
packages  two  decks  below  the  base  of  the  Bast.  Bith  winds 
coaing  over  the  bow  the  aeasured  tilt  angles  are  typically  only 
about  +7°,  indicating  that  the  ship's  distortion  of  the  Bean 
flow  is  not  enough  to  upset  the  dissipation  aethod.  However, 
winds  aore  than  30°  to  starboard  or  60°  to  port  were  found  to 
have  greatly  perturbed  high  frequencies,  which  proved  to  be  very 
unfortunate  because  the  ship  drifted  for  aany  hours  in  such 
winds.  In  addition  the  foreaast  location  received  a  great  deal 
of  spray  so  that  the  nicrobeads  broke  and  the  Ly a an -alpha 
windows  became  dirty  soon  after  the  first  encounter  with  heavy 
seas.  As  a  conseguence  very  little  temperature  data  and  no 
huaidity  data  are  available  froa  the  ship.  The  spray  also 
caused  a  great  deal  of  pitting  in  the  leading  edges  of  the 
propellers,  which  they  can  fortunately  tolerate.  In  a 
subseguent  experiment  (JASIH  1978) ,  the  saae  sensors  were 
aounted  forward  of  the  bow  of  the  7S  Heteor  so  that  the  wind 


FIGURE  10  A:  Ocean  Weather  Station  "PAP A",  50°H,  1«5°f,  and  the 
route  of  the  veatherahips. 

B:  The  instruaentation  on  the  foreaast  of  CCGS  Quadra. 
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carried  the  spray  away  and  a  great  deal  of  both  temperature  &i;d 
humidity  data  were  collected. 

While  on  CCGS  QUADRA,  the  dissipation  system  Ion-pass 
filtered  and  recorded  the  signals  from  the  following  sensors:  a 
rod  and  two  microbead  thermistors  from  the  HAT,  the  Gill 
anemometer,  a  Lyman-alpha  humidiometer  and  the  dewpoint  system. 
The  two  Gill  velocity,  the  two  microbead  and  the  Lyman-alpha 
signals  were  all  band-pass  filtered.  Dt  was  set  at  5  minutes, 
allowing  NI  to  be  5800  and  up  to  56  days  of  data  to  be  stored  on 
the  two  cassettes.  The  Reynolds  flux  system  was  included  to 
provide  velocity  and  scalar  spectra.  It  was  set  up  as  on  the 
tower,  but  the  Brady  arrays  were  replaced  by  the  dewooint  and 
Lyman-alpha  signals.  The  wind  speed  limit  was  always  set  to  at 
least  8  m/s  so  that  flux  records  were  taken  throughout  most  of  a 
patrol.  With  these  switch  settings  it  was  possible  to  turn  the 
systems  on  prior  to  sailing  and  to  retrieve  the  data  cassettes 
upon  return  seven  weeks  later.  A  sea  surface  "bucket" 
temperature  and  the  atmospheric  pressure  are  available  from  the 
ship's  three  hourly  meteorological  observations. 


CHAPTEB  4  BEYN3LDS  F10X  HEASPBEHENTS  FBOB 
THE  BEDFOBD  STABLE  TOWBB 

4. 1  Introduction 

The  Beynolds  flux  data  set  froa  the  Bedford  tower  consists 
of  196  aoaentua  flux  runs  with  winds  up  to  20  a/s  and  the 
aajority  of  stabilities  in  the  range  -0.4  <Z/L<  0.1.  In  aost 
cases  the  fetch  is  unliaited,  but  winds  froa  all  directions, 
except  those  which  put  the  sensors  in  the  wake  of  the  BIO  thrust 
(froa  the  east),  are  allowed,  so  soae  fetches  are  as  short  as 
10  ka.  The  runs  are  restricted  to  the  5  a/s  or  greater  winds 
necessary  to  keep  the  tilted  propeller  in  its  linear  operating 
range,  which  also  ensures  that  the  aeasurenent  height,  Z^13  a, 
is  in  the  "constant  flux"  layer.  It  is  very  gratifying  to  find 
that  in  each  case  where  siaultaneous  dissipation  data  exists, 
192  runs,  the  band-pass  data  confira  the  existence  of  a  -5/3 
region  in  the  velocity  spectra.  Other  runs  are  rejected  because 
the  Gill-w  and  Gill-u  signals  do  not  track  each  other  properly. 

All  196  runs  have  been  considered  for  sensible  heat  flux 
calculations,  because  even  5  a/s  winds  are  sufficient  to  flush 
the  HAT.  However,  during  aany  of  these  runs  the  teaperature 
data  are  not  available,  because  either  the  aicrobead  was  broken 
or  very  cold  air  drove  its  signal  off  scale.  A  few  runs  are 
also  rejected  because  the  Bean  air  teaperatures  froa  the  rod  and 
aicrobead  theraistors  do  not  agree  to  within  ±0.1°C.  In 


addition  aany  aore  runs  have  not  been  processed  due  to  what  is 
believed  to  be  the  sensitivity  of  a  salt  contaminated  aicrobead 
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to  humidity  fluctuations  as  reported  by  Schmitt  et  al. .  1978. 
This  behavior  is  recognized  by  relatively  little  variance  in  the 
temperature  spectrum  belov  n=0.01  and  in  many  of  these  cases, 
but  not  all,  the  dissipation  data  reveal  that  the  temperature 
spectrum  is  not  falling  as  steeply  as  -5/3.  only  60  of  the  196 
runs  have  been  found  to  satisfy  criteria  for  temperature  flux 
calculations. 

All  the  Beynolds  flux  results  are  tabulated  in  the  Appendix 
and  are  referred  to  as  runs  T1  to  T196. 


4.2  The  Stability  Parameter  Z/L 

A  stability  parameter  Z/L  (  at)  is  calculated  from  eguation 
2.38  for  each  of  the  196  runs.  The  surface  temperature,  TSFC, 
is  approximated  by  the  dissipation  system's  sea  temperature 
probe,  TS2A,  which  is  assumed  to  have  fallen  linearly  over  a  16 
hour  period  that  was  not  recorded  (runs  T90-T93) .  The  mean  air 
temperature,  TZ,  at  the  measurement  height,  Z,  usually  comes 
from  the  rod  thermistor  of  the  dissipation  system,  but  over  the 
unrecorded  gap  and  during  periods  when  the  rod's  signal  either 
has  erratic  behavior  (runs  T47-T55,  T74-T85  and  T138-T149)  or  is 
off scale,  the  flux  system's  recordings  of  the  microbead  are 
used.  Unfortunately,  during  runs  T102-T110,  T117-T122  and 
T131-T133  neither  temperature  sensor  was  operational  and  it  is 
necessary  to  use  the  meteorological  observations  from 
Shearwater.  In  figure  11B,  Z/L(aT)  is  plotted  against  the  more 


FZSUBZ  11  Coaparison  of  the  lost  complete  expression  for  the 
stabilit7  paraaeter  Z/L(n*v<vt>)  to: 

1:  an  approxiaate  expression  Z/L(q*,43)> 

B:  the  balk  estiaate  Z/L(AT). 
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exact  Z/l(u*,<wt>)  (equation  2.36)  for  the  60  teaperature  runs. 
The  tvo  calculations  tend  to  agree  on  average  and  seldcn  differ 
by  nore  than  ±0.05,  but  occasionally  the  difference  is 


substantial  (nore  than  0.2).  Since  u*  is  calculated  for  all  the 
runs  it  can  be  used,  instead  of  <0>,  to  calculate  a  stability 
parameter  Z/L (u *, AT)  from  equation  2.37.  However  this  estimate 
does  not  agree  as  well  with  Z/L(u*,<vt>)  as  shown  in  figure  111, 
where  systematic  departures  from  a  1:1  relationship  are  evident. 
Since  <wt>  is  only  available  from  the  temperature  runs  it  often 
must  be  approximated  by  CT  <0>  .AT,  but  evidently  the  associated 
error  is  partially  compensated  by  the  error  in  replacing  u*z 
with  CO  <0>z.  The  bulk  estimate  Z/L  ( AT)  is  to  be  used 
exclusively,  because  it  is  the  best  estimate  of  stability  that 
is  always  available,  even  though  it  may  not  be  very  accurate  for 
an  individual  run. 


4.3  Turbulence  Spectra  Ind  Cospectra 

The  spectra  of  the  fluctuating  velocity  components  and 
fluctuating  temperature,  ^(f) ,  provide  a  means  of  evaluating  the 
performance  of  the  Reynolds  flux  system  (section  3.3). 

To  find  a  value  of  E  for  each  quantity  integrated  by  the  13 
method,  equation  3.5,  normalized  spectra  and  cospectra,  H^(n), 
are  established  by  averaging  over  all  available  runs,  ill 
spectral  values  and  their  corresponding  natural  frequencies, 
n*fZ/<0>  (186  per  run),  are  first  calculated.  The  4 (n)  are  then 
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multiplied  by  f  to  produce  a  variance  preserving  plot  against 
log (n) .  Hext  the  f  f (n)  are  non-dimensionalized  by  dividing  by: 
n*2  in  the  case  of  velocity  spectra  and  cospectra,  (ot) 2  for 
temperature  spectra  and  <vt>  «-Ku*t*  for  w,t  cospectra. 
Finally,  the  discrete  nornalized  spectral  values  from  H  runs  in 
a  particular  stability  range  are  band  averaged  over  n,  such  that 
Alog(n)  remains  constant,  giving  a  mean,  taken  to  be  Hf(n),  and 
a  standard  deviation,  o’,  for  each  band.  A  plot  of  Hf(n) 
vs.  log  (n)  should,  according  to  similarity  theory,  display  a 
universal  fora,  depending  on  stability.  The  normalizing  factors 
u*2»-<uw>,  <tt>«(ot)2  and  <wt>  are  found  for  each  run  by 
integrating  fuw,  ft  and  fwt  from  n*0.004  and  are  therefore  as 
much  as  10%  small,  but  the  error  is  independent  of  wind  speed. 
There  is  a  great  deal  of  scatter  in  the  f(f)'s,  not  only  from 
run  to  run,  but  between  nearby  frequencies  of  the  same  run,  due 
to  the  inherent  variability  of  the  Fourier  coefficients.  The 
latter  effect  is  reduced  by  averaging  over  the  NG  groups  of  a 
flux  run.  However,  it  is  not  reduced  further,  by  averaging  over 
Fourier  bands,  because  this  greatly  increases  the  band-vidth, 
which  should  be  kept  as  narrow  as  possible  in  order  to  keep 
l+(n)  representative  of  its  natural  freguency,  especially  at  low 
frequencies  where  there  are  only  a  few  points  in  each  band.  As 
a  consequence,  or  is  extremely  large  and  not  indicative  ot  the 
variability  of  the  mean,  uf(n).  Since  the  H  runs  are 
independent,  ah*oy/W  is  taken  as  an  estimation  of  the  standard 
deviation  of  the  mean.  This  estimation  assumes  Gaussian 


statistics,  which  should  be  approached  with  a  large  number  of 
runs,  and  that  the  mean  and  or  are  good  measures  of  the 


population  statistics,  which  requires  that  NP,  the  number  of 
points  in  a  band,  be  such  larger  than  M.  This  latter  assunption 
is  not  strictly  satisfied  as  NP  approaches  M,  but  because  each 
run  continues  to  contribute  at  least  one  point,  the  statistics 
should  remain  nearly  Gaussian  with  Gv&O’/'TE  a  reasonable 
approxiaation.  Similarly,  when  NP  becoaes  less  than  M,  each 
point  coaes  from  a  different  run  and  oto-o// NP  is  assumed. 

In  the  following  noraalized  plots  N^  (n)  from  each  band  is 
plotted  in  the  middle  of  the  log(n)  band  and  shown  by  a  diamond 
with  vertical  bars  extending  up  and  down  1  oil.  In  the 
logarithmic  plots  the  means  (squares)  are  plotted  and  solid 
lines  of  -2/3  slope,  indicating  f  (f )  proportional  to  f-5/3,  are 
drawn.  The  runs  are  soaetiaes  split  into  stable  and  unstable 
groups  which  are  averaged  and  plotted  separately.  There  is  no 
atteapt  to  average  over  saaller  stability  ranges  because  the 
large  majority  of  runs  span  only  a  narrow  range  and  because  of 
the  large  uncertainty  in  Z/L(AT). 


Telocity  Spectra 

The  noraalized  spectra  of  the  downstream  velocity 
component,  N^u  (n) ,  are  shown  in  figure  12  and  there  is  a  marked 
dependence  on  stability.  The  peak  of  the  spectrum  of  the  stable 
runs,  figure  12B,  occurs  at  a  natural  freguency  more  than  a 
decade  higher  than  that  of  the  unstable  runs  (figure  12&) ,  whose 
spectrum  in  turn  has  a  greater  proportion  of  its  energy  at  lower 
frequencies.  The  spectral  points  below  n«10-*  cone  from  the 
highest  wind  speed  runs  and  their  average  over  the  natural 
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frequency  band,  log (n) =  -3.21  to  -3.00  is  about  0.75  in  both  the 
stable  and  unstable  case.  It  does  appear,  therefore,  that  a 
spectral  gap  is  emerging  between  the  fluctuating  motion  and  the 
mean  flow.  Some  platform  motion  is  expected,  but  it  is  not 
evident  in  the  spectra,  perhaps  because  it  is  obscured  by 
averaging  over  large  frequency  bands.  The  log- log  plots  of 
figure  12  display  a  similar  shape  and  stability  dependence  as  do 
the  over  land  spectra  of  HcBean,  1971.  In  contrast,  the 
measurements  over  land  of  Kaimal  et  al.,  1972,  suggest  sharper 
peaks. 

Contributions  to  the  bands  between  n=0. 1  and  n=0.4  come 
from  both  the  slow  and  the  fast  samples  depending  on  wind  speed 
and  the  small  pluses  on  figure  12&  represent  the  band  averages 
using  only  the  slow  samples.  These  are  plotted  at  the  average 
value  of  log  (n)  and  not  at  the  band  center  as  are  the  overall 
averages.  In  this  range  (0.1  <n<  0.4)  the  fast  sampled  spectrum 
does,  on  average,  match  the  more  statistically  certain  spectral 
values  found  from  the  slow  samples.  There  is  no  evidence  of 
aliasing  in  the  slow  samples  corroborating  the  argument  that  it 
is  compensated  by  the  1  second  time  constant  low-pass  filter. 
Similar  matching  is  not  done  in  figure  12B  because  the  points 
fall  near  the  peak  of  the  spectrum. 

It  happens  that  all  three  Gill-u  band-pass  filters  may  be 
utilized  under  all  conditions  encountered  at  the  tower  and  at 
"PAPA".  The  logarithmic  plots  of  figure  12  show  that  both  the 
stable  and  unstable  spectra  begin  to  display  a  -5/3  region  at 
about  a  natural  frequency  n*0.2.  ibove  n*1.0,  the  Hyguist 


frequency  when  <0>*18  n/s,  the  sfa  (n)  fall  above  the  -2/3  lines 
(fall  less  rapidly  than  -S/3) ,  reflecting  the  expected  spectral 
distortion  due  to  aliasing.  Figure  8  suggests  that,  at  Z*13a, 
the  output  of  the  0.4  Hz  band-pass  filter  contains  frequencies 
below  the  -5/3  range  with  12  n/s  and  higher  winds,  which  sets 
the  low  frequency  cut-off  of  Bc(f)  at  about  (0.2  12a/s/  13s)  - 
0.2  Hz.  However,  figure  8  also  indicates  that  the  stress 
calculated  fron  the  0.4  Hz  filter  is  only  about  5%  lower,  on 
average,  than  that  calculated  fron  the  other  filters,  for  winds 
up  to  20  s/s.  Therefore,  utilizing  all  three  filters  to  find 
the  stress  should  introduce  a  systematic  error  less  than  2%  and 
in  return  inprove  the  statistical  certainty  of  the  estisate. 
The  effective  low  frequency  cut-off  of  the  0.4  Hz  filter  and 
prewhitener  cosbination  is  n  <0>/Z  3  (0.2  20s/s/  13s)  *0.3  Hz, 
the  3db  down  frequency  of  Hc(f).  The  tower  winds-  are  always 
less  than  20  s/s  and  at  "PAPI",  the  wind  plus  ship  speed  is 
always  below  (0.3Hz  22s/  0.2)  *  33  m/s,  therefore,  the 
dissipation  calculations  always  use  the  data  from  all  three 
band-pass  filters. 

Of  course,  filters  cannot  be  used  in  the  dissipation  sethod 
if  they  pass  frequencies  contaninated  by  the  ship's  notion.  The 
horizontal  velocity  spectrun,  f  ^Q(f),  of  figure  13,  is  an 
average  of  four  Beynolds  flux  records  taken  during  very  rough 
seas  and  slightly  unstable  conditions  at  "P1PA".  The  spectral 
values  are  averaged  over  bands  of  log(f)  in  the  Banner  described 
for  the  noraalized  spectra.  The  nunber  of  runs,  H,  is  4,  so  the 
vertical  bars  extend  tl  ft  *  cr/2 ,  where  o~  is  the  standard 
deviation  about  the  Bean  f  lQ(f)  of  a  log  (f )  band.  The  22  b/s 
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winds  sake  n  and  f  very  nearly  equivalent  and  the  spectral  is, 
accordingly,  nearly  the  sane  as  figure  12 A  except  between  f*0.1 
to  0.3  Hz  where  the  notion  of  CCGS  Quadra  is  conspicuous. 
Thankfully,  even  under  such  extrene  conditions,  very  little  ship 
notion  should  be  passed  by  the  0.4  Hz  band-pass  filter 
prewhitener  combination.  Nonetheless,  its  output  while  cn  a 
ship  is  always  checked. 

Figure  14  shows  the  nornalized  vertical  velocity  spectrum, 
vfw(n).  The  peaks,  near  n  equals  1,  are  barely  reached  with  the 
3Hz  fast  sanpling  rate,  but  there  appears  to  be  a  shift  to  a 
higher  frequency  in  the  stable  case,  which  is  in  accord  with  the 
over  land  studies  previously  cited.  The  spectral  shapes  near 
the  peaks  aust  be  distorted  by  the  aliasing  of  a  relatively 
large  amount  of  high  frequency  variance,  sone  of  which  is 
inevitably  lost  because  there  is  only  partial  correction  for 
sensor  response.  Again  there  is  no  evidence  of  platform  motion. 
At  low  frequencies  ah  is  very  small  and  the  low  spectral  levels 
and  the  observed  stability  dependence  have  also  been  found  in 
other  studies,  indicating  that  the  horizontal  velocity  is  being 
properly  renoved  from  the  tilted  propeller  signal.  In  section 
3.5  the  Gill-w  response  investigation  suggested  that  the  -5/3 
range  of  ^w(f)  begins  above  ns  (0.8Hz  13m/12n/s)  *0.9,  and  this 
is  consistent  with  the  logarithnic  plots  of  figure  14. 
Therefore,  sone  restrictions  have  to  be  inposed  when  the  Gill-w 
band-pass  filter  outputs  are  used  in  the  dissipation  aethod.  On 
the  tower,  for  exanple,  the  0.4,  0.8  and  1.6  Hz  outputs  should 
definitely  not  be  used  when  the  wind  speed  exeeds  (0.3  Hz  13n/ 
1.0)  *  4  a/s,  8  a/s  and  16  n/s,  respectively. 
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FIGURE  14  normalized  vertical  velocity  spectra  from  averages 
of  f  ♦v(n)/u**  over: 

A:  108  unstable  runs  B:  88  stable  runs. 

Vertical  lines  are  ±1  on. 
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The  Temperature  Spectrum 

The  normalized  spectrum  of  the  temperature  fluctuations, 
Hft  (n) ,  from  all  60  temperature  runs,  is  shown  in  figure  15. 
Individual  plots  for  averages  over  the  27  unstable  and  33  stable 
runs  are  not  presented  because  the;  are  not  very  different  and, 
with  so  few  runs,  not  very  statistically  certain.  However,  the 
majority  of  runs  have  ]Z/L|  <  0.1,  so  figure  15  may  be 
representative  of  "near  neutral"  stratification.  Accordingly, 
the  spectral  shape  is  very  much  like  HcBean's  (1971)  near 
neutral  case,  however  the  generalized  temperature  spectrum  of 
Kaimal  et  al. .  1972,  indicates  a  sharper  peak  at  a  higher 
frequency.  An  average  over  the  log(n)  -  -3.25  to  -3.0  frequency 
band  gives  a  mean  of  0.15,  which,  despite  a  large  standard 
deviation,  supports  the  existence  of  a  spectral  gap  that 
seemingly  would  emerge  if  more  high  wind  speed  runs  were 
available.  The  two  pluses  plotted  near  n=0.2  and  n=C.3  are  the 
averages  of  1524  and  402  points,  respectively,  from  the  slow 
samples  only  and  their  average  fits  the  high  frequency  portion 
of  the  temperature  spectrum  quite  well.  The  rise  cf  the  higher 
frequency  plus  may  be  a  consequence  of  combining  the  stable  and 
unstable  runs,  or  due  to  a  little  aliasing,  but  when  the 
microbead  is  suspected  of  responding  to  humidity  fluctuations 
the  peak  of  the  spectrum  is  found  above  n*0.1  and  the  inclusion 
of  only  a  small  amount  of  such  data  could  likely  be  the  source 
of  this  feature. 


f*Z/<U> 


FIGUBE  15  The  normalized  temperature  spectrum  from  averages 
of  f  ♦tfiO/tot)*  over  all  60  temperature  runs 
(pluses  from  slow  samples  only) .  Solid  line  has 
slope  -2/3  and  vertical  lines  extend  ±  1  cm. 
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The  logarithmic  plot  does  not  exhibit  an;  evidence  cf  a 
-5/3  region,  which  the  San  Diego  results  of  Phelps  and  Pond, 
1971,  show  to  begin  at  about  n»0.6.  Here  the  spectres  is 
distorted  by  aliasing  before  it  can  develop,  however  one  of  the 
criteria  for  selecting  the  60  temperature  runs  is  that  the 
ratios  of  the  band  pass  filters  reveal  the  existence  of  a  >5/3 
region  above  n*0.6.  The  3db  down  frequencies  of  the  coabined 


prewhitener  and  0.4,  0.8 

and 

1.6 

Hz  teaperature  band- 

pass 

filters  are  about  0.3, 

0.55 

and 

1.0  Hz  respectively. 

The 

iaplication  is  that  the  use 

of 

the 

teaperature  data  in 

the 

dissipation  sethod  ought  to  be  restricted  in  the  following 
Banner:  the  0.4  Hz  to  winds  less  than  (Z  0.3HZ/0.6)  3  6  n/s  on 
the  tower  and  10  n/s  on  the  ship,  the  0.8  Hz  to  speeds  less  than 
(Z  0. 55HZ/0.6)  ^  12  a/s  on  the  tower  and  20  a/s  on  the  Quadra 
and  the  1.6  Hz  to  <U>  below  (Z  1.0Hz  /0.6)  &  20  a/s  on  the  tower 
and  37  n/s  on  the  ship.  Unluckily,  at  the  higher  wind  speeds 
only  data  froa  the  1.6  Hz  filter  are  useful,  so  there  is  no 
check  on  the  >5/3  region. 

The  u.w  Cospectrum 

The  nornalized  u,w  cospectrua,  ifuw(n),  froa  the  stable 
runs  is  significantly  different  froa  that  found  in  the  unstable 
case,  as  shown  by  figure  16.  The  over  land  results  of  both 
HcBean  and  Hiyake  (1972)  and  Kaiaal  al.  (1972)  are  in 
excellent  agreeaent  with  these  over  sea  spectra.  Vhen 
calculating  2  for  the  13  nethod  of  integrating  u,w  and  w,t 
cospectra  (page  51),  the  H^(n)  are  averaged  over  bands  of 
nlog(n)  *  0.25.  ?or  clarity,  the  bands  in  figures  16  and  17  are 
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FIG DBS  16  Normalized  unstable  (*)  and  stable  (B)  u*w 
cospectra  ±1  ala,  frc»  averages  of  f  tuv(n)/u«. 
Integration  under  the  sclid  curves  gives  E  for  the 
13  aethod. 
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0.40,  bat  the  solid  carves  trace  oat  the  approzisate  areas 

integrated.  Froa  its  peak  at  a  natural  fregaency  n*0.03  the 

unstable  cospectrua,  figure  161,  falls  steeply  to  lover 

fregacncies  leaving  about  1%  of  the  total  covariance  belov 

n»io-J.  the  aore  gentle  fall  to  higher  freguencies  is  usual, 

bat  here,  as  in  figures  12  and  14,  aliasing  distorts  the  shape 

above  n*1.  The  highest  fregaency  point  is  plotted  at  n»4,  the 

band  center,  but  the  average  log  (n)  of  this  band  is  at  u*3. 2, 

which  is  near  the  high  fregaency  cut-off  of  the  integration  over 

the  Alog(n)  *  0.25  bands.  The  total  area  under  the  carve  was 

found  to  be  1.06  tines  the  area  froa  n»0.004.  The  stable 

cospectroa,  figure  16B,  shows  aore  covariance  at  higher 
* 

feguencies  with  a  peak  at  about  n»0. 2.  There  is  no  covariance, 
on  average,  below  n*0.002,  but  individual  runs  often  display 
significant  aaounts,  both  positive  and  negative.  The  total  area 
under  the  solid  curve  is  only  about  1.005  tines  the  area  froa 
n«. 004.  The  ratio  B,  required  for  the  13  aethod  of  integrating 
the  cospectrua,  turns  out  to  be  a  function  cf  stability  and 
Et^l.06  and  Es*1.005  are  realistic  values  for  unstable  and 
stable  cases,  respectively. 

Xb£  JLti  SMSaiM 

The  noraalized  v,t  cospectrua,  tffvt  (n) ,  behaves  in  a 
siailar  fashion  as  Hfuv  and  also  agrees  with  the  over  land 
aeasureaents.  The  unstable  cospectrua,  figure  17&,  displays  a 
broader  peak  at  a  slightly  lover  fregaency  than  does  the  stable, 
figure  17B.  Again  there  is  a  greater  proportion  of  the 
covariance  at  the  lover  freguencies  in  the  unstable  case,  where 
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conditions.  Integration  under  the  solid  curves 
gives  2  for  the  13  nsthod. 


the  total  araa  under  the  solid  carve  is  aboat  1.10  tines  the 
anoant  fron  n»0.004.  In  the  stable  case  this  ratio  is  only 
1.04.  the  suggestion  is  that  the  13  aethod  of  integrating  w,t 
cospectra  ought  to  use  Eu*1. 10  (unstable)  and  2s«1.04  (stable), 
fith  only  33  stable  and  27  unstable  tenperature  runs  these 
cospectra  and  integrations  have  aore  uncertainty,  as  reflected 
by  the  larger  estinates  of  the  standard  deviation  of  the  aean, 
than  found  for  H^u*.  Huaidity  sensitivity  of  the  teaperature 
sensor  is  difficult  to  diagnose  froa  the  cospectrua,  because  it 
results  in  f  fwt(n)/<wt>  being  very  siailar  to  the  average 
stable  case,  figure  17B. 


4.4  Turbulence  Statistics 


The  spectra  and  cospectra  froa  each  run  are  integrated  froa 
f«0. 00065  Hz  (the  II  aethod)  ,  to  give  estiaates  of  the 
statistical  guantities  oh,  oV,  ov ,  art  and  the  u, v  correlation 
coefficient  r(uw)  ■  <uv>/(oru  oV)  .  The  aean,  standard  deviation 
and  wind  speed  dependence  of  sone  of  these  noraalized  guantities 
are  presented  in  table  II,  however  there  aay  be  soae  stability 
dependence  in  these  results.  It  is  evident  froa  figure  12  that 
this  aethod  of  integrating  fu(f)is  likely  to  underestiaate  <ru  by 
aore  than  13*  at  <0>»6  a/s.  This  effect  decreases  with 
increasing  wind  speed  and  accounts  for  at  least  half  of  the 
observed  increase  in  on/<0>  with  <0>.  &  siailar  behavior  is 
expected  for  OV  and  which  also  have  significant  low 
fregueney  contributions.  The  aean  and  scatter  of  OV/<0>  and 
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TABLE  II  Turbulent  velocity  statistics  fion  the  196  Reynolds 
flux  runs.  The  leans  1  1  standard  deviation  fron 
Saith  and  Banke,  1975,  are  shown  for  comparison. 

Oli/<0>  are,  therefore,  quite  reasonable.  The  values  cf  oli/<0> 
and  0W/u*  are  rather  less  than  those  of  Saith  and  Banke,  1975, 
perhaps  because  the  aliased  frequencies  are  not  completely 
corrected  for  sensor  response.  Figure  16  shows  that  the  II 
method  should  never  underestimate  crtiw  by  more  than  5%  and  that 
the  loss  of  high  frequency  covariance  should  be  less  than  half 
the  amount  of  variance  lost  by  ov.  The  combined  error  in  oliw 
may,  therefore,  be  about  half  the  sum  of  the  oil  and  OV  errors, 
which  is  consistent  with  the  correlation  coefficient  r(uw)  being 
similar  to  previously  reported  values. 

Figure  18,  shows  the  means  of  oh/u*,  ov/u*  and  ot/T* 
(T*«Xt*  is  used  to  conform  with  HcBean,  1971)  band  averaged  over 
ranges  of  Z/L.  Only  39  of  the  temperature  runs  are  used  in  the 
averaging,  because  many  runs  are  near  neutral  where  large  ot/T* 
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TABLZ  III  Standard  deviations  of  the  turbulence  statistics 
about  the  stability  band  Beans  plotted  in  figure  18. 

values  occur  as  a  result  of  the  "noise"  in  ot  discussed  by 
HcBean.  lith  so  few  teaperature  runs  available  ot/T*  at 
2/1*0.077  is  plotted  even  though  only  6  runs  fall  in  the  band. 
Otherwise,  the  stability  ranges  are  selected  so  that  at  least  10 
runs  fall  into  each  band.  For  clarity,  sone  standard  deviations 
froa  the  band  averaging  at  snail  vanes  of  I Z/L | ,  are  presented 
in  table  IXI.  Typically,  the  standard  deviations  about  the  aean 
Ct/T*  values  in  figure  18  are  about  0.5.  The  nagnitude  of  the 
standard  deviations  in  table  III  are  conparable  to  the  scatter 
in  HcBean* s  results.  oV/u*  exhibits  the  saallest  stability 
dependence  and  the  least  scatter,  perhaps  because  low  freguency 
contributions  to  0V  are  ainiaal.  ilaost  all  the  averages 
plotted  in  figure  18  fall  within  the  scatter  cf  HcBean*  s  (1971) 


FIGURE  18 


Von-diaensional  turbulence  statistics  band  averaged 
over  stability.  See  table  III  for  standard 
deviations  about  plotted  leans. 
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plots  and  the  observed  differences  are  not  unexpected.  His 
variances  are  coapnted  as  integrals  from  n=0.01  to  10,  giving 
even  saaller  Oil' s,  but  better  tftr's  than  the  II  aethod. 
Accordingly,  his  tfb/u*  values  tend  to  be  lever  and  his  ffV/u*'s 
higher  than  the  corresponding  aeans  in  figure  18.  Hovever,  the 
vind  speed  ranges  of  the  two  studies  are  considerably  different 


(<0>  <8  m/s  for 

all  of 

HcBean' s  runs). 

so 

any  vind 

speed 

dependencies. 

such  as 

shown  in  table 

Hr 

complicate 

the 

comparison.  In 

addi tion. 

both  u*  and  T* 

are 

calculated 

from 

integrals  over  different  frequencies.  Despite  these  problems, 
oV/u*  and  crt/T*  in  figure  18  are  generally  in  excellent 
agreement  vith  HcBean.  Measurements  of  the  sensible  heat  and 
moisture  fluxes  were  used  by  HcBean  to  find  Z/L.  The  overall 
stability  dependence  of  all  the  statistical  quantities  is  very 
similar  in  both  studies,  vhich  lends  credence  to  the  belief  that 
Z/L  (-OT)  is,  on  average,  a  good  estimate  of  the  stability 
parameter. 

There  do  not  appear  to  be  any  unexpected  discrepancies 
between  the  statistical  quantities  and  previous  results.  The 
loss  of  low  frequency  covariance  can  hopefully  be  avoided  vith 
the  12  or  13  methods  of  integrating  cospectra.  The  <J V  values 
are  probably  too  small,  but  because  this  is  due  to  frequencies 
above  1.5  Hz  there  should  be  no  serious  consequences  felt  by 
<uv>  and  <vt>.  In  conclusion  the  sensors  and  Beynolds  flux 
system  seem  to  be  performing  as  expected  and  quite  capable  of 
providing  reliable  estimates  of  the  momentum  and  sensible  heat 


fluxes 
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4.5  The  Flaxes  Of  Honentus  lad  Sensible  Beat 

Three  aethods  of  integrating  cospectra  over  the  uncertain 
low  frequencies  were  discussed  in  section  3.3  and  fornulated  in 
equations  3. 5.  II  includes  all  frequencies  fron  the  slow 
sanples,  12  adds  contributions  fron  the  lower  frequencies 
represented  by  the  group  neans  and  13  integrates  fron  n=0. 004  to 
which  a  stability  dependent  factor,  E,  is  applied  to  account  for 
the  low  frequencies  not  integrated.  The  ratios  of  <uw>  fron  the 
different  nethods  are  averaged  over  2  n/s  wind  speed  bands  and 
tabulated  in  table  IT.  The  13:11  ratio  is  expected  to  converge 
systesatically  to  1.0  fron  greater  values  as  the  wind  speed 
increases.  For  the  unstable  runs  this  should  occur  at  about  13 
n/s  when  the  integral  fron  f*0. 00065  Bz  begins  to  cover  all  the 
natural  frequency  range  of  the  nornalized  cospectrun  (figure 
161).  With  Eu*1.06  this  does  occur,  but  above  14  n/s  the  ratio 
is  again  greater  than  1  possibly  because  these  runs  are  nearer 
to  neutral  conditions  than  the  average  and  require  a  lower  Eu. 
Bowever,  the  overall  average  of  1.00  suggests  that  Eu  should  be 
increased  to  reflect  that  the  II  nethod  sonetines  does  niss  sone 
of  the  covariance.  The  uncertainty  in  Z/L  nakes  stability 
adjustnents  to  Eu  in practical,  but  these  would  affect  the  flux 
by  such  less  than  the  neasurenent  error.  Keeping  Eu»1.06  seens 
reasonable  and  the  overall  12:13  ratio  of  0.99  indicates  that 
using  this  factor  includes  the  average  contribution  to  the 
covariance  fron  the  group  neans,  which  in  fact  usually  reduces 
the  downward  aonentun  flux  as  shown  by  11:12  >1.0. 
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TAB1E  IV  Comparison  of  the  different  methods  of  integrating 
the  u,v  cospectrum. 

At  U=11  m/s  the  II  method  begins  to  integrate  over  the 
entire  stable  normalized  u,v  cospectrum  (figure  16B) .  When 
averaged  over  all  runs  above  11  m/s  the  13:11  ratio  is  1.00  and 
the  12:13  ratio  0.98,  shoving  Es=  1.005  to  be  appropriate.  At 
lover  speeds  the  13  method  does  not  include  lov  freguency 
contributions  vhich  reduce  the  dovnvard  flux,  making  13:11  >1 
and  12:13  <1.  These  positive  contributions  to  -<uv>  vere 


balanced  in  the  normalized  cospectrum  by  negative  contributions 
(dovnvard  flux)  from  the  higher  vind  speed  runs.  One  particular 


run  in  the  6  a/s  band  gives  13:11  *  2.13  and  there  are  not 
enough  other  runs  to  balance  it  off.  Excluding  this  run  froa 
the  band  gives  X3:X1  <1.  Xt  is  likely  then  that  the  overall 
13:11  would  be  reduced  if  aore  runs  were  available  for 
averaging.  One  stable  run  in  the  8  a/s  band  has  such  a  large 
positive  contribution  . froa  the  group  Beans  that  the  total 
integral,  <uv>,  becoaes  positive  with  12:13*  -1.56.  without 
this  one  run  both  the  I2:X3  band  average  and  overall  average 
increase  to  0. 96,  so  that  aost  of  the  covariance  froa  the  group 
aeans  is,  on  average,  included  in  the  X3  aethod.  The  decrease 
of  13:11  with  wind  speed  in  the  stable  case,  suggests  that  if 
sore  runs  and  a  better  aeasure  of  Z/L  were  available  Es  should 
be  aade  a  function  of  stablility. 

The  different  aethods  of  integrating  the  w,t  cospectrua  are 
coa pared  in  table  7.  The  strange  results  in  the  12  a/s, 
unstable  band  are  caused  by  run  Till  whose  group  aean 
contribution,  -0.018  °Ca/s,  is  larger  in  nagnitude  than  and  of 
opposite  sign  to  the  .11  integral,  0.015  °Ca/s.  Xf  this  one  run 
is  excluded,  the  band  average  and  overall  average  I1:X2  becoae  a 
sore  reasonable  1.06  and  0.94,  respectively,  with  the 
corresponding  X2:I3  ratios  going  to  0.99  and  1.11.  For  the  low 
speed  unstable  runs,  X1:X2  less  than  one,  suggests  that  there  is 
a  significant  contribution  to  the  w,t  covariance  at  frequencies 
below  f*0.00065  as  do  the  results  of  BcBean  and  Biyake,  1972. 
Without  high  wind  speed  runs  it  is  not  possible  to  extend  the 
noraalized  cospectrua  (figure  171)  below  n*0. 001,  where  it  is 
still  fairly  large.  The  cospectrua  was  aade  to  drop  off  to  0  at 
n*0.0002  so  that  Eu*1. 10,  which  is  a  coaproaise.  The  BcBean  and 
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TABLE  7  Comparison  of  the  different  methods  of  integrating 
the  wrt  cospectrum. 


Hiyake  results  indicate  a  more  rapid  dropr  but  the  overall  12:13 
ratio  of  1.06  hints  that  Eu  should  be  larger.  The  wind  speeds 
during  unstable  conditions  were  always  too  low  for  the  II  method 
to  integrate  over  the  entire  normalized  cospectrum  making  13:11 
>1  at  all  speeds  and  the  overall  average  1.04. 


At  all  wind  speeds  the  II  method  includes  the  entire  stable 
cospectruo,  figure  17B,  and  accordingly  the  13: II  ratio  shows  no 
systematic  trend  with  wind  speed  and  its  variability  attests  to 
the  existence  of  ran don  low  frequency  contributions  to  the 
sensible  heat  flux  which  are  snoothed  by  the  13  nethod.  It  is 
suspected  that  more  runs  at  the  lower  speeds  would  give  an 
average  ratio  of  1.0,  the  sane  as  above  15  m/s.  Again  there  are 
sone  very  anonalous  runs  with  large  group  Bean  contributions  as 
shown  by  the  12:13  band  averages  of  0.52  and  6.28.  Without  such 
runs  the  13  method  again  accounts  for  most  of  the  average  group 
aean  contribution.  In  view  of  the  uncertainty  in  table  V  and  in 
figure  17B,  caused  by  the  lack  of  runs,  a  value  of  ££=1.04  is 
acceptable. 

Tables  IV  and  V  exhibit  evidence  of  large  random 
contributions  to  the  fluxes  from  the  uncertain  low  frequencies 
causing  a  great  deal  of  scatter  in  the  fluxes  calculated  from 
the  II  and  12  methods.  To  avoid  the  resulting  scatter,  the  13 
method  will  be  adopted  as  the  means  of  integrating  the  cospectra 
of  all  the  runs.  The  choices  for  fuw  of  Eu=1.06  and  Es=1.005 
and  for  0wt  of  Eu=1. 10  and  Es=1.04  have  some  uncertainty,  but 
they  are  reasonable  compromises  and  do  not  appear  to  cause  any 
significant  systematic  errors  or  apparent  trends.  The  integrals 
will  be  denoted  by  <uw>FL0X  and  <wt>FL0X  and  the  derived 
parameters  by  u*FLUX,  t*FL0X  and  so  on. 


The  Reynolds  flax  results  are  presented  in  figures  19  and 
20  and  compared  with  the  dissipation  nethod  in  Chapter  5.  A 
more  complete  view  of  the  behavior  of  the  fluxes  is  offered  in 
Chapter  6  by  the  aore  extensive  dissipation  data  set.  Por  the 
present  the  neasured  Mind,  UZ,  temperatures,  TZ  and  TSFC, 
<uw>FL0X,  <wt>FL0X  and  Z/L(AT)  have  been  put  into  eguations  2.4 
and  2.8  to  give  the  roughness  lengths,  Zo  and  Zot,  and  into  2.11 
and  2.10  to  give  a  vind  speed,  010  an  air  teaperature,  T10,  and 
a  drag  coefficient,  CIO,  at  10  meters.  A  neutral  drag 
coefficient,  CDN,  is  derived  from  Zo  with  equation  2.12.  The 
plot  of  CDN  vs.  010,  figure.  19,  looks  identical  to  the  BIO  tover 
results  in  Smith,  1979,  from  which  a  regression  of  120  near 
neutral  CIO  values  on  010,  gives  0.44  *  0.063  010  =  103C10  as 
compared  to  0.46  ♦  0.069  010  =  103CDN  from  figure  19.  Since 
there  is  nearly  an  egual  partition  between  stable  (triangles) 
and  unstable  (crosses)  runs  in  figure  19,  a  regression  of  CIO  on 
010,  0.43  ♦  .069  010  =  103C10,  is  not  very  different.  The 
higher  coefficients  at  the  higher  vind  speeds  are  commonly 
observed,  but  overall  these  values  are  distinctly  smaller  than 
those  at  similar  wind  speeds  in  Garratf's,  1977,  review,  in 
figure  19  the  stable  (triangles)  and  unstable  (pluses)  data  do 
not  separate  into  distinguishing  patterns.  Throughout  the  small 
stability  range  found  over  the  sea,  average  stability  effects 
appear  to  be  small. 

In  figure  20  <wt>F10Z  is  plotted  against  010 (TSEA-T10)  for 
the  52  temperature  runs  with  (oT|  >  0.5°C,  so  that  a  line  from 
any  point  to  the  origin  has  a  slope  equal  to  CT10,  equation 
2.10.  The  solid  line  represents  the  parameterization  of  Friehe 
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716 QBE  20  <wt>  ws.  010  (TS7C-T 10)  is  °C  and  a/s  for  the  52 

teaperatore  runs  with  I at |  >0.5°C.  Stability 

ranges:  A,  z/L  >0.05;  X,  -0.1  <Z/L<  0.05;  +  , 

-0.2  <Z/1<  -0.1;  Z/L  <-0.2.  Solid  line  is  froa 

Triehe  and  Schaitt  (1976). 
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and  Schaitt,  1976,  which  fits  quite  well  for  -10  <010  aT< 
25°Ca/s.  The  data  do  not  support  an  increase  in  CT10  above  this 
range,  which  Friehe  and  Schaitt  suggest  on  the  basis  of  Saith 
and  Banke's,  1975,  aeasureaents  on  the  beach  at  Sable  Island. 
The  BIO  tower  results  (Saith,  1979)  span  | <0>  aT| <  150°Ca/s  and 
indicate  slightly  higher  CT10  values  than  do  either  Friehe  and 
Schaitt  or  figure  20,  in  both  the  unstable  and  stable  cases. 
Saith  (1979)  finds  10*<wt>=  3.2+  1.10  010  AT,  for  aT>0  and  -0. H 
0.83  010  .AT,  for  AT  <0,  from  regressions  of  <wt>  on  010  at.  ill 
three  studies  show  the  stable  coefficient  to  be  snaller  than 
CT10  in  unstable  stratification.  There  is  not  a  great  deal  of 
scatter  in  figure  20,  except  froa  two  runs  (plotted  as  squares), 
during  which  very  warn  air  noved  over  a  cold  sea.  Such 
conditions  greatly  influence  the  average  sensible  heat  flux  and 
are  discussed  in  section  6.3,  when  the  corresponding  dissipation 
data  is  presented.  Near  neutral  runs  (|Z/L|  <0.05),  plotted  as 
crosses,  are  much  the  sane  as  the  aore  stable  (triangles) ,  but 
the  aore  unstable  (pluses,  Z/L  0.05  to  0.25  and  diaaonds,  0.25 
<Z/I)  seea  to  give  saaller  CTIOs,  however  there  are  far  too  few 
data  to  be  conclusive  as  there  are  wind  speed,  fetch  and  other 
effects  to  consider. 
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FIOX  MS  DISSIPATION  METHODS 

S. 1  Introduction 

the  Bedford  tower  experiment  provides  an  excellent 
opportunity  to  establish,  by  conparison  with  both  the  Reynolds 
flux  and  BIO  results,  a  dissipation  aethod  that  is  valid  over  a 
wide  range  of  open  sea  conditions.  The  dissipation  systes 
recorded  data  during  192  Reynolds  flux  sosentus  runs  and  all  60 
tesperature  runs  and  the  results  are  tabulated  with  those  of 
their  corresponding  flux  runs  in  the  Appendix.  The  stability, 
Z/L (AT)  ,  and  €  give  a  u*DISS  *  (<uw>DISS)  from  each  of  the 
four  nethods  of  nanipulating  the  turbulent  kinetic  energy 
equation  and  velocity  profile,  equations  2.28.  Sisilarly 
<wt>DISS  is  obtained  both  by  including  and  excluding  the 
stability  effect  on  the  tesperature  profile,  equations  2.32. 
The  "best"  sosentus  and  sensible  heat  flux  dissipation  methods 
are  to  be  detersined  by  conparison  with  u*FL0X  and  <vt>FL0X  fros 
the  direct  eddy  correlation  seasuresents.  However,  the 
dissipation  and  Reynolds  flux  systess  are  not  entirely 
independent,  because  they  share  the  saae  sensors.  In  order  to 
cosplete  the  intercosparison  of  the  sethods,  u*DISS  is  also 
checked  with  eddy  correlation  seasuresents  cf  u*  fros  the  BIO 
systes. 
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It  is  essential  to  the  inter cob par ison  that  the  runs  be  as 
nearly  simultaneous  as  possible.  The  dissipation  runs  are 
chosen  to  begin  between  0  and  4  minutes  before  the  start  of  the 
flux  run  and  in  most  cases  the  starts  are  within  2  minutes. 
They  last  for  56  minutes  if  the  flux  run  consists  of  four,  13.5 
minute,  consecutive  groups  and  for  44  minutes  in  the  few  cases 
that  only  three  groups  comprise  the  run.  These  times  become  the 
averaging  periods,  <  >.  The  dissipation,  €,  is  taken  as  the 
average  of  the  three  individual  values  obtained  by  substituting 
the  average  power,  <P>,  from  each  Sill-u  band-pass  filter,  into 
equation  3.7.  The  separate  values  never  differ  from  their 
average  by  more  than  15X  and  a  difference  of  more  than  10X  is 
found  in  only  10  runs.  With  both  (fu(fc)  fc®#5)  and  u*2 
proportional  to  e2#5  (equation  2.23),  their  average  values  differ 
from  those  derived  from  a  single  band-pass  filter  by  less  than 
10X  and  usually  by  less  than  7%.  The  deviation  in  (fu(fc)  fc®#5) 
is  probably  due  to  the  spectrum's  fluctuations  about  a  -5/3  line 
and  the  average  €  should  be  a  good  measure  of  the  molecular 
dissipation.  Values  of  the  dissipation  of  temperature 
fluctuations  are  calculated  from  equation  3.6  using  only  filters 
that  lie  entirely  in  the  -5/3  region  of  ft(f).  Thus,  Nt  is 
often  an  average  of  only  1  or  2  separate  estimates  and  therefore 
may  not  be  as  reliable  as  6.  In  order  for  the  calculated 
(ft  (fc)  fc®#5)  values  to  agree  with  individual  band-pass  filters 
to  within  10X,  the  average,  Nt,  must  not  differ  from  each 
individual  estimate  by  more  than  10%  (equation  2.23).  Equation 
2.33  shows  that  10%  deviations  in  both  ft  (f )  and  ^u(f)  will 
produce  a  10%  deviation  in  the  calculated  <wt>. 
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5.2  The  Hoaentua  Flax 

At  neutral  stability  all  four  noeentua  flax  aethods  are 
equivalent  and  in  non-neutral  conditions  aethods  2,  3  and  4 
siaply  adjust  the  neutral  approxiaation,  aethod  1,  by  a  function 
of  2/1  #  figure  1.  u*DISS1,  equation  2.24,  is  the  siaplest 
calculation  and  it  is  plotted  against  u*P10X  in  figure  21. 
There  is  generally  good  agreeaent  between  the  two  calculations 
froa  the  192  simultaneous  runs,  for  which  2/1  is  usually  between 
-0.45  and  0.20.  A  202  deviation  in  the  u*  estiaations  froa  a 
1:1  relationship  is  indicated  by  the  dashed  lines,  which  satisfy 
the  equation 

l*-yi  /[<*+y)/2]»  0.2  ,  (5.i) 

with  x»u*FI0X  and  y»u*DISS.  In  view  of  the  errors  in  both 
aethods,  deviations  of  this  aagnitude  are  expected.  Points  at 
the  higher  u*  values  rarely  fall  outside  the  dashed  lines.  It 
appears,  therefore,  that  the  neutral  dissipation  aethod  provides 
a  very  good  estiaate  of  aoaentua  fluxes  greater  than  about 
u*a«o. 16  (a/s)2,  which  occur  at  wind  speeds  above  about  11  a/s. 
The  saaller  fluxes  span  a  greater  stability  range  and  the 
cluster  of  points  lying  above  the  upper  dashed  line,  with  u*<0.4 
a/s,  coae  froa  the  aost  stable  runs.  Apparently,  the 
assumptions  of  dissipation  aethod  1,  cause  a  systeaatic  error  in 
stable  conditions,  which  tends  to  sake  u*DISS1  significantly 
greater  than  n*F10X. 
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Coaparison  of  a*  in  a/s  froa  the  neutral 
dissipation  aethod  and  the  Reynolds  flax  aethod  for 
all  192  siaultaneoas  Bedford  tower  runs.  Dashed 
lines  indicate  a  20%  deviation  froa  the  solid  1:1 
line. 
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Hethod  1  should  be  valid  over  a  range  of  "near  neutral" 
stabilities  where  the  buoyancy  and  two  vertical  divergence  terns 
of  the  turbulent  kinetic  energy  equation  and  the  effects  of  Z/L 
on  the  velocity  profile  are  either  snail  or  tend  to  cancel  one 
another.  The  extent  of  such  a  regiae  is  investigated  by 
grouping  the  sinultaneous  runs  according  to  Z/L(^T)  and  seeing 
over  what  range  u*DISS1  and  u*FL(JX  agree  on  average.  On  the 
stable  side  it  is  observed  to  extend  to  Z/L*0.05  as  shown  by 
figure  221.  The  solid  triangles,  representing  the  runs  with 
0.04  <Z/L<  0.05,  are  still  showing  reasonable  agreenent.  For 
Z/L>0. 10  every  run  gives  u*DISS1  greater,  often  by  nore  than 
20%,  than  u*FLUX  and  a  stability  correction  that  reduces  n*DXSS 
is  necessary.  On  the  unstable,  side,  figure  22B  shows  that  down 
to  Z/L  *-0. 10  use  of  the  neutral  equation  is  acceptable.  The 
agreeaent  is  evident  in  the  runs  with  -0.1  <Z/L<  -0.05  (solid 
triangles).  Buns  in  the  -0.3  <Z/L<  -0.1  range  are  not  included 
because  u*DISSl  tends  to  be  snaller  than  u*FL0X  by  an  average  of 
about  10%.  However,  in  non  unstable  conditions  the  two 
calculations  again  tend  to  agree.  1  combination  of  figures  221 
and  22B  shows  that  in  a  "near  neutral"  regiae,  -0.1  <Z/L<  0.05, 
u*DISS1  gives  a  very  good  estiaate  of  the  aoaentua  flux  without 
an  explicit  knowledge  of  the  stability.  Extension  of  the 
unstable  liait  down  to  at  least  Z/Ls~0.4,  introduces  very  little 
error.  This  is  a  useful  result,  because  open  ocean  conditions 
are  often  within  this  range  and  because  an  accurate  Z/L  is  not 
always  available.  It  is  only  in  rather  stable,  Z/L>0.05,  and 
perhaps  in  very  unstable  stratification,  that  the  neutral 
approxiaations  cause  appreciable  errors. 
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Fortunately,  estiaates  of  Z/L (4T)  are  available  for  all  the 
siaoltaneous  runs.  Under  stable  conditions  the  buoyancy  tera  of 
the  turbulent  kinetic  energy  equation  (2.19)  is  a  sink,  so  there 
is  actually  sore  production  than  is  lost  through  dissipation 
alone  and  the  larger  production  gives  larger  u*  estiaates. 
lnother  effect  of  stable  stratification  is  to  reduce  the  aaount 
of  flux  associated  with  a  given  velocity  profile,  so  the  u* 
estiaates  become  smaller.  Method  4  (equation  2.27)  incorporates 
the  buoyancy  tern  and  aethod  3  (2.26)  the  effect  on  the  profile, 
while  aethod  2  (2.25)  does  both,  with  the  profile  effect 
doainating.  Thus,  in  stable  conditions,  u*DISS4>  u*DXSS1> 
u*DISS2>  u*DISS3,  as  is  clearly  shown  by  figure  1.  It  is 
observed  that  in  the  stable  runs  u*DXSS1  tends  to  be  greater 
than  u*F1UX,  so  only  aethods  3  and  2  can  provide  the  proper 
adjustaent.  The  u*  estiaates  froa  both  these  aethods  are 
plotted  against  u*FLUX  for  all  88  stable  runs  in  figure  23. 
Both  dissipation  aethods  appear  to  give  equally  good  agreeaent, 
establishing  that  the  local  production  nearly  equals 
dissipation,  on  average,  with  both  the  divergence  and  buoyancy 
teras  of  equation  2. 19  being  snail.  However,  aethod  2  is 
preferred  because  it  does  account  for  buoyancy,  its  adjustaent 
is  snaller  and  its  assuaption,  that  the  two  divergence  teras 
balance,  is  supported  by  HcBean  and  Elliot's,  1975, 
observations. 
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FIGUBE  23  Coaparisoa  of  88  stable  Beynolds  flux  runs  to 
1:  dissipation  aethod  2  (u*  in  a/s) 

B:  dissipation  aethod  3  (u*  in  a/s). 
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During  unstable  conditions  the  roles  of  buoyancy  and  the 
stability  modified  profile,  on  the  dissipation  estiaates  are 
reversed,  with  u*DISS3>  u*DXSS1>  u*DISS4.  Figure  1  shows  that 
in  method  2  the  effect  on  the  profile  still  dominates  down  to 
Z/L*-0. 1 ,  when  .the  buoyancy  finally  becomes  an  important  source 
of  turbulent  kinetic  energy.  It  about  Z/L*-0.4,  u*DISS2  becomes 
less  than  u*DXSS1.  Throughout  the  unstable  range  of  the  runs, 
u*DISS2  differs  from  u*DISS1  by  less  than  10*  and  gualitatively 
it  is  the  only  method  that  provides  the  appropriate  stability 
refinements.  Hethod  3  clearly  gives  far  too  large  a  u*  in  the 
more  unstable  cases.  The  u*  estiaates  from  methods  4  and  2  are 
plotted  against  u*FLUX  for  the  104  simultaneous  unstable  runs  in 
figure  24.  The  agreement  of  method  2  is  excellent  and  the 
scatter  is  very  evenly  distributed  about  the  1:1  line.  This 
result  indicates  that  the  divergence  terms  also  tend  to  cancel 
in  unstable  conditions  as  suggested  by  BcBean  and  Elliot's 
findings.  Because  u*DISS1  also  gives  a  reasonable  estimate 
throughout  -0.4  <Z/L<  0.05,  the  effect  of  stability  on  the 
profile  must  nearly  balance  the  buoyant  production.  Pond 
Si  41- .  1971,  calculated  u*  up  to  0.3  m/s  using  method  4,  which 
is  seen  to  be  acceptable,  but  not  applicable  to  higher  upvalues. 


sssia  *n 


FIGOEE  24  Conparison  of  104  unstable  Reynolds  flux  runs  to 
k:  dissipation  aethod  2  (u*  in  n/s) 

B:  dissipation  aetbod  4  (u*  in  a/s). 
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Whenever  Z/X  is  available  it  is  evident  that,  ovez  the 
stability  range  of  the  intercoaparison,  method  2  gives  the 
"best**  estimates  of  a*,  henceforth  u*DISS  and  <av>DISS  will  be 
calculated  using  this  method.  Figure  25  is  a  plot  of  u*DISS  vs 
u*FLOX  froa  the  192  siaultaneous  runs  and  a  comparison  with 
figure  21  illustrates  the  overall  effect  of  the  stability 
correction.  In  the  region  0.15  <u*FL0X<  0. 4a/s,  the  reduction 
of  u*DISS  froa  the  aost  stable  runs  greatly  improves  the 
agreement,  but  u*DISS  still  tends  to  be  greater  than  u*FLDX.  is 
a  consequence,  the  average  of  regressions  of  u*DISS  and  u*FL0X, 


u*DISS  »  0.96  U*F10X  ♦  0.025  n/s. 


has  a  positive  offset  and  the  overall  average  u*DISS: u*FLUX 
ratio,  1.03  (standard  deviation  0.10),  is  greater  than  1.00. 
The  two  techniques  differ  by  at  aost  28%  and  usually  by  less 
than  20%  in  u*,  which  is  about  the  aaount  of  scatter  expected, 
i  20%  error  in  <uv>DISS  and  a  non-coapensating  20%  error  in 
<uv >PLUX  give  a  20%  deviation  in  u*  estiaates.  It  is  doubtful 
that  the  u*  agreeaent  would  be  so  good  if  a  major  systematic 
error  had  arisen  froa  the  propeller  response  correction  because 
it  is  of  fundaaental  importance  to  u*DISS,  but  only  secondary  to 
u*PLUX.  Conversely,  the  non-cosine  behavior  of  the  Gill 
propellers  is  a  potential  source  of  substantial  error  in 
Beynolds  flux  measurements,  but  not  to  the  dissipation 
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estiaates,  thus  the  good  agreeaent  indicates  that  it  too  is 
being  treated  properly.  Similarly,  it  appears  as  if  the  low 


frequency  covariances  are  being  handled  satisfactorily. 
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intercoaparison  of  u*  in  a/s  froa  the  "best" 
dissipation  (2)  and  the  "best"  Reynolds  flax  (13) 
aethod  for  all  192  sinultaneous  runs. 
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T1BLE  71  Batio  of  dissipation  to  Reynolds  flax  estiaates  of 
the  aoaentun  flax  band  averaged  over  2  a/s  wind 
speed  intervals. 


The  aoaentaa  flax  and  drag  coefficient  are  proportional  to 
<uw>,  so  trends  in  the  <uw>DISS:<uw>FLOX  ratio  are  investigated 
by  averaging  the  ratio  over  wind  speed  and  stability  bands  in 
tables  71  and  711,  respectively.  In  order  that  the  band  aeans 
are  not  andaly  weighted  by  individual  runs  with  atypical  ratios, 
only  those  runs  whose  ratios  fall  within  ±2  standard  deviations 
about  a  complete  band  average  are  used  to  calculate  the  aeans, 
standard  deviations  and  ranges  in  tables  71  and  711.  In  table 


71,  the  overall  average  of  1.05  and  standard  deviation,  <r,  of 
0.17,  are  quite  acceptable,  i  r  of  17%  is  coaparable  to  the 
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TABLE  VII  Batio  of  dissipation  to  Reynolds  flux  estimates  of 
the  nonentun  flax  band  averaged  over  stability 
ranges. 


expected  error  in  each  nethod  and  the  aean  is  less  than  o/3 
greater  than  the  desired  overall  average  of  1.00,  despite  the 
large  ratios  that  occur  daring  lover  winds  when  stability  has 
its  greatest  range.  Otherwise  there  is  no  systenatic  trend  with 
wind  speed,  although  the  range  of  the  ratio  seeas  to  decrease 
with  increasing  speed.  At  the  higher  wind  speeds  where 
stability  should  always  be  near  neutral  the  band  aeans  are  about 
1.00.  Therefore,  it  seeas  likely  that  any  wind  speed  dependency 
of  the  drag  coefficient  observed  in  dissipation  results  would 
also  be  found  in  corresponding  eddy  correlation  aeasoreaents. 
Table  VII  shows  that,  as  expected,  the  stable  runs  produce  the 
largest  average  ratios.  Figure  1  shows  that  switching  to 
dissipation  nethod  3  would  not  greatly  alter  this  result  (4%  in 
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the  0.10  <Z/L<  0.20  band).  It  is  possible  that  the  higher 
ratios  are  caused  by  underestimating  <uw>FL0X  with  the  13 
■ethod.  is  Z/L  increases,  the  proportion  o f  covariance  at  high 
frequencies  and,  hence,  the  aaount  lost  through  incomplete 
sensor  response  corrections  also  increases,  but  the  effect 
should  not  exceed  5%.  It  has  been  shown  in  section  4.5,  that, 
on  average,  Es*1.005  treats  the  low  frequency  covariance 
adequately,  but  the  four  runs  that  give  such  a  large  average  to 
the  most  stable  band  of  table  Til  occurred  nearly  sequentially 
(runs  T174,  T175,  T176  and  T178),  so  the  band  does  not 
neccessarily  reflect  average  conditions.  Since  there  are  no 
fluctuating  temperature  data  from  these  runs,  Z/L(4l)  could  be 
underestimating  Z/L  and  causing  overestimates  of  u*DISS. 

In  order  to  complete  the  intercomparison  of  methods  it  is 
important  to  show  that  the  observed  agreement  of  figure  25  does 
not  depend  on  using  the  same  sensors.  Dr.  S.D.  Smith  of  the 
Bedford  Institute  of  Oceanography  has  kindly  allowed  same  of  his 
eddy  correlation  measurements  of  u*,  u*BIO,  from  the  tower 
(Smith,  1979)  to  be  compared  to  simultaneous  dissipation 
measurements.  The  runs  overlap  as  do  the  Beynolds  flux  runs  and 
are  typically  about  44  minutes  in  duration.  The  BIO  nark  6.4 
thrust  anemometer  operated  from  October  7  to  December  8,  1976. 
In  figure  26,  u*DXSS  is  plotted  against  u*BX0  for  all  20 
simultaneous  runs  from  this  period.  Os  average,  the  agreement 
is  excellent  and  the  scatter  is  no  larger  than  expected 
considering  the  addition  of  calibration  errors  and  the 
separation  of  the  anemometers  by  about  2m.  I  direct  comparison 
of  u*FL0X  and  u*BI0  values  is  not  possible,  because  simultaneous 
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FIGUBE  26:  Comparison  of  u*  is  m/s  tzom  the  dissipation  system 

and  froa  the  BIO  eddy  correlation  systes  fros  20 
runs  on  the  Bedford  tover. 


117 


neasureaents  ace  race,  however  the  dissipation  intercoapacison 
gives  a  favourable  indirect  one.  These  conclusions  ought  to  be 
qualified  by  noting  that  there  are  disagreeaents  with  the  aean 
winds  froa  later  BIO  data  which  are  still  to  be  resolved.  The 
u*  coaparison  of  these  data  shows  u*DISS  to  be  greater  than 
u*BIO  (often  by  sore  than  20%)  in  28  of  30  runs. 


5.3  The  Sensible  Heat  Flux 

The  sensible  heat  flux  is  estiaated  using  the  two 
dissipation  aethods  for  the  60  siaultaneous  temperature  runs  and 
these  are  compared  to  <wt>FL0X  in  figure  27.  The  dashed  lines 
satisfy  5. 1  with  x=<wt>FL0X  and  y«<vt>DISS,  indicating  a  20% 
deviation  of  the  actual  fluxes  froa  the  solid  1:1  line. 
<wt>DISS,  the  neutral  approxiaation,  is  clearly,  on  average, 
aore  negative  than  <wt>FI0X  (figure  27i) .  Method  2  applies  a 
rather  large  correction  for  the  influence  of  stability  on  the 
teaperature  profile,  but  figure  2  shows  the  adjustaent  to  be  in 
the  proper  sense  in  both  stable  and  unstable  conditions,  k 
regression  of  <wt>DISS2  against  <vt>FLUX  (figure  27B)  gives 

<wt>0ISS2  »  1.04  <vt>FL0I, 

with  a  correlation  coefficient  of  0.99.  The  agreeaent  is 
reaarkably  good  considering  the  errors  in  bgth  aethods  and  the 
sensitivity  of  <wt>DISS2  to  the  rather  uncertain  Z/L(4T).  These 
results  suggest  that  the  turbulent  transport  tera  in  the 
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temperature  variance  budget  is  itself  snail  and  that  it  cannot 
compensate  for  changes  in  the  local  production  of  temperature 
variance  due  to  stability.  In  the  range  -0.03  <z/L<  0.05  the 
correction  to  <vt>DISS1  is  less  than  10%  and,  if  necessary, 
method  1  could  be  used  to  estimate  the  sensible  heat  flux 
without  an  an  explicit  Z/L.  However,  method  2  gives  the  better 
estimate  and  will  always  be  used  to  calculate  <wt>DISS.  The  two 
points  in  the  fourth  guadrants  of  figures  27A  and  27B  illustrate 
an  inherent  difficulty  in  finding  small  fluxes  with  the 
dissipation  method.  <wt>DISS  is  calculated  frcm  equation  2.32 
and  the  sign  of  the  square  root  is  chosen  so  as  to  force  <wt>  to 
have  the  same  sign  as  .aT  *  TSFC-TZ.  However,  it  is  commonly 
observed  in  eddy  correlation  measurements  that  <wt>  is  between  0 
and  +0.05  °Cm/s  when  is  slightly  negative.  There  have  been 
very  few  measurements  of  large  sensible  heat  fluxes  over  the 
sea,  but  the  excellent  agreement  in  figure  27B,  when  the 
magnitude  of  the  flux  is  large,  strongly  suggests  that  this 
measurement  may  be  done  using  the  dissipation  method. 

The  low  and  high  frequency  temperature  fluctuations  affect 
the  sensible  heat  flux  calculations  in  the  same  manner  as  the 
velocity  variations  affect  the  momentum  flux.  Similar  arguments 
regarding  non-cosine  behavior,  sensor  response  and  low  freguengy 
covariance  lead  to  the  conclusion  that  the  temperature  elements 
of  both  the  Beynolds  flux  and  dissipation  systems  must  be 
performing  properly  in  order  to  achieve  the  observed  agreement. 
With  no  large  systematic  errors  in  evidence,  <vt>DXSS  and 
<wt>FLUX  should  both  be  giving  representative  estimates  of  the 
sensible  heat  flux. 
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As  previously  noted  many  temperature  runs  have  not  been 
processed  because  of  suspected  salt  contamination  of  the 
microbead.  In  a  few  of  these  cases  the  temperature  band-pass 
filters  unexpectedly  display  a  -5/3  drop  to  within  10X.  The 
band-pass  filter  check  is,  therefore,  a  necessary,  but 
insufficient  test  for  this  behavior.  This  is  unfortunate 
because  it  restricts  dissipation  temperature  measurements  to 
periods  when  the  temperature  spectrum  can  be  obtained  from  the 
Reynolds  flux  system.  A  possible  means  of  overcoming  this 
problem  would  be  to  average  the  output  of  another  band-pass 
filter  centered  at  about  ns0.02  for  comparison  with  the  outputs 
of  the  filters  in  the  -5/3  range.  It  is  also  interesting  to 
note  that  <wt>DISS  and  <wt>FL0X  are  in  good  agreement  in  some 
unstable  contaminated  runs,  albeit  both  seem  to  be  somewhat 
high. 

It  is  expected  that  the  largest  dissipation  errors  are  in 
the  constants  and  the  assumptions.  The  fact  that  there  are  no 
major  persistent  errors,  suggests  that  the  combined  errors  in 
the  constants  is  not  very  large.  This  is  supported  by  equations 
2.29  and  2.33  which  show  <wt>DISS  to  be  proportional  to 

(  KZ  )  *A*  (Bt ' )  “V*  U*DISS  . 

If,  for  example,  the  Kolmogoroff  constant,  Bt*,  was  greatly  in 
error,  u*DISS  would  be  expected  to  be  more  accurate  than 
<wt>DISS  and  hence,  in  better  agreement  with  Reyn^.  is  flux 
calculations.  The  results,  figures  25  and  27b,  do  not  support 
this.  Similarly,  K  and  Z  appear  to  be  reasonable.  Of  course 
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there  coaid  be  aajor  offsetting  errors,  or  the  constants  nay  be 
highly  variable,  bat  it  seess  lore  likely  that  it  is  the 
assaaptions  which  are  the  aost  ascertain  aspect  of  the 
dissipation  method.  It  is  also  probable  that  aach  of  the 
scatter  in  the  a*  and  <vt>  intercoaparisons  originates  with  the 
Beynolds  flax  values. 

On  average,  the  Beynolds  flax  and  dissipation  aethods  give 
nearly  the  saae  aoaentan  and  sensible  heat  flaxes.  In  addition, 
both  the  dissipation  and  Beynolds  flax  systeas  appear  to  be 
functioning  as  expected.  With  the  possible  exception  of  a*  from 
the  aost  stable  cases  (Z/L>  .05),  the  estiaates  of  <nw>  and 
<vt>,  by  both  techniques,  should  be  reliable  and  free  of  aajor 
systeaatic  errors.  It  is,  therefore,  possible  to  have 
confidence  in  the  dissipation  system  when  it  is  operating  on  a 
ship  where  the  Beynolds  flax  aethod  is  not  practicable. 
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CHAETER  6  DISSIPATION  MEASUREMENTS  FROM  THE  BEDFORD 
STABLE  TOWER  AND  CCGS  QUADRA 

6. 1  Introduction 

The  dissipation  system  has  provided  a  great  deal  more  data 
from  the  Bedford  tower  than  the  Reynolds  flux  system,  because  it 
continuously  records  and  has  proved  to  be  more  reliable.  In 
addition,  a  considerable  amount  of  dissipation  data  has  been 
collected  at  "PAPA".  In  total  1086  hours  of  momentum  flux 
measurements  from  the  tower  and  505  hours  from  the  veathership 
are  found  to  satisfy  a  variety  of  criteria  for  data  reliability. 
Only  237  hours  of  tower  data  and  23  hours  from  CCGS  Quadra  are 
found  to  be  suitable  for  sensible  heat  flux  calculations.  There 
are  several  long  periods  of  continuous  recording  which  make  it 
possible  to  investigate  the  time  histories  of  the  fluxes,  winds 
and  temperatures. 

Only  data  from  the  horizontal,  Gill-u,  propeller  are  used 
to  find  the  dissipation,  €,  and  hence  momentum  flux.  Because  of 
the  greater  uncertainty,  it  is  not  worthwhile  including  the 
moderate  wind  speed  measurements  from  the  tilted,  Gill-w, 
propeller  that  are  available  from  the  few  periods  when  the 
horizontal  propeller  failed.  The  tower  analysis  is  limited  to 
wind  speeds  greater  than  4.0  m/s.  At  lower  winds  the  Gill-u 
propeller  is  still  in  its  linear  regime,  but  the  sensor  response 
corrections  become  very  large  and  it  is  not  impossible  for  the 
thickness  of  the  "constant  flux"  layer  to  fall  below  the 
anemometer  height.  On  CCGS  Quadra  the  measurement  height  is 
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nearly  twice  as  high,  so  the  limit  is  set  at  8.0  s/s.  1  farther 
restriction  requires  the  stability  to  be  in  the  range  -0.6  <Z/L< 
0.15.  At  Z/L*0.15  dissipation  momentum  flax  calculations 
require  a  correction  for  stability  of  about  25X  (figure  1) . 
This  adjustment  is  supported  nainly  by  the  Reynolds  flux 
intercom parison,  where  Z/L  is  greater  than  0.1S  for  only  three 
runs,  and  also  by  BcBean  and  Elliott's,  1S75,  aeasureaents,  that 
include  only  three  stable  runs  with  the  highest  at  Z/l»0.12.  It 
is,  therefore,  dangerous  to  extrapolate  these  findings  beyond 
Z/1*0. 15,  where  the  correction  becoaes  even  larger.  On  the 
unstable  side  there  are  21  runs  in  the  dcBean  and  Elliot  study 
and  104  runs  in  the  section  5. 2  intercoapariscn,  with  the  aost 
unstable  at  Z/L=-0.3  and  -0.45,  respectively.  At  Z/L*-0.6  the 
stability  adjustment  is  still  less  than  10X  so  it  should  be  all 
right  to  extend  dissipation  method  2  at  least  this  far.  This 
criterion  effectively  sets  the  outer  liaits  of  sensible  heat 
flux  calculations,  because  they  require  €  froa  Gill-u.  The 
stability  corrections  to  <wt>DISS  are  1.8  at  Z/L*-0.6  and  0.75 
at  Z/L=0.15  (figure  2)  and  although  these  are  substantial,  they 
are  strongly  supported  by  the  work  of  Wyngaard  and  Cot£,  1971, 
which  consists  of  runs  froa  Z/L*-1.1  to  Z/L*0.1,  and  the 
intercoaparison  of  section  5.3.  Teaperature  data  are  also 
rejected  if  there  is  any  evidence  of  salt  contaaination  of  the 
aicrobead  or  if  the  rod  and  bead  thernistors  do  not  agree,  on 
average,  to  within  ±0. 1°C. 
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In  processing  a  dissipation  momentum  [ temperature ]  run  fron 
either  the  ship  or  tower,  20  ainute  averages  of  all  Gill-u 
[teaperature  aicrobead ]  band-pass  filters  in  the  -5/3  range  of 
fu(f)  [*t(f)]  are  used  to  find  independent  values  of  the 
dissipation  of  turbulent  kinetic  energy  [teaperature 
fluctuations]  from  eguation  3.7  £3.6].  These  individual  values 
are  averaged  together  to  give  the  6  [Ht]  necessary  for  the 
calculation  of  the  ooaentua  flux  [sensible  heat  flux]  fron 
eguation  2.28  [2.32].  However,  if  any  separate  band-pass  value 
differs  from  €  [Ht]  by  Bore  than  15%  [10%]  the  run  is  rejected. 
This  criterion  ensures  that  the  spectral  values  at  the 
freguencies  used  fall  within  10%  of  the  average  -5/3  slope  of 
the  input  spectrum  and  that  fluxes  calculated  fron  individual 
filters  are  within  10%  of  the  value  found  fron  the  average.  The 
20  minute  averages  could  be  used  to  produce  tine  series  of  the 
fluxes  and  related  parameters,  but  in  practice  hourly  averages 
are  employed  because  it  is  felt  the  longer  averaging  tine 
increases  the  reliability  of  the  dissipation  nethod.  In  order 
to  obtain  the  hourly  averages,  three  seguential  20  ninute  flux 
estimates  are  averaged  and  then  used  to  calculate  the  related 
paraneters.  Hereafter,  these  hourly  averages  fron  the 
dissipation  system  will  not  be  specifically  denoted,  but 
referred  to  simply  as  <uw>,  <wt>,  u*  t*,  CD,  CT  and  so  on.  For 
dynamic  flux  calculations  the  air  density  £  is  found  from 
eguation  2.35,  using  the  atmospheric  pressure  from  the 
meteorological  observations  from  Shearwater  and  CCGS  Quadra. 
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BIGOBE  28  The  neutral  drag  coefficient  as  a  function  of  wind 
speed  froa: 

i:  1086  hoarl;  averages  froa  the  Bedford  tower. 

B:  505  hourly  averages  froa  the  CC6S  Quadra. 


6.2  Bulk  Aerodynamic  Parameterization  Of  The  Momentum  Flux 

Hourly  averages  of  u*  and  OZ  and  Z/L  are  used  to  calculate 
010  from  equation  2.  11  and  the  neutral  drag  coefficient ,  CDN, 
from  2.8  and  2.12.  These  are  plotted  separately  for  the  tower 
and  ship  data  in  figure  28.  There  are  very  few  extreme  values 
and  most  points  are  obscured  by  the  concentration  about  a  mean 
CDN  value.  Between  010  *4  to  11  m/s  the  scatter  is  uniform  with 
only  a  few  values  of  103CDN  outside  the  range  0.65  to  1.5.  At 
higher  wind  speeds  a  tendency  for  higher  values  develops  and 
there  appears  to  be  less  scatter.  Above  010=20  m/s  there  are 
few  points,  so  these  data  may  not  fully  reflect  average 
conditions.  Fortunately,  in  the  region  8  <010<  18  m/s,  there  is 
a  very  good  overlap  of  tower  and  ship  results  which  allows  the 
two  situations  to  be  compared. 

In  figure  29,  CDN  is  averaged  over  2  m/s  intervals  of  010, 
provided  there  are  at  least  10  runs  in  a  band,  and  the  means 
plotted  with  vertical  bars  extending  up  and  down  1  standard 
deviation.  For  clarity,  tower  results  (triangles)  are  just  to 
the  left  of  band  center  and  the  ship  points  (pluses)  to  the 
right.  Over  the  range  8.5  <010<  20.5  m/s,  the  band  averaged 
CDN's  from  one  platform  are  seen  to  differ  by  less  than  one 
standard  deviation  (but  usually  less  than  1/2  a)  from  the  other 
platform's  mean.  No  systematic  errors,  such  as  in  the  choice  of 
Z,  appear  to  have  been  introduced  in  moving  the  dissipation 
system  from  the  tower  to  CCGS  Quadra.  The  drag  coefficient 
agreement  also  establishes  that  the  tower  site  is,  in  fact, 
representative  of  open  sea  conditions,  at  least  up  to  20  m/s 
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7Z60BZ  29  Coaparison  of  the  ship  (pluses)  and  tovar 
(trianglas)  nantral  drag  coofficients.  Vortical 
bars  extend  *  1  or  about  plotted  band  averages. 
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winds.  This  allows  the  tower  and  ship  results  to  be  combined 
into  a  single  data  set  of  nearly  1600  hours  from  about  16  months 
of  operation  between  September  1976  and  April  1978. 

£§!<&  Effects 

The  tower  data  of  figure  28A  include  sone  runs  fron  United 
fetches,  whose  inclusion  in  the  overall  data  set  nust  be 
justified.  Op  to  about  10  n/s,  CON  fron  the  tower  does  not  vary 
appreciably  with  wind  speed,  therefore  data  in  this  range  are 
used  to  investigate  the  influence  of  fetch  on  the  drag 
coefficient.  The  results  are  shown  in  table  Till,  where  the 
neutral  drag  coefficient,  CDN,  fron  neasurenents  about  13n  above 


Fetch 

(hn) 


Number 

of 

hours 


Mean  10*CDN 
±1  standard 
deviation 


103CDN 

Bin 


103CDN 

aaz 


Aziauthal 
range 
(°  True) 


10  -  20 


200 


1.14  ±.18 


0.75 


2.03 


253  -  393 


20  -  100 


54 


1.10  ±.22 


0.73 


1.87 


246  -  253 
33  -  65 


100  -  200 


85 


1.13  ±.24 


0.64 


1.76 


235  -  246 


Onliaited 

tower 


263 


1.13  ±.22 


0.62 


1.75 


67  -  235 


Onliaited 
tower ♦ship 


291 


1.14  ±.21 


0.62 


1.75 


Tower 

all  fetches 


590 


_I_ 


1.13  ±.21 


0.62 


2.03 


0  -  360 


TABLE  Till  Tariation  of  the  neutral  drag  coefficient  with 
fetch  for  winds  between  4  and  10  a/s. 


the  sea,  does  not  exhibit  any  dependency  on  fetches  gceatec  than 
10  ka  and,  in  view  of  the  overall  scatter,  the  leans  are  very 
consistent.  This  conclusion  depends  on  using  CDN,  since  the 
aean  CIO  froa  the  10-20  ka  fetch  range  is  17X  larger  than  froa 
the  nnliaited  fetch  runs,  'presnaably  because  the  offshore  Kinds 
tended  to  be  aore  unstable.  The  indication  is,  that  for  winds 
below  10  a/s,  fetch  effects  are  not  iaportant  if  the  aspect 
ratio,  Z/fetch,  is  at  least  as  saall  as  10~3.  There  is  also  an 
indirect  inplication  that  the  surface  roughness,  Zo,  does  not 
depend  -strongly  on  surface  wave  parameters  that  are  not  fully 
developed  by  10  ka  during  10  a/s  winds.  The  data  coapiled  by 
Siegel  (1964,  p  216)  ,  for  fetch  United  waves,  show  the 
significant  wave  height  and  the  phase  speed  froa  a  1000  ka  fetch 
(fetch  g/<U>*“10S)  to  be,  respectively,  10  and  5  tines  greater 
than  expected  with  a  10  ka  fetch.  The  necessity  of  extracting 
the  wind  speed  dependency  froa  CDS  coaplicates  extending  the 
investigation  above  lOa/s,  but  as  the  aspect  ratio  stays  the 
sane,  there  should  not  be  a  sudden  fetch  dependency.  Tower  data 
froa  all  fetches  are,  therefore,  grouped  together  in  all 
subsequent  analyses.  The  better  statistics  resulting  froa  the 
additional  data  should  aore  than  coapensate  for  any  saall  fetch 
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Effects 


Equation  2. 12  defines  CDN  as  a  function  only  of  2o  and  2.8 
iaplies  CIOCDN  for  2/L<0,  C10<CDN  for  2/L>0  and  C10=CDN  at 
neutral  stability.  Table  IX  indicates  that  CDN  and  hence  2o  are 
independent  of  2/L  in  the  unstable  case.  it  2/L<-0. 3,  CIO 
values  are  clearly  higher  than  the  average,  but  reduction  to 
neutral  stability  gives  no  systeaatic  trend  to  the  aean  unstable 
CDN's,  which  are  reasonably  consistent  in  view  of  the  scatter. 
Extraneous  effects  sees  to  influence  the  stable  runs  to  a  auch 
greater  degree.  Nhen  reduced  to  neutral,  the  average  stable 
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The  stability  dependency  of  the  drag  coefficients 
froa  the  4  <  010  <  10  a/s  data  of  table  Till. 
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CDS  *s  are  larger  than  the  overall  mean  and  they  seen  to  increase 
with  stability.  However,  the  large  drag  coefficients 
responsible  for  this  feature  are  found  to  be  associated  with 
rapid  changes  in  wind  direction  following  a  stors.  Such  a 
sequence  of  events  say  possibly  influence  parameters  inportant 
to  Zo  and  aerely  coincide  with  stable  stratification.  A  tiae 
series  shoving  the  development  of  these  high,  low  wind  speed, 
stable  drag  coefficients  is  presented  later  (figure  31).  The 
data  suggest  that  the  stability  dependence  of  the  drag 
coefficient  is  well  described  by  siailarity  theory  with  Zo 
independent  of  Z/L  and  that  sea  surface  conditions  are  the 
source  of  the  trend  in  CDH  with  Z/L  found  in  table  IX. 

aisfl  Dependency 

In  figure  30  the  1591  neutral  coefficients  froa  all  the 
data  are  band  averaged  over  2  a/s  intervals  and  plotted  at  the 
010  mean  (vertical  bars  show  t  1  or).  The  Bean  CIO  values  are 
siailar,  because  of  the  almost  equal  numbers  of  stable  and 
unstable  runs.  1  saooth  curve  (concave  up)  with  CDH  increasing 
with  010  could  be  Bade  to  fit  these  points.  However,  below 
about  lOa/s,  the  slight  increase  of  the  average  CDN's  should  be 
adequately  described  by  a  constant.  The  higher  wind  speed 
points  show  a  aore  rapid  rise  of  CDH  with  010,  that  is 
approxiaately  linear  from  10  to  26  a/s.  In  order  to  quantify 
this  behavior,  all  618  hourly  CDH  values  with  010<10  a/s  are 
averaged  and  all  973  with  U1Q>10  a/s  are  regressed  against  010 
(correlation  coefficient  0.74).  All  the  band  averages  are  very 
well  described  by  the  resulting  solid  line  of  figure  30: 
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PIGOHE  30  The  neutral  drag  coefficient  averaged  over  vind 
speed  bands.  Vertical  bars  indicate  ±  1  o'  and  the 
nunber  of  points  in  band  is  shown  below  each 
average.  Lines  show  the  Chaznock  representation 
with  °< *  0.0144,  0.41  (dashed)  and  egnations  6.1 

(solid) . 
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10*  CON  -  10*  CIO  *  1.1*  4  <  010  <  10  m/s 

(6.1) 

10*  CDN  -  0.49  ♦  .065  010  10  <  010  <  27  s/s 

10*  CIO  *  0.46  ♦  .068  010  , 

where  the  CIO  results  are  included  for  comparison,  the  exact 
numerical  values  of  6.1  depend  slightly  on  the  somewhat 
arbitrary  choice  of  a  lower  wind  speed  limit  for  the  regression. 
Conveniently,  with  the  choice  of  10  m/s,  the  line  segments 
happen  to  match.  The  behavior  of  CDH  described  by  6.1  is 
distinctly  different  from  the  dashed  curve  of  figure  30,  which 
is  the  Charnock  representation  with  •<  *  0.0144,  0.4  1  (or 

•0*0.016,  K*  0.4),  as  suggested  by  sarratt  (1977).  Equations 
6.1  fit  Smith's  (1979)  eddy  correlation  data  from  the  Bedford 
tower  site  extremely  w.ell.  For  winds  below  10  m/s,  his  average 
10*CDH  is  1.11  from  14  runs  (1976-1978)  and  about  1.24  from  11 
near  neutral  runs  (1968-1969).  1  regression  of  all  120  runs 

(1976-1978)  between  6  and  22  m/s  yields  0.44  ♦  0.063  010  » 

10*C  10.  However,  the  Sable  Island  results  of  Smith  and  Banke 
(1975)  are  significantly  higher,  suggesting  that  the  Sable 
Island  site  with  its  surf  zone  and  perhaps  other  shallow  water 
locations:  Lough  Beagh  at  8  to  15  meters  (Sheppard  &£  a)... 
1972),  Lake  Flevo  at  4  meters  (Hieringa,  1974),  the  10m  Caspian 
Sea  location  of  Kitaigorodski  $£  al.  (1973)  and  the  Spanish 
Banks  site  of  Hiyake  si  (1970  1)  and  Heiler  and  Burling 

(1967) ,  may  not  be  representative  of.  the  open  ocean  situation. 
Unfortunately,  almost  all  the  measurements  above  14  a/s  (and 
aany  at  lover  speeds) ,  that  were  available  tc  Garratt  case  from 
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these  shallow  water  locations.  A  number  of  open  sea 
measurements  contributed  to  his  10  to  11  m/s  band  average,  which 
very  nearly  lies  on  6. 1.  It  also  seems  likely  that  excluding 
the  results  of  Sheppard  et  al.  would  give  band  averages  in  the  3 
to  7  m/s  .range  that  would  be  adequately  described  by  6.1.  In 
the  range  7  to  10  m/s,  the  band  averages  of  Garratt  generally 
run  about  20%  (less  than  1  O' )  higher  than  equations  6.1. 

Open  ocean  measurements  used  by  Garratt  (1977)  are  in 
better  agreement.  From  the  Argus  Island  tower,  a  site  similar 
to  the  Bedford  tower,  the  overall  average  of  69,  near  neutral, 
Reynolds  flux  CD's  at  7.5m  is  1.24x10-3  for  wind  speeds  from  4 
to  10  m/s  (De  Leonibus,  1971).  The  corresponding  value  of  CIO 
is  about  1.17x10-3.  Brocks  and  Krugermeyer  (1970)  present  CIO's 
from  152,  near  neutral,  15  minute  profiles  over  the  Baltic  and 
North  Seas  and  the  Equatorial  Atlantic  measurements  (787,  near 
neutral,  10  minute  profiles)  of  Hoeber  (1969).  Eguations  6.1 
are  within  the  scatter,  but  generally  lower  by  10  to  15S.  The 
overall  average  103ci0  from  Hoeber  (5  to  11  m/s)  is  1.23±0.25 
and  his  band  averages  agree  very  well  with  those  cf  figure  30, 
except  between  5  and  6  m/s  where  an  average  of  only  11  points  is 
1.48.  With  2/3  of  his  profiles  unstable,  CDN  may  average  a  few 
percent  less.  Averaging  the  Baltic  and  North  Sea  data  (4  to  12 
m/s)  gives  103cl0  *  1.30±0. 18,  but  there  is  a  slight  trend  with 
wind  speed.  The  Bass  Strait  site  of  Hicks  (1972  A)  ought  to  be 
representative  of  the  open  sea,  however  his  wind  speeds  are 
reduced  by  a*  to  account  for  surface  drift.  If  the  eddy 
correlation  drag  coefficients  are  reduced  by  about  7%,  they  then 
conform  to  the  CDN's  of  figure  28  and  the  30  values  between  3 
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and  8  a/s  than  average  1.13x10-*,  however  6  runs  frca  8  to  10 
a/s  average  about  1.4x10-*.  Garratt  also  uses  the  18  eddy 
correlation  runs  (3  to  11  a/s)  £roa  Hasse  (1970),  quoting  an 
average  of  10*C10  =  1.21  ±20%. 

The  two  data  sets  used  by  Garratt  that  reaaln  to  be 
discussed  are  open  sea  studies  froa  5/7  FLIP,  however,  they  are 
not  independent.  Because  of  flow  distortion  due  to  FLIP, 
Paulson  et  a],.  (1972)  use  an  eapirical  correction  factor  chosen 
so  that  their  average  <uw>  (froa  profiles)  and  the  average  of 
corresponding  eddy  correlation  aeasurenents  froa  Pond 
et  al.  (1971) ,  are  equal.  The  latter  data  (20  runs)  give  an 
average  CD  at  about  8  aeters  of  1.52x10-*  in  winds  between  4  and 
8  a/s  and  slightly  unstable  conditions.  The  equivalent  average 
CDB  is  expected  to  be  about  1.3x10-*  with  a  20%  variability. 
Accordingly,  Paulson  al.  find  10*CDH  =  1.32  for  19  runs. 

It  appears  as  if  the  aajor  source  of  the  disagreeaent 
between  the  data  of  figure  28  and  the  data  used  by  Garratt 
(1977)  is  that  the  latter  includes  aeasurenents  froa  onshore  and 
shallow  water  (less  than  15  aeters)  sites.  Equations  6.1  are 
consistent  with  aost  of  the  deep  water  (aore  than  50  aeters) 
drag  coefficients  to  within  possible  aeasureaent  error.  The 
discrepancies  suggest  a  slightly  higher  drag  coefficient  at  wind 
speeds  below  10  a/s.  In  contrast,  recent  aeasureaents  (10 
hours)  in  Bass  Strait  (Antonia  §£  al. .  1978)  in  winds  between  5 
and  10  a/s  and  -0. 1  <Z/L<  0  give  an  average  10*CD  at  5  aeters  of 
1.05  and  1.25  (about  0.9  and  1.1  at  10  aeters)  froa  the  Beynolds 
flux  and  dissipation  aethods,  respectively.  Since  a  constant 


COM  is  a  good  fit  to  most  individual  data  sets,  it  is  likely  to 
remain  a  good  description  of  the  overall  open  ocean  situation 
below  about  10  n/s.  The  "best"  constant  to  use  is  debatable, 
but  it  should  be  in  the  range  1.1x10“3  to  1.3x10-3.  The  average 
103C10*s  (neutral  or  near  .neutral)  from  all  the  discussed  open 
ocean;  eddy  correlation,  dissipation  and  profile  aeasureaents  in 
winds  less  than  10  to  12  m/s  are  about;  1.17  (178  runs),  1.14 
(626  hours)  and  1.24  (958  profiles,  188  hours),  respectively. 
The  overall  average  is  1.20,  but  only  1.17  if  the  results  from 
Hoeber  and  from  Brocks  and  Kriiggermeyer  are  weighted  by  1/6  and 
1/4,  respectively,  to  make  them  comparable  to  hourly  values.  At 
the  higher  wind  speeds,  the  vast  majority  of  the  open  sea  drag 
coefficients  come  from  Smith  (1979)  and  figure  28,  which  are 
both  well  described  by  6. 1.  it  is  an  interesting  feature  of 
6.1,  that  an  extrapolation  to  50  m/s  fits  the  hurricane  and  wind 
flume -data  compiled  by  Garratt  (1977)  as  well  as  a  continuation 
of  the  Charnock  line  of  figure  30. 

The  Variability  Of  CDN  And  7,0 

It  may  be  possible  to  describe  turbulence  and  the  drag 
coefficients  with  a  variable  Zo,  although  this  is  expected  to  be 
difficult  and  there  are  no  measurements  of  surface  conditions 
available.  The  wind  speed  dependence  of  CDN  above  lOm/s  shows 
that  a  constant  Zo  is  not  applicable.  A  Charnock  representation 
with«K  =0.0123  (K=0.40),  passes  within  ±1  standard  deviation  of 
all  the  band  means  of  figure  30,  but  its  predicted  low  wind 
speed  behavior  is  not  observed.  A  Charnock  line  with  «<  *0.008 
fits  the  region  from  010  *  11  to  17  m/s  fairly  well,  but  then 
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falls  off  such  too  quickly  to  account  for  the  values  up  to  26 
a/s.  The  CDH * s  to  50  a/s  in  Garratt  (1977) ,  require  «<  to  be 
about  0.016.  COM  is  very  nearly  a  constant,  0.00114,  for  4 
<01 0<  10  a/s,  as  is  the  average  of  CIO.  This  behavior  nay  be  a 
consequence  of  highly  variable  sea  surface  conditions  tending  to 
aake  Zo  increase  with  decreasing  wind  speed  and  balancing,  on 
average,  the  tendency  of  lover  u*  values  tc  sake  it  saaller. 
Above  10  a/s,  the  important  paraaeters  seea  to  change  character, 
so  that  Zo  varies  with  wind  speed,  producing  the  observed  near 
linear  rise  in  CDH  with  010.  It  aay  be  possible  to  describe  Zo 
by  a  Charnock  foraulation  with  an  experimentally  determined 
function  of  various  surface  paraaeters. 

To  explore  the  influence  of  the  sea  surface  on  the 
turbulence  above  it,  it  is  useful  to  examine  the  time  histories 
of  the  winds  and  fluxes.  Ar  example,  with  hourly  averages,  is 
run  D33C,  figure  31.  The  drag  coefficient  is  seen  to  vary  quite 
smoothly  with  time,  such  that  the  random  variability  about  a 
running  average  over  a  fev  hours  is  only  about  10*.  The 
increase  with  wind  speed  is  evident  from  hour  6  to  21.  In 
addition,  CDH  is  seen  to  behave  very  differently  on  the  rising 
wind  than  it  does  when  the  wind  speed  is  dropping.  The  falling 
wind  is  also  accompanied  by  an  almost  complete  reversal  in 
direction  and  it  is  this  situation  that  gives  rise  to  the  high 
drag  coefficients  at  low  wind  speeds.  The  picture  of  lower  than 
average  drag  coefficients  on  the  rising  wind  and  higher  ones  on 
the  falling  wind  is  also  supported  by  runs  D29D  and  D31F, 
figures  32  and  33.  Denman  and  Hiyake  (1973)  observed  that  drag 
coefficients  tended  to  increase  on  the  leading  side  of  a  storm. 
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then  either  renain  constant  or  decrease  slightly.  This  is  also 
a  very  apt  description  of  the  tine  series  in  figures  31,  32  and 
33.  The  scaling  arguaents  of  section  2.1  predict  that  the 
difference  between  the  10a  and  surface  aoaentua  flax  should  be 
greater  (giving  lower  coefficients) ,  on  the  falling  wind. 
Apparently  the  increase  in  Zo  in  this  situation  is  aore  than 
enough  to  overcone  this  loss.  It  is  possible  to  speculate  that 
Zo  depends  on  surface  paraaeters  that,  as  the  wind  speed  falls, 
reaain  near  their  "old"  values  produced  by  the  previously  higher 
winds.  If  °<  were  dependent  on  the  total  wave  generating  force 
and  hence  010-Cw  (section  2.3) ,  where  Cw  is  the  conponent  of  the 
wave  speed  in  the  i  direction,  it  would  be  very  sensitive  to 
sudden  changes  in  wind  direction,  because  these  would  tend  to 
aake  Cw  saaller  and  perhaps  even  negative.  It  appears  likely 
that  sone  of  the  observed  scatter  in  drag  coefficients  is  due  to 
the  influence  of  the  past  history  of  the  wind  as  "reneabered"  by 
the  sea  surface.  With  access  to  wave  spectra.  Denaan  and  Biyake 
(1973)  concluded  that  the  drag  coefficient  was  dependent  on  the 
nature  of  the  wave  field  to  the  order  of  20%. 

In  order  to  ezaaine  their  influence  on  the  average  CDS, 
tower  data  froa  various  wind  conditions  are  presented  in  table 
X.  Unfortunately,  above  15  a/s  there  are  too  few  data  for 
aeaningful  averages.  Winds  of  "constant  speed"  and  "constant 
direction"  are  allowed  to  have  4  n/s  and  30  degree  ranges 
respectively.  Data  used  in  the  "after  direction  changes" 
category  are  froa  the  six  hours  following  a  change  in  direction 
of  at  least  a  60  degrees  in  less  than  two  hours.  Peak  winds  are 
included  in  the  rising  winds.  The  suggestion  is,  that  either  a 
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|  WIND  SPEED  BANDS  (M/S)  1 


Wind 

conditions 
(nnnber  of  hoars) 


Beg cession 
and 

correlation 


4.5-7 


7-9 


9-11 


11-15 


All  tower  data 
(1086) 


79  +.043  010 
.59 


1.11 

(176) 


1.13 

(264) 


1.17 

(259) 


1.33 

(247) 


Constant  speed 
and  direction 
(253) 


,71  ♦. 049  010 
.60 


1.03 

(31) 


1.  13 
(68) 


1.14 

(79) 


1.31 

(59) 


Fising  winds 
(167) 


64  +.054  0  10 
.70 


1.01 

(22) 


1.11 

(33) 


1.14 

(31) 


1.25 

(49) 


Bising  wind 
steady  direction 
(81) 


60  +.055  010 
.67 


1.08 

(3) 


1.11 

(15) 


1.10 

(18) 


1.24 

(29) 


Falling  wind 
(111) 


.94  +.036  010 
.51 


1.40 

(3) 


1.25 

(12) 


Falling  wind 
steady  direction 
(80) 


1. 05+..025  010 
.36 


1.31 

(2) 


1.28 

(8) 


1.28 

(24) 


1.29 

(18) 


1.38 

(49) 


1.34 

(29) 


After  direction 
changes 
(60) 


1.18 

(22) 


1.34 

(12) 


1.61 

(16] 


TABLE  X  Mean  103CDN  of  wind  speed  bands  in  different  vind 
conditions.  Bracketted  values  are  the  naaber  of 
hoars  contained  in  an  average  or  regression  of  CDN 
against  010. 


direction  change  or  a  falling  wind  tends  to  produce  higher  than 
average  drag  coefficients.  Table  X  also  indicates  that  if 
observations  were  available  only  daring  rising  winds  with 
constant  direction,  as  nay  be  the  case  in  sone  data  sets,  the 
average  CDN  value  for  7  <010<  11  i/s,  would  be  about  4%  lower 
than  given  by  6.1.  The  observations  of  CDN  during  rising  winds 


show  a  systematic  rise  with  wind  speed  and  a  regression  of  CDH 
against  010  gives  10*  CDH*  0.64*  0.054  010  with  a  correlation 
coefficient  of  0.70.  Thus,  if  only  these  data  were  available 
one  could  easily  conclude  that  the  average  drag  coefficient 
increases  linearly  with  wind  speed  even  below  10  a/s.  Snch  a 
trend  nay  also  emerge  even  if  data  with  constant  speed  and 
direction  are  included.  Similar  effects  are  expected  at  the 
higher  speeds  too,  so  in  order  to  arrive  at  a  realistic  mean 
drag  coefficient  for  all  wind  speeds  it  is  imperative  that 
measurements  be  made  during  all  likely  wind  conditions.  It  is, 
therefore,  possible  that  drag  coefficient  formulations  nay  vary 
with  location  and  time  of  year,  with  the  break  from  a  constant 
to  a  linearly  increasing  CDH  being  quite  variable. 

*  The  tower  data  is  also  sorted  into  monthly  groupings  in 
table  XI.  There  is  no  trend  in  the  September  data  and  the 
average  CDH's  are  lower  than  those  of  the  other  months.  There 
is  a  reasonably  significant  trend  in  the  October  data  throughout 
the  entire  wind  speed  range.  The  high  CDH's  at  low  winds  are 
evident  in  the  December  results,  while  Barch  seems  tc  reflect 
average  conditions,  table  Z.  The  results  of  experiments 
conducted  at  different  times  of  the  year  may,  therefore,  be 
considerably  different,  possibly  because  wind  conditions  vary 
throughout  the  year.  It  would  be  desirable  to  include  data  from 
every  month  in  any  drag  coefficient  formulation.  Unfortunately, 
no  measurements  from  Hay  through  iugust  are  available  from 
either  the  tower  or  the  weathership,  but  the  high  wind  speed 
formulation  should  not  be  affected  as  these  are  not  the  months 
with  the  greatest  winds. 


WIND  SPEED  BANDS  (M/S) 


Hind 

conditions 
(number  of  hoars) 


Regression 

and 

correlation 


5-7 


7-9 


9-11 


11-15 


Septeaber 

(241) 


1.05+.0C1  010 
.008 


1.07 

(94) 


1.04 

(79) 


1.05 

(47) 


1.06 

(21) 


October 

(131) 


.83  *.047  010 
.65 


1.11 

(28) 


1.16 

(25) 


1.30 

(28) 


1.45 

(34) 


December 

(251) 


,78 


+.044  010 
.62 


1.22 

(20) 


1.14 

(42) 


1.14 

(53) 


1.33 

(73) 


March 

(271) 


77  +.044  010 
.65 


1.12 

(24) 


1.  18 
(49) 


1.  17 
(96) 


1.31 

(70) 


TABLE  XI  Mean  lO’CDN  of  wind  speed  bands  fcr  different 
aonths.  Bracketted  values  are  the  nuaber  of  hoars 
contained  in  an  average  or  regression  of  CDN 
against  010. 


The  Empirical  Drag  Coefficient 

The  investigation  of  the  tover  data  has  shown  that  neither 
stability  nor  fetch  (greater  than  10  ka)  have  a  significant 
effect  on  the  average  neutral  drag  coefficient.  It  has  been 
demonstrated  that  aach  of  the  variability  in  CDN  measurements, 
not  only  between  individual  data  points,  but  between  entire  data 
sets,  aay  be  due  to  the  influence  of  the  surface  wave  field. 
However,  it  appears  as  if  the  aeasureaent  of  the  important 
surface  parameters  could  possibly  be  more  complicated  than 
direct  moaentum  flux  measurements.  An  important  parameter  that 
can  be  practically  incorporated  into  a  neutral  drag  coefficient 
foraulation  is  the  wind  speed. 


I 


The  empirical  COM  formulation  of  6.1  gives  bulk  estimates 


of  the  momentum  flux  from  the  wind  speed,  UZ  and  stability, 
Z/L (AT) ,  following  the  procedure  outlined  in  section  2.3.  An 
iterative  technique  is  used  to  solve  equation  2.15  for  010.  For 
the  first  pass  010*=  0Z/  [1+  0.1  K(Z,Z/L)  ]  is  assumed  (COM*/2/  K 
-0.1),  giving  a  CON'  from  6.1.  on  subsequent  iterations  (never 
more  than  3)  010'  is  found  from  2.15  using  the  CON*  cf  the 
previous  pass  until  successive  010*  values  change  by  less  than 
1%.  CON'  from  the  final  pass  is  substituted  into  2.14  to  give  a 
CD,  which  gives  the  momentum  flux  through  2.9.  Bulk  estimates 
are  to  be  compared  to  direct  estimates  in  later  time  series  that 
include  data  with  -1.0  <Z/L<  0.3.  The  stability  need  not  be 
included  in  the  010  calculation,  because  the  maximum  stability 
of  the  ship  data  is  only  0.08  and  on  the  tower,  where 
stabilities  are  higher,  the  reduction  is  only  frcm  13  to  10 
meters,  therefore,  (section  2.3)  the  worst  error  is  less  than  2% 
in  010  and  about  1%  in  CD  and  <uw>.  The  error  in  Z/L  (is!) 
(section  4.1)  introduces  an  error  in  CD  of  about  5%.  A  2%  error 
in  0Z  becomes  4%  in  <uw>.  The  large  standard  deviations  in 
figure  30  and  the  range  of  previous  measurements,  suggest  that 
any  COM  formulation  has  at  least  a  10*  uncertainty.  The 
combined  error  in  bulk  estimates  from  equation  2.9  are, 
therefore,  greater  than  15%  as  are  the  errors  in  direct 
estimates.  Although  individual  bulk  estimates  may  differ  from 
direct  dissipation  or  Beynolds  flux  measurements  by  50%  or  more, 
there  should  be  reasonable  agreement  on  average,  albeit  the 
necessary  averaging  period  may  vary  with  time  and  place. 
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In  all  the  tiae  series  presented,  figures  31  to  34,  the 
■oeentna  flax  or  Be yn olds  stress,  (in  I/b*  oz  pascals.  Pa)  is 
calculated  froa  the  dissipation  aethod  (pluses) ,  froa  the 
Beynolds  flux  aethod  (squares)  and  froa  the  bulh  for aula  of  6. 1 
(solid  line).  Bun  0290,  figure  32,  is  a  typical  Septeaber 
situation.  During  the  unstable  conditions  of  the  rising  wind 
(to  hour  22)  the  bulk  foraula  gives  a  total  aoaentna  input  of 
2.23  H-hr/a*,  which  is  34%  higher  than  found  froa  the 
dissipation  calculations.  Upon  reaching  the  peak  winds  the  two 
aethods  are  in  auch  better  agreeaent  and  the  deviation  over  the 
whole  run  is  reduced  to  18%.  A  front  is  passing  the  tower  in 
D33C  (figure  31)  and  a  siailar  situation  occurs,  but  with  the 
rising  and  falling  wind  effects  balancing,  the  overall  deviation 
is  only  1%.  In  0290  the  falling  wind  and  changing  direction 
effects  are  not  great  enough  to  sake  the  average  stress  froa  the 
two  aethods  equal.  It  would  be  expected  that  these  influences 
are  strongly  felt  after  the  passage  of  low  pressure  systeas, 
because  as  centers  pass  the  winds  quickly  fall,  then  regain 
velocity  as  the  direction  changes  after  which  they  drop  off 
slowly.  i  low  pressure  center  is  passing  to  the  west  of  the 
tower,  froa  south  to  north,  in  D31F,  figure  33.  Unfortunately 
the  earlier  winds  froa  the  east  were  obstructed  by  other 
instruaentation  on  the  tower,  stability  is  always  very  nearly 
neutral  or  slightly  stable  and  should  not  have  been  a  factor  in 
causing  the  large  observed  stresses  that  give  a  total  aoaentua 
input  over  the  period  shown  of  6.5  B-hr/a*,  which  is  12%  higher 
than  given  by  the  bulk  calculations,  i  run  froa  Barch  1977, 


D38B,  figure  34  shows  that  all  three  aethods  can  be  in  excellent 
agreement  over  a  long  period  of  tise  provided  the  wind  is 
reasonably  steady.  Over  this  period  of  nearly  two  days,  the 
emulative  momentum  inpat  from  the  dissipation  calculations  is 
less  than  3%  higher  than  that  foand  from  the  balk  formula. 

In  general  the  fieynolds  flux  method  verifies  the 
dissipation  calculations,  but  an  exception  is  found  at  hour  28 
of  run  D29D,  figure  32.  The  flux  run,  T61,  is  an  average  of 
only  3  groups  (40  minutes)  starting  IS  minutes  after  the 
beginning  of  the  hour  28  dissipation  average.  In  this  15  minute 
interval  there  was  a  sudden  4  a/s  increase  in  wind  speed,  which 
the  hourly  averaging  smooths.  The  simultaneous  dissipation  run 
does  agree  with  T61  (Appendix) ,  illustrating  that  identical  time 
intervals  are  required  for  intercoaparisons.  The  flux  run. 
Til 2,  at  hour  14  of  run  D33C  (figure  31)  is  also  found  to  be  in 
much  better  agreement  with  its  simultaneous  dissipation  run  than 
with  the  hourly  average.  Bun  D33C  discloses  a  sampling  problem 
associated  with  the  tower  operation  of  the  Reynolds  flux  system. 
Host  of  the  low  wind  speed  runs  were  collected  in  September  when 
the  drag  coefficients  tend  to  be  low.  Later,  when  the  more 
variable  winds  seem  to  produce  some  very  high  coefficients  at 
wind  speeds  below  10  a/s,  the  Beynolds  flux  system  was  prevented 
from  recording  by  the  wind  speed  limit  setting.  Thus  the 
Beynolds  flux  data  set  (figure  19)  tends  to  be  biased  toward 
snail  CO*  s,  relative  to  6.1,  at  the  low  wind  speeds  and  to 
exhibit  a  trend  with  wind  speed  over  the  whole  range  of 
sea sure sen ts.  Above  the  anxious  wind  speed  limit  setting  (12 
a/a)  6.1  fits  the  data  of  figure  19  very  well. 


TIME  (HOURS) 


FIGURE  33  Boaentaa  flax  tiae  series  frci  ran  D31F.  Starting 
tiae  is  4:00  GUI  October  20,  1976.  The  solid  line 
represents  the  balk  estiaates  of T. 
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It  appears  that  the  bulk  aerodynamic  method  gives  good 
estimates  of  the  total  momentum  input  or  average  stress  over 
periods  of  a  few  days  or  more.  The  bulk  estimate  should  improve 
if  measured  or  approximated  stability  conditions  are  included. 
It  has  been  shown  that  in  steady  winds  the  hourly  bulk  estimates 
are  good  measures  of  the  stress,  but  with  more  variable  winds 
they  may  consistently  differ,  over  periods  up  to  a  day,  from 
dissipation  estimates.  With  comparable  error  in  both  methods, 
it  is  not  obvious  which  is  the  more  accurate.  It  is  conceivable 
that  the  varying  winds  are  producing  systematic  errors  in  one  or 
more  of  the  dissipation  method's  assumptions  and  causing  the 
discrepancies.  However,  if  this  occurred  it  is  likely  that  the 
overall  scatter  of  figure  28  would  be  greater  than  the  scatter 
in  drag  coefficients  calculated  from  other  methods,  which  is  not 
observed.  It  is  felt  that,  since  at  least  some  of  the  scatter 
in  measured  drag  coefficients  at  the  same  wind  speed  is  real, 
dissipation  measurements,  although  subject  to  random  errors, 
ought  to  follow  changes  in  the  real  stress  more  closely  than 
bulk  estimates.  A  bulk  formula  based  on  long  term  averages  nay 
not  be  strictly  applicable  to  short  tern  phenoiena,  such  as  wave 
development  during  a  rising  wind  or  the  deepening  of  an  inlet's 
upper  layer  following  a  switch  from  down  to  strong  up  inlet 
winds.  In  such  cases,  it  may  be  possii le  to  devise  an 
appropriate  formulation  cabable  of  providing  good  hourly  stress 
values  with  minimal  measurement  error.  A  possible  means  of 
incorporating  variable  wind  effects  would  be  to  allow  the  drag 
coefficient  to  depend  on  the  past  history  of  the  wind. 


150 


6.3  Balk  Aerodynaaic  Parameterization  Of  The  Sensible  Beat  Flax 

's.. 

The  sensible  heat  flax  is  paraaeterized  in  terns  cf  a 

surface  -  air  teaperatnre  difference.  The  aean  air  teaperatare 

at  10a,  T10,  is  obtained  froa  the  average  teaperatare  at  the 

aeasareaent  height,  TZ,  through  equation  2.11.  At  the  Bedford 

tower  the  surface  teaperatare,  TSFC,  is  approxinated  by  a  sea 

teaperature,  TSSA,  aeasured  about  10a  below  the  aean  sea  level, 

where  the  heating  and  cooling  of  the  surface  by  radiation  and 

heat  exchanges  with  the  ataosphere  aay  not  be  followed  exactly. 

On  CC6S  Quadra  a  surface  "bucket"  teaperatare  was  recorded  as 

part  of  the  3-hourly  meteorological  observations.  These  are 

interpolated  to  give  hourly  TSSA  values,  which  aay  be  a  problem, 
• 

because  the  ship  was  steaming  through  large  sea  surface 
teaperatare  gradients  during  the  one  run  that  useful  teaperature 
aeasureaents  were  recorded.  At  any  time  then,  TSSA  may  differ 
substantially  from  a  representative  surface  teaperature.  As  a 
precaution  that  the  Stanton  nuaber,  CT  (equation  2.9) ,  is  not 
too  adversely  affected,  data  used  in  the  paraaeterization  are 
restricted  to  conditions  with  |TSEA-TZ|  >  1.0°C,  even  though  the 
flux  aeasurements  aay  still  be  reliable. 

Figure  35  is  a  plot  of  <wt>  against  010  -AT  (AT*TSEA-T10) , 
for  all  129  hours  of  unstable  teaperature  data,  and  a  regression 
line  (the  average  of  <wt>  against  010  AT  and  010  AT  against 
<vt>) .  The  23  hours  of  data  froa  the  weathership  (plotted  as 
pluses)  show  considerable  scatter,  but  no  systeaatic  departure 
froa  the  tower  results  (triangles) .  The  regression  is 
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FIGURE  35  Paraaeterizaticn  of  the  sensible  heat  flax  in  °Cb/s 
in  unstable  stratification.  Triangles  represent 
tover  data  and  pluses  are  frci  CCGS  Quadra. 


<wt>  *  0.00100  010  AT  ♦  0.0029  °Cm/s 
103  CT1 0  »  1.00  ♦  2.9  (°Ca/s)  /  (010  AT)  ,  (6.2) 

with  a  correlation  coefficient  of  0.86.  This  is  almost 
identical  to  the  foranla  given  by  Friehe  and  Schmitt,  1976,  for 
0<  010  AT  <25  °Ca/s  only.  Bestricting  the  parameterization  to 
data  with  |<wt>|  >0.004°Cm/s,  010  >4.0  a/s  and  |aT|  >1°C 
eliminates  low  heat  flnz  and  small  010  AT  situations,  bat  a 
positive  heat  flax  is  still  predicted  at  010  AT  »  0,  although 
the  value  of  0.002  °Ca/s  given  by  ?riehe  and  Schmitt,  who  used 
mostly  small  heat  flax  data,  is  probably  more  realistic.  The 
additional  data  at  higher  sensible  heat  flaxes  indicate  that  a 
single  parameterization  is  acceptable  from  010  at  =0  to  possibly 
more  than  100  °Cm/s.  Again  this  result  is  in  accord  with  the 
BIO  tower  data  of  Smith,  1979,  bat  the  average  coefficient  is 
aach  lower  than  found  by  Smith  and  Banke  (1975)  from  a  limited 
number  of  Sable  Island  measurements. 

The  sensible  heat  flux  time  series  from  run  D33C,  figure 
36,  shows  that  seemingly  very  small  Stanton  numbers  arise  during 
a  rapid  increase  in  air  temperature,  just  after  AT  first  goes 
negative.  Fortunately,  Beynolds  flux  measurements  (squares)  are 
also  available  over  this  period  and  their  good  agreement  with 
the  dissipation  method  verifies  the  treatment  of  equation  2.21. 
Temperature  spectra  from  the  Beynolds  flux  recordings  do  not 
display  any  evidence  of  salt  contamination,  so  it  is  felt  that 
the  <wt>  values  from  this  period  are  reasonably  accurate. 
However,  it  is  possible  that  they  are  significantly  less  than 
their  surface  values.  The  loss  of  sensible  heat  flux  with 
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PIGOBB  36  Use  series  of  the  sensible  beat  flux  fros  run  D33C 
(figure  31).  line  is  fros  4:40  GHT  December  7# 
1976.  Hs  is  fros  the  dissipation  (pluses)#  Reynolds 
flux  (squares)  and  bulk  (solid  line)  setbods. 
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height  was  discussed  in  section  2.1  and  if  the  0. 5°C/hour 
heating  were  due  solely  to  the  vertical  flux  divergence,  then  a 
40%  loss  in  flux  at  13b  would  be  expected.  Qualitatively  the 
low  CTN  values  observed  can  be  explained  by  speculating  that 
<wt>  was  measured  above  the  "constant  sensible  heat  flux"  layer, 
since  if  this  were  the  case,  the  percentage  of  the  loss  should 
decrease  as  the  magnitude  of  the  surface  flux  increases  and  no 
loss  should  be  evident  after  the  peak  temperature.  However, 
there  is  probably  a  great  deal  of  heating  due  to  horizontal 
advection,  which  would  mean  less  loss  of  flux  with  height. 
There  are  other  situations  of  rising  temperature  which  do  not 
seem  to  be  affected  to  such  a  great  degree  so  possibly  the 
measured  <wt>« s  do  not  differ  from  the  surface  flux  by  as  much 
as  the  low  CTN's  would  suggest.  knot her  possibility  is  that 
TSEA  at  a  depth  of  about  10m  is  not  tracking  the  surface 
temperature.  In  the  12  hours  prior  to  run  D33C  there  was  very 
little,  wind  and  the  air  temperature  was  below  -15°C,  giving  a 
huge  AT  of  about  20°C.  It  appears  as  if  the  surface 
temperature,  TSFC,  may  have  become  much  colder  than  TSEA  during 
this  time,  because  as  the  air  warmed  and  the  wind  increased, 
presumably  mixing  the  unstable  water  column,  TSEA  actually 
decreased  from  3.9  to  3.5°C.  It  is  not  impossible  that  the 
surface  heating  produced  by  the  later  high  air  temperatures  was 
prevented  from  mixing  down  to  the  sea  temperature  probe  by  the 
then  stable  upper  water  column.  In  this  situation  TSEA  could 
become  less  than  TSFC  by  enough  to  produce  the  observed  behavior 
of  033C  and  to  upset  the  parameterization  of  the  sensible  heat 
flux.  There  was  only  one  other  similar  situation  at  the  tower 
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and  even  lower  Stanton  nnabers  were  calculated,  but  these  data 
were  rejected  on  the  basis  of  abnorsal  temperature  spectra. 
Despite  the  preceding  arguments,  lower  than  average  CT  values 
■ay  occur  in  the  circumstances  described,  but  they  are  not 
likely  to  occur  over  the  open  sea,  so  this  period  of  D33C  is  not 
included  in  the  parameterization  and  analysis  of  the  stable 
sensible  heat  flux. 

Figure  37  is  a  plot  of  -<wt>  against  -010  (TSEA-T10)  from 
131  hours  of  stable  stratification  with  the  suspect  D33C  results 
plotted  as  crosses.  The  regressions  of  -<wt>  against  -0 10  AT 
and  -010  AT  against  -<wt>  for  the  123  triangles  have  a 
correlation  coefficient  of  0.93  and  when  averaged  yield 

<wt>  ■  0.00075  010  AT  ♦  0.0020  °Cm/s 
103  CT10  *  0.75  ♦  2.0  (°Cm/s)  /  (010  oT)  ,  (6.3) 

which  is  also  plotted  on  figure  37.  igain  a  positive  heat  flux 
of  about  0.002  ®Cm/s  at  0  AT*0  is  indicated  even  though  there  is 
no  data  with  |0  AT|<  5°Cm/s.  Inclusion  of  the  suspect  D33C  data 
lessens  the  slope  to  about  0.00065.  Friehe  and  Schmitt  only 
considered  data  with  0  AT>  -15°Cm/s  and  their  suggested  formula 
differs  from  6.3  by  only  10%  at  U  AT  *  -15°Cm/s.  There  are  only 
a  few  data  points  beyond  0  AT*  -35°Cm/s  that  contribute  to  6.3, 
which,  therefore,  cannot  be  expected  to  be  representative  of 
■ore  negative  0  AT*s.  The  BIO  tower  results  (Smith,  1979) 
contain  more,  large  negative  heat  flux  runs  and  suggest  a  larger 
CT10  of  about  0.00083.  Although  the  BIO  regression  differs  from 
6.3  by  about  13%  at  0  AT  *  -100  °Cm/s,  the  data  points  overlap 


Paraae terization  of  the  sensible  heat  flax  (°Ca/s) 
is  stable  stratification.  Crosses  are  froa  D33C. 
Solid  line  is  a  regression  of  the  triangles  only. 
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very  veil.  &  single  parameterization  of  the  sensible  heat  flax 
in  stable  stratification  should  adequately  describe  all  the 
scasureaents,  but  it  appears  to  differ  froa  the  unstable  case. 

It  is  felt  that  the  uncertainty  in  010  <oT  is  perhaps 
greater  than  in  <wt> ,  so  band  averaging  is  carried  out  over 
bands  of  <vt>  and  the  results  are  tabulated  in  table  XII.  The 
band  averages  are  reasonably  well  described  by  6.2  and  6.3  over 
the  entire  range  of  measurement,  thus  linear  fits  to  the  data  of 
figures  35  and  37  are  appropriate.  The  ratio  of  the  averages 
gives  a  CT10  for  each  band  and  the  unstable  values  are  clearly 
higher  than  the  stable  ones.  The  individual  points  are  far  too 
fev  and  scattered  for  any  trends  to  be  seen.  Por  example,  the 
103CT10  value  of  1.21  in  the  0.05  to  0.07  °Cs/s  band,  beccaes 
1.04  if  010  ^.T  is  increased  by  only  1/2  a  standard  deviation. 

StabU&Y  Niseis 

The  stability  dependence  of  the  neutral  Stanton  nuaber  is 
exaained  in  figure  38,  vhere  CTN ,  equation  2.12,  is  plotted 
against  Z/L.  The  solid  lines  represent  the  overall  averages: 

123  stable  103  CIS  «  0.69 

129  unstable  10>  CTH  -  1.08  ,  (6.4) 

about  which  the  standard  deviations  are  0.16  and  0.36, 
respectively.  There  is  not  a  great  deal  of  data  on  which  to 
base  definitive  conclusions,  but  it  appears  as  if  CTH  and  Zot 
are  reasonably  independent  of  Z/L  in  stable  and  unstable 
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TABLE  XII  Parameterization  of  the  sensible  heat  flux  by  band 
averaging  U  AT  over  ranges  of  <vt>. 

stratification  separately.  However,  the  discontinuity  in  the 
mean  CTN  at  Z/L«0,  implies  a  dramatic  change  in  Zot.  This 
feature  could  be  incorporated  into  the  theory  by  relating  Zot  to 
a  parameter  that  changes  character  abruptly  at  neutral 
stability,  but  that  otherwise  has  very  little  stability 
dependence  as  is  indicated  by  the  relative  constancy  of  the 


average  CTN  away  from  Z/L=0. 

There  is  a  great  deal  of  scatter  in  figure  38  and  it  should 
be  reiterated  that  nicrobead  contanination  by  salt  particles 
coaid  only  be  checked  at  the  few,  sonetines  infrequent,  tines 
for  which  spectra  were  available  fron  the  Reynolds  flux  systen. 
It  is  therefore  possible  that  sone  contaninated  data  has  escaped 
notice  and  is  responsible  for  sone  of  the  high  CTN  values  at 
Z/L<0.  Similarly,  because  tenperature  and  humidity  are 
oppositely  correlated  in  stable  stratification,  contaminated 
data  nay  be  responsible  for  some  low  CTN's  at  Z/I  >0.  The  high 
CT 10  values  following  hour  40  of  run  D38B  (figure  39)  are  guite 
likely  a  result  of  bead  contanination,  because  a  Reynolds  flux 
run  at  about  hour  44  was  rejected  on  the  basis  of  a  lack  of  low 
freguency  variance  in  the  temperature  spectrum  and  the  run  at 
hour  41  was  a  borderline  case. 

Bind  Speed  Effects 

Equation  2.12  relates  CTN  to  Zo,  indicating  that  it  should 
increase  with  wind  speed  above  10  n/s,  if  Zot  does  not 
counteract  the  increase  in  Zo.  The  situation  may  be  complicated 
by  the  observed  dependency  of  Zot  on  Z/L,  which  correlates  with 
wind  speed.  The  effect  of  wind  speed  on  CTN  is  investigated  for 
the  stable  and  unstable  case  in  tables  XIII  and  XIV 
respectively.  In  the  stable  case  there  is  no  indication  of  an 
increase  in  CTN  with  wind  speed.  In  fact  above  10  n/s  where  CDN 
begins  to  increase,  there  is  a  hint,  albeit  not  significant, 
that  CTN  decreases.  In  table  XI?  the  neans  are  scattered,  but 
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TABLE  XIII  Averaged  neutral  Stanton  nmber  as  a  function  of 
wind  speed  in  stable  stratification  only. 

again  there  see ns  to  be  no  obvious  wind  speed  dependency. 
However,  nearly  all  of  the  data  in  table  XIV  above  18  a/s  cone 
froa  CCGS  Quadra,,  which  did  not  really  yield  enough  data  at  the 
lower  wind  speeds  to  fully  interconpare  the  ship  and  tower 
results.  There  is  a  suggestion  in  table  XIV  that  the  unstable 
CTH  increases  froa  DIO  «  10  to  18  a/s.  Larger  CTH  values  are 
reported  over  shallow  water  by  Erancey  and  Garratt,  1978,  who 
find  an  increase  with  wind  speed  given  by  10»CTH  *  0.083  010  ♦ 
0.48.  The  contention  that  CTH  increases  with  010  is  apparently 
supported  by  run  033C,  figure  36,  where  the  increase  in  0Z  froa 
10  to  13  a/s  (hours  7  to  10)  is  accoapanied  by  a  sharp  rise  in 
CTH  froa  about  1.0  to  1.4x10-*.  The  rise  in  CTH  with  winds 
above  10  a/s  appears  to  be  greatly  perturbed  by  the  transition 
to  stable  stratification  and  although  the  following  stable  data 
are  suspect,  the  rise  seeas  to  subsequently  continue.  Ho 
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TABLE  XIV  Averaged  neutral  Stanton  number  as  a  function  of 
wind  speed  in  unstable  stratification  only. 


definite  conclusion  is  possible  because  of  the  lacX  of  data  and 
the  strong  influence  of  stability.  For  example,  if  the  suspect 
D33C  data  were  included  in  the  14  to  18.2  m/s  range  of  table 
XIII,  the  band  average  would  decrease  giving  a  smaller  CTN  at 
the  highest  speed,  even  though  these  additional  CTN's  increase 
with  wind  speed. 


These  results  do  not,  therefore,  rule  out  the  possibility 
that  CTN  follows  the  wind  speed  dependency  of  Zo.  However, 
because  this  trend  is  not  supported  by  the  Quadra  results, 
because  of  the  limited  amount  of  data,  some  of  which  is  possibly 
contaminated  (run  D38B) ,  and  because  the  stability  effects  are 
not  fully  understood,  the  neutral  Stanton  number  is  not 


formulated  as  a  function  of  wind  speed.  It  is  an  interesting 
observation  that,  on  average,  CTN  for  Z/L<0,  0.00108,  is  very 
nearly  equal  to  CDN  for  4. 0<  010  <10  i/s,  0.00114,  iaplying  Z 6s 
Zot*  0.0058  ca.  It  is  possible  that  Zot  remains  constant  above 
10  a/s  vith  CTN  increasing  due  to  Zo.  If  this  is  trne,  CTN 
conld  be  obtained  froa.an  experimentally  determined  Zot  and  a  Zo 
given  from  a  CON  formulation,  using  equation  2.12. 

Method  Comparison 

The  bulk  estimates  of  the  sensible  heat  flux  shown  by  the 
solid  lines  of  figures  36  and  39,  are  calculated  from  the 
stability  dependent  CTN  given  by  6.4.  As  outlined  in  section 
2.3,  the  neutral  Stanton  number  is  converted  to  a  CT  at  the 
measurement  height,  wind  speed  and  stability.  The  error  in  bulk 
estimates  found  from  <wt>  «  ct  OZ  AT,  equation  2.9,  is  about  2% 
from  OZ,  perhaps  10X  from  CTN  and  only  5X  from  Z/L,  because  the 
bulk  stability  estimate  is  always  available.  In  addition 
considerable  error  is  introduced  through  AT,  which  may  sometimes 
be  in  error  by  0.5°C.  Even  vith  a  large  AT  of  5°C,  the  error  is 
10*,  making  the  bulk  sensible  heat  flux  calculations  more 
uncertain  than  those  of  the  momentum  flux.  iith  a  large 
possible  error  and  problems  with  the  data,  it  is  difficult  to 
compare  the  bulk  and  dissipation  estimates  in  the  time  series. 
The  Beynolds  flux  calculations  (squares)  lend  credibility  to  the 
dissipation  estimates  (pluses) ,  because  of  the  generally 
excellent  agreement  between  the  two  methods.  Figure  39  shows 


the  sensible  heat  flux  time  series  from  run  D38B  and  (like  the 
momentum  flux  of  figure  34)  the  bulk,  eddy  correlation  and  the 
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FIGURE  39  Sensible  heat  flux  tine  series  fron  run  D38B 
(figure  34).  Tiae  is  froa  16sC0  GUI  Harch  12,  1977 
and  symbols  are  the  sane  as  in  figure  36. 
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dissipation  estiaates  of  the  sensible  heat  flax  continue  to 
agree  for  a  considerable  length  of  tine.  On  the  whole  the 
sensible  heat  flax  parameterization  is  conpatible  with  the 
results  of  Friehe  and  Schoitt,  1976,  and  Saith,  1979,  so  a 
combination  of  all  available  data  should  give  a  CTN  formulation 
capable  of  providing  reasonable  averages  (over  a  few  days  or 
aore)  of  sensible  heat  flux  values,  depending  on  the  accuracy  of 
at.  However,  as  with  the  momentum  flux,  there  nay  be  short 
lived  situations,  as  possibly  seen  in  D33C,  where  a 
parameterization,  valid  for  long  tern  averages,  does  not 
strictly  apply.  It  is  also  possible  that  erronous  aT  values  are 
causing  the  difference  between  the  bulk  and  dissipation 
calculations  between  hours  2  and  10  of  D33C,  figure  36  and,  as 
has  been  discussed,  between  hours  12  and  24. 


CHAPTER  7 


The  experimental  program  described  in  this  thesis 
successfully  measured  the  momentum  and  sensible  heat  fluxes  over 
the  sea  at  winds  between  4  and  26  m/s.  As  hoped,  a  great  many 
hours  of  momentum  flux  data  were  suitable  for  analysis,  but  much 
less  sensible  heat  flux  and  no  moisture  flux  data  were  found  to 
be  reliable.  Success  depended  chiefly  on  the  performance  of  the 
sensors  and  on  the  establishment  of  the  dissipation  method  as  a 
viable  means  of  measuring  the  fluxes  of  momentum  and  sensible 
heat. 


The  velocity  sensor  worked  very  well,  tut  there  were  some 
problems  with  the  temperature  measurements  and  no  humidity 
sensor  was  found  to  be  suitable  for  remote  operation  in  a 
salt-air  environment.  The  Gill  twin  propeller-vane  anemometer 
operated  for  periods  of  more  than  a  month  in  adverse  conditions 
without  servicing.  It  provided  the  fluctuating  horizontal  and 
vertical  velocities  to  the  Reynolds  flux  system  and  responded  to 
the  lower  frequencies  of  the  downstream  velocity  spectrum's  -5/3 
region  sufficiently  for  the  molecular  dissipation  to  be 
inferred,  although  the  propeller  responses  first  had  to  be 
determined.  The  distance  constant  was  found  to  depend  on  the 
type  and  weight  of  propeller,  the  wind  speed  and  the  angle  of 
attack.  The  humidity  sensitivity  of  salt  contaminated 
microbeads  was  recognized  by  the  lack  of  lew  freguency  variance 
in  the  temperature  spectrum,  but  only  sometimes  were  other 
characteristics,  such  as  the  absence  of  a  -5/3  region,  observed. 


167 


This  behavior  regains  a  major  problen  vith  the  renote  operation 
of  this  type  of  sensor.  The  response  of  the  aicrobeads  seesed 
to  be  limited  by  its  protective  enclosure  and  was  described  by  a 
distance  constant  of  about  0. 90m,  which  is  adequate  for  both  the 
Beynolds  flux  and  dissipation  nethods. 

Beynolds  flux  aeasurenents  from  the  Bedford  tower  were 
shown  to  be  realistic  by  conparisons  with  spectra,  cospectra  and 
turbulence  statistics  fron  previous  studies.  In  addition,  the 
drag  coefficients  and  Stanton  nunbers  were  generally  comparable 
to  the  results  of  Smith,  1979,  however  there  is  probably  a  bias 
to  small  drag  coefficients  at  low  wind  speeds  as  a  result  of  the 
Beynolds  flux  sampling.  Universal  shapes  for  the  velocity 
spectra  and  cospectra  in  both  the  stable  and  unstable  cases, 
were  found  from  averages  over  the  196  momentum  flux  runs.  Hith 
only  60  temperature  runs  available,  there  was  considerable 
uncertainty  in  the  normalized  temperature  spectrum  and  w,t 
cospectrum.  The  integration  of  all  cospectra  began  at  n*0.004, 
then  the  unstable  and  stable  fuw  {f )  and  the  unstable  and  stable* 
♦wt  (f )  integrals  were  multiplied  by  1.06,  1.005,  1.10  and  1.04, 
respectively,  in  order  to  account  for  the  lower  frequencies. 
This  method  was  found  to  preserve  covariance,  on  average,  and  to 
reduce  the  scatter  in  eddy  correlation  measurements  caused  by 
the  uncertain  low  frequency  contributions  to  the  fluxes. 
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Reliable  Reynolds  flux  estimates  were  needed  for  comparison 
with  simultaneous  dissipation  calculations  and  the  dissipation 
nethod  was  shown  to  give  essentially  the  same  results,  on 
average.  The  agreement  between  the  two  methods  was  found  to  be 
best  when  the  magnitude  of  the  fluxes  was  large.  In  all  but  the 
most  stable  stratification  (Z/L>0.05),  u*DISS1 ,  the  neutral 
dissipation  method,  which  does  not  require  an  explicit  stability 
parameter,  was  found  to  be  in  quite  good  agreement  (to  within 
about  20%)  with  eddy  correlation  values  of  u*.  The  agreement 
between  the  two  techniques  improved,  particularly  the  sensible 
heat  flux  calculations,  when  the  stability  modification  of  the 
logarithmic  profiles  and  the  buoyant  production  were 
incorporated  into  the  dissipation  method  (u*DISS2  and 
<wt>DlSS2) .  These  corrections  involved  the  stability  parameter, 
for  which  a  bulk  estimate  Z/L(^T)  was  shown  to  be  a  reasonable 
approximation,  on  average.  A  linear  regression  gave 

u*DISS  =  0.96  u*FL0X  ♦  0.025  m/s  , 

where  the  positive  offset  results  from  a  tendency  for  <uw>DISS 
to  be  greater  than  <uw>F10X  by  more  than  30%,  in  the  most  stable 
runs  (Z/L>0. 10) .  In  near  neutral  conditions  (-0.45  <Z/L<  0.05) 
at  all  wind  speeds  the  momentum  flux  calculations  agreed  to 
within  4%,  on  average.  The  agreement  between  the  sensible  heat 
flux  calculations  was  very  good,  with  a  regression  giving 
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The  Bedford  tower  experiment  established  that  reliable 
dissipation  estimates  of  both  the  momentum  and  sensible  heat 
flnxes  and  the  balk  estimates  of  Z/L  could  be  obtained  from  the 
CCGS  Quadra.  k  favourable  comparison  of  ship  and  tower 
dissipation  drag  coefficients  showed  the  Bedford  tower  to  be 
essentially  an  open  ocean  site,  which  allowed  the  combined  1591 
hoars  of  momentum  flax  and  260  hours  of  sensible  heat  flux 
measurements  to  be  considered  as  a  single  open  ocean  data  set. 

The  dissipation  data  showed  the  neutral  drag  coefficient  to 
depend  on  wind  speed,  as  approximated  by 

10aCDN  ■  1.1*  4  <  010  <  10  m/s 

1 03CDN  »  0.49  ♦  0.065  010  10  <  010  <  26  m/s- 

Below  10  m/s  the  variability  of  CON  with  wind  speed,  fetch,  and 
stability  was  minimal,  <  5%  on  average.  Time  series  of 

calculated  CDN*s  displayed  about  a  10%  randou  fluctuation.  Bulk 
estimates  of  the  momentum  flux  were  obtained  from  the  above 
formulation  by:  first,  shifting  the  measured  wind  speed  to  10m 
with  an  iterative  technique  involving  CON;  second,  calculating 
CON;  third,  shifting  CON  to  CD  (the  drag  coefficient  at  the 
measurement  height,  Z,  wind  speed,  0Z,  and  stability  Z/L (AT)  ) ; 
fourth,  applying  the  bulk  aerodynamic  formula  -<uw>  *  CD  0Z*. 
It  was  shown  that  the  neglect  of  Z/L  would  produce  a  minor  error 
in  the  calculated  010  and  CON,  but  that  the  CO  found  in  step  3 
could  be  affected  by  as  much  as  20%.  However  this  error  was 
reduced  to  5%  by  using  the  bulk  stability  estimate.  Errors  in 
the  CON  formulation,  0Z  and  Z/L  could  add  up  to  15%  in  the  bulk 
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estimates,  which  were  seen  to  give  a  good  aeasnre  of  the 
aoaentum  flax  averaged  over  a  fev  days  or  sore  and  a  good  hourly 
average  when  the  wind  was  steady.  Over  periods  of  ap  to  a  day, 
the  bulk  and  dissipation  calculations  were  seen  to  consistently 
differ  by  as  ouch  as  30%.  These  discrepancies  were  found  to  be 
associated  with  varying  winds,  with  the  dissipation  estimate 
being  smaller  on  the  rising  wind  and  larger  on  the  falling  wind 
or  after  a  change  in  wind  direction.  The  conclusion  drawn  was 
that  the  surface  roughness  and  hence  drag  coefficient  depend  on 
surface  paraaeters  which  are  a  product  of  both  past  and  present 
winds. 


A  constant  neutral  10a  drag  coefficient  was  found  to  be  an 
adequate  description  of  alaost  all  recent  aeasureaents  over  deep 
water,  throughout  the  wind  speed  range,  4  to  10  a/s.  The  value 
of  the  constant  varied  from  about  1.1x10-*  to  1.3x10-*.  This 
behavior  has  been  found  by  the  eddy  correlation,  dissipation  and 
profile  nethods,  but  aeasureaents  froa  onshore  or  shallow  water 
sites  often  displayed  a  distinctly  different  wind  speed 
dependency.  A  reasonable  coaproaise  over  the  open  ocean  would 
be  to  use  10*CDN  -  1. 2  in  the  bulk  aerodynaaic  method  for  winds 
up  to  11  a/s.  At  U10  x  11  a/s,  the  high  wind  speed  regression 
of  this  study,  10*CDN  =  0.49  ♦  0.065  010  =  1.2  .  Since  this 
regression  fits  the  only  other  large  set  of  open  sea  high  wind 
speed  data  (Smith,  1979),  it  should  give  satisfactory  CDH's  at 
the  higher  wind  speeds  (to  at  least  26  a/s) . 
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It  was  necessary  to  paraaeterize  the  sensible  heat  flax 
differently  in  stable  and  unstable  stratification.  The  majority 
of  the  data  were  in  the  ranges  -40  <010  at<  60°Cb/s,  -0.3  <2/L< 
0.1  and  6  <0 10<  18  a/s,  for  which  either 

103<wt>  *  1.00  010  at  ♦  2.9  °Cn/s  stable 
10*<wt>  *  0.75  010  aT  ♦  2.0  °Cn/s  unstable 

or  103CTH  ®  1.08  unstable 
103CTH  *  0.68  ,  stable 

described  the  results,  (lore  aeasureaents  of  CTN  were  needed  in 
order  to  confirs  or  deny  a  wind  speed  dependency,  which  was 
suggested  by  soae  of  the  data.  Bulk  estiaates  of  the  sensible 
heat  flux  were  obtained  with  <wt>  *  CT  02  AT,  where  CT  was  found 
froa  Z/L,  COB,  CD  and  the  stability  dependent  CTH.  The  errors 
were  coaparable  to  those  of  the  bulk  aoaentua  flux  estiaates, 
except  for  the  uncertainty  in  AT,  which  could  have  added 
considerable  error  depending  on  its  aagnitude.  It  was  noted 
that  all  available  data  are  reasonably  consistent,  aaking  a  good 
paraaeterization  of  the  sensible  heat  flux  feasible. 
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0.321 

0.324 

0.294 

T68 

18 

06 

29/  9 

11.0 

232 

12.8 

15.5 

-.11 

0.370 

0.352 

T69 

19 

48 

29/  9 

10.2 

226 

13.3 

14.3 

-.06 

0.356 

0.321 

T7C 

23 

00 

29/  9 

11.8 

222 

13.9 

13.7 

-.01 

0.392 

0.410 

T71 

2 

00 

30/  9 

13.6 

236 

14.  1 

13.2 

0.01 

0.420 

0.437 

T72 

5 

00 

30/  9 

10.  1 

243 

13.4 

12.8 

0.01 

0.352 

0.329 

T73 

11 

30 

30/  9 

8.  3 

239 

12.4 

13.1 

-.07 

0.270 

0.256 

0.059 

0.075 

T74 

2 

23 

10/10 

13.9 

194 

18.8 

16.2 

0.04 

0.523 

0.496 

-.276 

-.246 

T7  5 

5 

23 

10/10 

17.  1 

191 

19.4 

16.3 

0.04 

0.638 

0.670 

-.368 

-.374 

T76 

8 

23 

10/10 

12.7 

224 

16.5 

16.2 

-.00 

0.476 

0.470 

-.055 

-.064 

T77 

11 

23 

10/10 

15.8 

223 

17.1 

16.2 

0.01 

0.660 

0.580 

-.043 

-.067 

T78 

12 

53 

10/10 

17.5 

232 

16.6 

16.1 

0.00 

0.714 

0.682 

T79 

14 

01 

10/10 

17.8 

234 

15.6 

16.0 

-.01 

0.762 

0.676 

T80 

17 

23 

10/10 

19.2 

238 

13.8 

15.  1 

-.02 

0.786 

0.755 

T81 

20 

23 

10/10 

16.6 

254 

12.2 

14.  2 

-.04 

0.641 

0.624 

T82 

22 

23 

10/10 

14.0 

268 

11.5 

13.6 

-.05 

0.523 

0.504 

T83 

23 

23 

10/10 

13.6 

271 

11.  1 

13.3 

-.06 

0.442 

0.495 

T84 

2 

23 

11/10 

12.3 

265 

10.1 

12.4 

-.07 

0.396 

0.417 

T85 

5 

23 

11/10 

10.5 

269 

9.2 

12.5 

-.  14 

0.401 

0.373 

T86 

16 

01 

14/10 

7.7 

177 

13.3 

7.7 

0.30 

0.  173 

0.203 

-.027 

-.052 

T87 

17 

17 

14/10 

10.  2 

174 

13.6 

8.3 

0.18 

0.318 

0.297 

-.112 

-.105 

T88 

18 

35 

14/10 

11.4 

179 

13.2 

10.2 

0.08 

0.355 

0.376 

-.162 

-.  158 

T89 

20 

17 

14/10 

6.3 

242 

12.  1 

10.7 

0.09 

0.224 

0.237 

190 

8 

23 

15/10 

16.  1 

259 

7.6 

10.5 

-.05 

0.625 

T9 1 

11 

23 

15/10 

15.7 

263 

7.9 

10.1 

-.04 

0.692 

T92 

12 

29 

15/10 

14.4 

263 

7.7 

9.7 

-.05 

0.515 

T9  3 

14 

23 

15/10 

14.6 

259 

8.2 

9.3 

-.03 

0.559 

T9  4 

17 

23 

15/10 

13.7 

259 

10.5 

9.0 

0.02 

0.546 

0.491 

T95 

19 

00 

15/10 

13.2 

260 

11.4 

8.8 

0.04 

0.535 

0.496 

-.186 

-.187 

T96 

2 

10 

27/11 

12.4 

199 

7.9 

5.7 

0.04 

0.496 

0.441 

-.  181 

-.  168 

T97 

2 

10 

28/11 

11.3 

215 

8.1 

5.4 

0.07 

0.377 

0.416 

-.173 

-.229 

T98 

0 

00 

- i 

3/12 

11.7 

123 

- 1 

5.4 

4.6 

i 

0.02 

- 1 

0.393 

0.411 

-.052 

-.059 
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1  ■■  1 

— — 

1 

- 

BOR 

TIBS 

DAIS 

02  °TBOE  TZ 

IS  PC 

Z/L 

a*  (i/s) 

10<wt> 

GMT 

1976 

FLUX 

DISS 

FLUX 

DISS 

T9  9 

6 

00 

3/12' 

12.4 

216  6.9 

4.6* 

0.05 

0.424 

0.447' 

-.180 

-.185' 

T100 

9 

00 

3/12 

17.9 

250  3.1 

4.6 

-.02 

0.826 

0.752 

T101 

12 

00 

3/12 

18.6 

264  -2.5 

4.6 

-.09 

0.733 

0.737 

T102 

15 

00 

3/12 

18.5 

274  -9.7 

4.6 

-.19 

0.746 

0.724 

T103 

18 

00 

3/12 

16.9 

278  -11. 

4.6 

-.26 

0.672 

0.690 

T 104 

21 

00 

3/12 

15.0 

270  -12. 

4.5 

-.33 

0.552 

0.588 

T10S 

0 

00 

4/12 

16.8 

267  -11. 

4.5 

-.25 

0.735 

0.690 

T106 

3 

00 

4/12 

19.3 

269  -9.6 

4.5 

-.17 

0.817 

0.806 

T107 

6 

00 

4/12 

15.9 

266  -10. 

4.4 

-.27 

0.633 

0.599 

T108 

9 

00 

4/12 

16.0 

261  -8.9 

4.3 

-.22 

0.551 

0.559 

T109 

12 

00 

4/12 

14.4 

264  -7.8 

4.3 

-.26 

0.500 

0.495 

T110 

15 

00 

4/12 

10.7 

261  -5.8 

4.2 

-.34 

0.352 

0.414 

Till 

15 

00 

7/12 

11.7 

169  2.5 

3.5 

-.04 

0.418 

0.380 

0.179 

0.200 

T112 

18 

00 

7/12 

14.6 

139  4.1 

3.5 

0.01 

0.478 

0.519 

T113 

21 

00 

7/12 

17.  1 

143  5.8 

3.6 

0.03 

0.660 

0.651 

-.131 

-.125 

Til  4 

0 

00 

8/12 

18.6 

153  7.6 

3.5 

0.04 

0.746 

0.736 

-.385 

-.352 

T115 

3 

00 

8/12 

18.9 

159  7.8 

3.6 

0.04 

0.728 

0.729 

-.435 

-.516 

T116 

6 

00 

8/12 

13.3 

170  7.9 

3.6 

0.08 

0.518 

0.506 

-.372 

-.452 

T117 

18 

00 

9/12 

14.6 

300  -11. 

3.6 

-.31 

0.601 

0.598 

T118 

21 

00 

9/12 

18.3 

302  -14. 

3.5 

-.27 

0.640 

0.739 

T119 

0 

00 

10/12 

17.  1 

302  -17. 

3.4 

-.33 

0.675 

0.695 

T120 

3 

00 

10/12 

15.6 

292  -17. 

3.4 

-.38 

0.656 

0.607 

T121 

6 

00 

10/12 

16.9 

289  -15. 

3.3 

-.31 

0.487 

0.592 

T122 

9 

00 

10/12 

13.5 

278  -14. 

3.3 

-.45 

0.463 

0.485 

T123 

2 

00 

11/12 

10.5 

211  5.2 

3.3 

0.06 

0.362 

0.348 

-.069 

-.089 

T124 

8 

00 

11/12 

12.4 

220  5.9 

3.4 

0.06 

0.483 

0.462 

-.159 

-.152 

T125 

11 

00 

11/12 

15.3 

230  6.3 

3.5 

0.04 

0.661 

0.578 

-.228 

-.176 

T 128 

14 

00 

11/12 

16.0 

231  6.5 

3.5 

0.04 

0.606 

0.613 

-.251 

-.186 

T127 

17 

00 

11/12 

17.2 

291  1.8 

3.5 

-.03 

0.645 

0.673 

T128 

20 

CO 

11/12 

16.0 

290  -1.0 

3.5 

-.08 

0.565 

0.591 

T129 

23 

00 

11/12 

11.9 

308  -2.6 

3.5 

-.21 

0.378 

0.433 

T130 

14 

00 

12/12 

12.6 

201  1.6 

3.4 

-.06 

0.446 

0.405 

0.265 

0.296 

T131 

5 

00 

14/12 

14.2 

317  -14. 

3.2 

-.35 

0.531 

0.541 

T132 

8 

00 

14/12 

14.6 

312  -16. 

3.2 

-.43 

0.534 

0.543 

T133 

23 

00 

14/12 

13.2 

218  -7.0 

3.0 

-.27 

0.493 

0.456 

T 134 

2 

00 

15/12 

18.2 

220  2.2 

3.0 

-.01 

0.769 

0.735 

T135 

5 

00 

15/12 

17.5 

219  3.9 

2.9 

.0.01 

0.693 

0.681 

-.023 

-.047 

T136 

8 

00 

15/12 

15.9 

213  4.8 

2.9 

0.03 

0.663 

0.602 

-.132 

-.082 

T137 

11 

00 

15/12 

1977 

15.6 

231  5.4 

2.9 

0.04 

0.613 

0.592 

T138 

0:48 

13/  3 

10.6 

233  1.9 

0.0 

0.06 

0.221 

0.264 

-.097 

-.146 

T139 

9 

48 

17/  3 

10.4 

236  1.8 

0.2 

0.05 

0.363 

0.355 

-.142 

-.  104 

T140 

12:48 

17/  3 

11.2 

240  1.6 

0.3 

0.03 

0.395 

0.387 

-.071 

-.097 

T 14 1 

15:  48 

17/  3 

11.0 

234  1.5 

0.2 

0.04 

0.356 

0.359 

-.079 

-.061 

T142 

18:48 

17/  3 

12.6 

259  1.0 

0.3 

0.01 

0.456 

0.456 

T143 

21:48 

17/  3 

12.8 

270  -0.3 

0.3 

-.02 

0.440 

0.464 

T144 

0 

:48 

18/  3 

12.6 

278  -1.3 

0.2 

-.04 

0.442 

0.457 

T145 

3:48 

18/  3 

10.5 

280  -1.9 

0.2 

-.09 

0.360 

0.387 

0.214 

0.253 

T146 

6:48 

18/  3 

11.3 

271  -2.4 

0.2 

-.09 

0.389 

0.371 

0.268 

0.237 

T147 

9:48 

18/  3 

11.7 

270  -2.2 

0.2 

-.07 

0.380 

0.410 

0.227 

0.279 

- H 

180 


BO  B 

time' 

GHT 

DATE 

1977 

UZ  ®TBOEl 

TZ 

TSFC' 

z/l' 

a*  (a/s)  10<at> 

FLUX  DISS  FLUX  DISS 

_ L 

T148 

12 

V 

48 

18/ 

— r 

3 

11.5 

282 

-1.7 

o.  2* 

-.06’ 

0.405 

- 1 - r 

0.421 

T149 

15 

48 

18/ 

3 

12-7 

254 

0.4 

0.3 

-.00 

0.398 

0.444 

T150 

6 

48 

19/ 

3 

17.3 

61 

-1.0 

0.3 

-.02 

0.694 

0.685 

T151 

9 

48 

19/ 

3 

12.8 

16 

-1.7 

0.2 

-.05 

0.458 

C.488 

T152 

12 

48 

19/ 

3 

16.4 

347 

-2.8 

0.2 

-.05 

0.646 

0.645 

T153 

15 

48 

19/ 

3 

17.6 

325 

-2.9 

0.3 

-.05 

0.685 

0.708 

T154 

18 

48 

19/ 

3 

14.  1 

323 

-2.8 

0.3 

-.08 

0.469 

0.530 

T155 

0 

48 

20/ 

3 

10.7 

296 

-1.9 

0.2 

-.08 

0.347 

0.405 

T156 

21 

48 

24/ 

3 

12.3 

358 

-0.7 

0.1 

-.03 

0.369 

0.425 

T157 

0 

48 

25/ 

3 

12.2 

356 

-1.1 

-0.1 

-.03 

0.359 

0.415 

T158 

3 

48 

25/ 

3 

11.3 

357 

-0.9 

-0.2 

-.03 

0.321 

0.369 

T159 

9 

48 

25/ 

3 

10.7 

353 

-0.5 

-0.2 

-.02 

0.321 

0.360 

1160 

12 

48 

25/ 

3 

10.3 

3 

0.3 

-0.2 

0.01 

0.290 

0.329 

T161 

15 

48 

25/ 

3 

10.6 

359 

1.5 

-0.2 

0.05 

0.310 

0.342 

T162 

18 

48 

25/ 

3 

10.5 

358 

2.2 

0.0 

0.07 

0.276 

0.341 

T163 

21 

48 

25/ 

3 

10.4 

356 

0.9 

0.0 

0.03 

0.312 

0.368 

T164 

3 

48 

26/ 

3 

10.7 

353 

0.1 

0.0 

-.00 

0.  303 

0.352 

T165 

6 

48 

26/ 

3 

10.5 

350 

0.1 

0.0 

0.00 

0.259 

0.333 

T166 

9 

48 

26/ 

3 

10.3 

353 

0.5 

0.0 

0.01 

0.267 

0.333 

T 167 

12 

48 

26/ 

3 

11.4 

352 

1.3 

0.0 

0.04 

0.328 

0.369 

T168 

15 

48 

26/ 

3 

10.6 

351 

0.9 

0.0 

0.02 

0.323 

0.365 

T 169 

18 

48 

26/ 

3 

12.2 

17 

1.4 

0.0 

0.03 

0.388 

0.424 

T170 

21 

48 

26/ 

3 

12.3 

14 

1.0 

0.0 

0.02 

0.395 

0.421 

T171 

3 

48 

27/ 

3 

10.8 

10 

1.4 

0.0 

0.04 

0.341 

0.371 

T172 

12 

48 

27/ 

3 

11.3 

4 

1.6 

0.0 

0.04 

0.315 

0.370 

T173 

15 

48 

27/ 

3 

10.4 

17 

2.4 

0.0 

0.07 

0.313 

0.350 

T174 

18 

48 

27/ 

3 

10.9 

3 

3.5 

0.0 

0.10 

0.266 

0.335 

T175 

21 

48 

27/ 

3 

8.0 

1 

2.7 

0.0 

0.13 

0.235 

0.287 

T176 

21 

48 

29/ 

3 

8.6 

172 

2.2 

0.0 

0.10 

0.232 

0.298 

T  177 

18 

48 

31/ 

3 

11.5 

168 

3.0 

0.0 

0.08 

0.367 

0.413 

T178 

21 

48 

31/ 

3 

7.8 

204 

3.  1 

0.0 

0.16 

0.200 

0.242 

T17S 

12 

48 

3/ 

4 

12.6 

147 

0.6 

0.0 

0.01 

0.451 

0.428 

T180 

15 

48 

3/ 

4 

14.2 

125 

1.2 

0.0 

0.02 

0.562 

0.531 

T181 

18 

48 

3/ 

4 

11.4 

170 

3.1 

0.0 

0.09 

0.407 

0.445 

T182 

3 

48 

4/ 

4 

19.  1 

288 

1.8 

0.0 

0.02 

0.626 

0.762 

T183 

6 

48 

4/ 

4 

18.0 

300 

-0.3 

0.0 

-.01 

0.684 

0.723 

T184 

9 

48 

4/ 

4 

16.2 

313 

-1.4 

0.0 

-.03 

0.610 

0.628 

T185 

12 

48 

4/ 

4 

15.2 

311 

-0.7 

0.0 

-.01 

0.511 

0.566 

T186 

15 

48 

4/ 

4 

13.5 

315 

0.4 

0.0 

0.01 

0.462 

0.519 

T187 

18 

48 

4/ 

4 

10.8 

321 

1.7 

0.0 

0.05 

0.311 

0.388 

T188 

21 

48 

4/ 

4 

10.9 

311 

1.8 

0.0 

0.05 

0.316 

0.374 

T189 

3 

48 

6/ 

4 

15.  1 

189 

6.0 

0.0 

0.09 

0.566 

0.573 

T190 

6 

48 

6/ 

4 

13.0 

209 

4.7 

0.0 

0.11 

0.482 

0.471 

T191 

9 

48 

6/ 

4 

10.7 

223 

1.3 

0.0 

0.04 

0.346 

0.399 

T192 

12 

48 

6/ 

4 

12.4 

224 

1.3 

0.0 

0.03 

0.466 

0.460 

T193 

15 

48 

6/ 

4 

12.4 

220 

1.2 

0.0 

0.02 

0.460 

0.448 

T194 

18 

48 

6/ 

4 

14.5 

223 

1.9 

0.0 

0.03 

0.545 

0.540 

T195 

21 

48 

6/ 

4 

13.2 

223 

1.5 

0.0 

0.03 

0.461 

0.494 

T196 

0 

48 

7/ 

4 

11.6 

217 

0.9 

0.0 

0.02 

0.406 

0.427 

